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Further calculations are reported on the problem of the distribution in size of cosmic-ray 
showers. In §2 it is shown that this distribution is completely determined when the average 
energy distribution of the particles in a shower is known. In §3 previous calculations of the 
fluctuation in size of showers have been revised and extended. The main result is that for the 
simplified model chosen (the so-called Furry model) the fluctuations are small and roughly 
equal to twice the Poisson value for all values of the thickness (see Table II; Fig. 3). In §4 
another simplified model is considered for which it is possible to take the ionization exactly 
into account. In §5 and §6 the calculations of the fluctuation are extended to the actual cosmic- 
ray problem. For one value of the initial energy and for one depth a numerical calculation has 
been made (see Table III). The result for the fluctuation is again a few times the Poisson value. 


§1. INTRODUCTION 


HE problem of the distribution in size of 

cosmic-ray showers has been treated by 
several authors but no satisfactory solution has 
yet been given. The question is to determine the 
probability P(Eo, N, x) that N particles emerge 
from a layer of matter of thickness x, when an 
electron of energy Epo falls normally upon it. 
Furry! succeeded in solving the problem for a 
special model in which the essential approxima- 
tion consisted of the neglect of the ionization. 
He found: 


P(N, x)=(1/N)[1—(1/N)}", (1) 


where N is the average number of particles, 

which in this case is equal to exp (Bx) and is 

independent of Eo. The quantity B is a material 
* Now at Amherst College. 


1W. H. Furry, Phys. Rev. 52, 569 (1937). See also B. 
Rossi and K. Greisen, Rev. Mod. Phys. 13, 240 (1941). 


constant. On the other hand Bhabha and Heitler? 
have asserted that the function P(N, x) will be 
essentially the Poisson distribution : 


P(N, x) =e-"(N)*/N! (2) 
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Fic. 1. The Furry and Poisson distributions 
(1) and (2) for N=10. 


? H. J. Bhabha and W. Heitler, Proc. Roy. Soc. 159, 432 
(1937). Compare also: H. Euler, Zeits. f. Physik 110, 450 
(1938), N. Arley, Proc. Roy. Soc. 168, 519 (1938). 
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This may be obtained by assuming no correla- 
tion between the different particles in a shower. 
The striking difference between (1) and (2) can 
be seen from Fig. 1. In the actual cosmic-ray 
case one has to take into account both the corre- 
lation between different particles and the ioniza- 
tion. It seems likely that the result will lie some- 
where between (1) and (2). 

In order to estimate which one of the two 
results is more nearly correct, especially when one 
varies Ey and x, an attempt was made in a 
previous paper* to calculate the fluctuation 
o =(N*)—(N)x? more exactly. The same simpli- 
fied model, as chosen by Furry, was considered, 
except that the ionization was now taken ap- 
proximately into account by the so-called cut-off 
method (see I, p. 350). The results were found to 
lie between the values which follow from (1) 


and (2), namely: 


(N?)y—(N) x2 = N(N — 1) (1a) 
for the Furry distribution‘ and: 
(NN?) — (N) a? = N (2a) 


for the Poisson distribution. It seemed to us 
necessary to refine and to extend these calcula- 
tions, since no estimates of error were made in I 
and since the dependence of the fluctuation on 
the thickness x remained unclear. The results are 
given in §3; in contrast to the values found in I, 
p. 358, we now find that the fluctuation is much 
smaller than the Furry value (1a), and that it is 
roughly twice the Poisson value (2a) for the 
interesting range of values of x (see Fig. 3). 
The reason for this discrepancy was traced to the 
fact that in I the calculation of the integrals by 
the method of steepest descent was not made 
with sufficient accuracy. The actual behavior of 
the fluctuation as a function of x is quite curious 
and is difficult to explain in a qualitative way. 
In §4 we have therefore considered another 
simplified model, which has much less similarity 


3 A. Nordsieck, W. E. Lamb, Jr., and G. E. Uhlenbeck, 
Physica 7, 344 (1940). In the following this paper will be 
quoted as I and we shall use the same notations as much as 
possible. 

‘ In the comparison the value of N is taken which follows 
from the multiplication curve with the ionization taken into 
account, instead of the value exp (Bx). This procedure was 
first proposed by C. G. and D. D. Montgomery [Phys. 
Rev. 53, 955 (1938) ]. It is hard to give a logical justification 
for this, but it is the best which one can do. 


G. E. URLERSBSECK 

with the cosmic-ray problem, but for which the 
influence of the ionization can be taken exactly 
into account. For this model one can show easily 
that without ionization the fluctuation is N —1, or 
practically the Poisson value. However, here 
again the influence of the ionization on the 
fluctuation (which can be computed exactly) is 
strange, and difficult to explain qualitatively. 

In §5 and §6 we have extended the calculations 
to the actual cosmic-ray equations. Again we 
have taken the ionization into account by means 
of the cut-off method. As was to be expected, 
the results for the fluctuation are quite similar to 
those for the Furry model. 

Before going into the details of these calcula- 
tions we shall show in §2 that the function 
P(Eo, N,x) is completely determined if one 
knows the average energy distribution F(o, E, x) 
of the particles in a shower. For the Furry model 
we shall write down a formal expression for the 
connection between these two functions. In the 
case of no ionization this leads again to (1) with 
N =exp (Bx). We have, however, been unable to 
use this connection when the ionization is taken 
into account, so that the problem of the actual 
shape of the distribution function P(N, x) re- 
mains unsolved. 


§2. A FORMAL EXPRESSION FOR P(N, x) 


As explained in I §3, all statistical questions 
regarding the formation of showers can be 
answered if one knows the probability of a given 
energy distribution of the particles which emerge 
after the thickness x. Let us assume first that 
the energy of any particle can only have the 
discrete values €1, €2, €3---. An energy distribu- 
tion is then specified by giving the numbers 
Ni, N2, N3--- which have these energy values. 
The function we seek is W(m, m2---;x), the 
probability that at depth x we have the distribu- 
tion n;. We shall call it the ‘“‘master function.” 
Every other statistical function may be obtained 
from the master function by taking suitable 
average values. For instance, the average number 
of particles of energy ¢; is given by: 


n(x) = Sn;W(n, neo: - (3) 


4% 


where the round summation signs will always 
mean a sum over all possible values that each n; 
may have. Analogously one can form the quad- 
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ratic averages: 
nin ;(x) = Snn;W(n,, me? x), (4) 


and so on. For a continuous energy variation 
n(x) will become the average energy distribution 
F(E,x)dE; analogously (nm, will become a 
function of two energy variables, which we de- 
note by F2(/,, Es, x)dE,dE2, and so on. 

The function P(N, x) is obtained from W by 
the operation: 


P(N, x) =S'W(my, no-++3 x), (5) 


where the primed summation sign S’ means that 
one has to sum over all values of m; with the 
restriction : 

Di nz=N. 
P(N, x) is completely determined when one knows 
all the moments (NV*(x))4, and one may write: 


P(N, x) => ay.(N*(x))w, (6) 


where the ay; are numerical coefficients, which 
clearly are independent of the distribution func- 
tion P(N, x). Since on the other hand: 


N(x)=Di tix), (N*(x))w=L (nin(x))m, (7) 


and so on, it is clear that in this way P(N, x) 
will be related to the energy distributions F(E, x), 
F.(Fk4, Eo, x), etc. 

To discuss this connection further, we will 
consider first the case of no correlation. When the 
probabilities of finding a particle in the different 
energy intervals are completely independent of 
each other, the master function must have the 
form: 


(7i1)"(vig)™- 75 . 
W(m, m2, ---)=A , (8) 
n!n!- ie 





since 7; will then be proportional to the proba- 
bility of finding a particle with the energy €¢;; 
the. normalization constant A must be deter- 
mined from: 


» (n;)™ 
SW=A [I > rr N exp (}-’;) =1, 
+ n=0 n;! 


so that A=exp (—N). From the definition (5) 
one then shows easily that P(N, x) becomes the 


5 This follows by considering (6) as the solution of the 
equations 


(N*(x))y= 2 N*P(N,x); k=1,2,3--- 
N=1 


Poisson distribution (2). The quadratic averages 
(4) and also the higher order averages can in this 
case all be expressed in terms of the 7;. One finds 
for instance: 


(nn j)nv = nNj+ 6; jN; 
(nn jn) =n iN Ny + 6; jn k + 5 54 n jn (9) 
+ Oxnih Rng t+ 5556 jun F 


and the generalization is clear. Consequently one 
can express (N*), as a polynomial of degree k 


in N: 
(N*)w =X bar(N)!. (10) 


Of the coefficients },; we will only need the 
property : 
a (—1)- 
° 2, Queby: = ————__, (11) 
k=! Ni(l—N)! 
which follows immediately by introducing (10) 
in (6) and remembering that P(N, x) is now the 
Poisson distribution (2). 

Consider now the actual case, in which the 
different energy intervals will be correlated be- 
cause of the splitting processes. When a particle 
“splits” it adds simultaneously to the number in 
each of two energy intervals and hence the 
probabilities for particles to be in these ranges 
are not independent. A general method would 
now be to derive from the continuity equation, 
which the function W has to fulfill [the so-called 
‘“‘master’’ equation, see I, Eq. (22) ], equations for 
all the average values 7i;, (nj), etc. However, 
since these equations are almost self-evident, we 
shall omit the formal derivations and in addition 
pass immediately to the limit of a continuous 
energy variation. One finds that the quadratic 
averages F,(E,, Es, x) and the higher order aver- 
ages F,(E,, E2---E.,x) are singular whenever 
two or more of the energy variables coincide. 
Analogously to (9) one then can decompose the 
functions F; into regular parts as follows: 


F(E, x)=K,(E, x), 
F.(E,, Es, x) = K,(Fi, Eo, x) 
+ 6(£,— E2)Ki(Fi, x), 
F;(E,, Es, Es, x) = K3(F,, Eo, Es, x) (12) 
+ 6(E.1— Ex) K2(Fi, Es, x) 
+6(E.— E3;)Ke(Ee, Ei, x) 
+6(E;—E;)K2(Es, Es, x) 
+ 6(E,— E.)6(E.— E3)Ki(Ei, x), 
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and so on. The functions K;,(E,---E.,x) are 
regular and symmetric in all the energy variables. 
Since the decomposition rules (12) are exactly 
the same as (9), except for replacing the products 
of average values 7; by the functions K;, one 
obtains: 


(N*(x))w= bxi:M (x), (13) 


where 
Mix) = [+++ f dBy dE, K(By- +E x) (14) 


and the coefficients },; are the same as in (10). 
Introducing (13) in (6) and using the relation 
(11) one gets: 


2 (—1)-" 
Mt bat — Bt ded. 5 
(N, #)= 2 n-th ) 


The average energy distribution F(E, x) 
=K,(E, x) fulfills a continuity equation of the 


form: 
OF /dx=LeF(E, x), (16) 


where Lg is a linear operator acting on the 
variable EZ. Equation (16) is an abbreviation for: 


OF E 
—=-F(E, ») f q(E, u)du 
0 


a OF 
+2 f duq(u, E)F(u, x)+8—, (16a) 
E OE 


where g(H, u)du is the probability per unit 
thickness that a particle of energy E splits into 
two particles of which one has the energy be- 
tween u and u+du, while the other has the re- 
maining energy E—u; q(E, u) =q(E, E—u) since 
one does not distinguish between the two par- 
ticles. Finally 8 is the average energy loss per 
unit thickness due to: ionization.* For the higher 
order averages K,(E,--+E,, x) one then finds in- 
homogeneous equations of the form: 


0K ,/dx= (Le; +Leo+ s¢ -Le,)K, 
+I],(E,---Ei, x). (17) 


The inhomogeneous part J; consists of a sum of 
l(t—1)/2 terms of which a typical one is: 


2q(E:+ Es, E,)Ki-i(F-it+Eo, E;- ° -E, x). (18) 
6 See I, §2; the following considerations would remain 
valid for an arbitrary ionization probability p(E, «). 


The other terms correspond to the other ways in 
which one can get the / energy variables occurring 
in K; by one splitting process from (/—1) energy 
variables occurring in K;_:. The interpretation of 
(17) is clear : for J;=0 Eq. (17) would be separable 
and the solution would be F(A, x) F(Es, x)--- 
X F(Ei, x); Eqs. (12) would therefore become 
identical with (9), and one would get the case of 
no correlation. The correlation is therefore due 
to the inhomogeneous part J;, and one may say 
that J; is the probability per unit thickness that / 
particles with energies E;, E2---E,; are produced 
from (/J—1) particles by one splitting process. 
Since one particle of energy Ep is falling in, one 
has to solve Eqs. (17) with the initial condition 


F,(E,- - - Ey, 0) =6(E,:— Eo) - - -6(E:— Ep). 


From (12) one sees that this means that F(E, 0) 
=6(E—E,) while all K, for />1 are zero for 
x=0. Writing for the solution of (16) with this 
initial condition F(Eo, E, x) it is clear that the 
product F(&:, £1, x)---F(é:, Zi, x) is not only 
a solution of the homogeneous part of (17) 
but that it is also the so-called fundamental 
solution. One then verifies easily that with 
K,(E,---E:, 0)=0 the solution of the inhomo- 
geneous Eq. (17) is given by: 


Ki=f arf vf dé,-- -d&,F(&, Bi, x—t)--> 
0 0 0 


X Fé, E,, x—t)Ii(é1- + - &1, t). (19) 


This is a recurrence relation for K,; it becomes 
physically obvious when one considers the physi- 
cal meaning of the K;. From (19), (14), and (15) 
one sees therefore that P(N,x) is completely 
determined when the average energy distribution 
F(Eo, E, x) is known. Unfortunately, these equa- 
tions seem to be of little use for the practical 
determination of P(N, x). One can still write: 


Mi(x)= f ef dé dtif dtl (y+ -+€:, 0) 
0 0 0 
XN (Ei, x—t)---N(E1,x—2), (20) 


but only for the case of no ionization (@=0) and 
“homogeneous” splitting probability [q¢(Z, «) 
= x(u/E)/E; see I, §2] can one obtain from this 
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a recurrence relation for the M;. In this case: 


N(E,, x) = @Bz 


= f dex. 


Substituting in (20) one finds: 


with 


M(x) =I(1+1)Be'® f dte—'B' M,_,(t). 
0 


Since M(x) = N(x) one then finds by induction: 


M(x) =le'42(1 —e-87)"-!, 


Substitution in (15) leads to the Furry for- 
mula (1). 


§3. CALCULATION OF THE FLUCTUATION 
FOR THE FURRY MODEL 


As in I we shall take g(E, u)=1/E (which 
makes B=1) and we shall take the ionization 
into account by means of the cut-off method. 


One finds for N and (N?)y the expressions (see I 
Eqs. (18) and (38)): 


= ds 
N(z, »)=— f —ep [ss+2x/(s+1)], (21) 


21 Ss 


(N2(z, x) w=N(z, »)+2f déI(z, x, &), 


ee? I'(s) P(t) 
I Z, X, rao, 
sii (2x)? sf fet T(s+t+2) 


2& 
xexp | +os+ —-—— 
st+t+1 


(22) 








1 1 
+iHe- i -———+- ~)|. 
s+1 t+1 
Here z=log (Eo/e), where € is the cut-off energy ;’ 
all the integrals must be taken along paths 
parallel to the imaginary axis and to the right of 


all singularities of the integrand. They can be 
computed approximately with the saddle point 


7 The connection between ¢ and 8 can be fixed by means 
of the relation: 


N(x)dx = Eo/8, 
/0 
which follows strictly from (16a) and which is also physi- 


cally obvious. Introducing (21) and carrying out the 
integrals one finds: 2e*—1= E»/B so that «28. 





TABLE I. N(x) for the Furry model with z=4.75. 


N N 





So x main term %e correction corrected 

1/z=0.21 0.00 1.084 —8.3 0. 996 

0.3 1.20 2.89 

0.4 2.20 5.57 

0.5 3.09 8.68 

0.6 3.95 11.93 —1.1 11.80 

0.7 4.80 15.1 

0.8 5.67 17.9 —1.7 17.6 

0.9 6.57 20.0 

1.0 7.50 21.15 —1.8 20.8 

1.1 8.48 21.25 

1.2 9.47 20.4 

io 10.54 18.5 —1.6 18.2 

1.4 11.6 16.1 

15 12.8 13.3 —1.5 13.1 

1.6 13.9 10.4 

1.7 15.2 7.96 —1.5 7.84 

1.8 16.4 5.84 

1.9 17.8 3.97 

2.0 19.1 2.65 —1.3 2.62 





| 
| 


method. The following improvements were made 
over the analogous calculations reported in I: 

All variable factors in the integrand were 
written into the exponent and included in the 
saddle point expansion. Or in other words: no 
part of the integrand was considered as “slowly 
varying.” 

2. The equations which determine the position 
of the saddle points for given values of x [in (21) ] 
or of x and & [in (22)] now become quite in- 
volved. To circumvent this difficulty we have 
turned the question around; convenient values 
of the saddle points were chosen and the corre- 
sponding values of x or of were then determined 
from the equations. 

3. In order to estimate the errors involved in 
the saddle point method we calculated for a few 
values of x a next approximation by extending 
the Taylor expansion around the saddle points up 
to the fourth-order terms. 

For N one obtains: 


exp (—x+ ¢(5o)) 





N= ~ 
[2x¢’ (so)? 
IV(so)  Sy’""(So)? 
ns sa | (23) 
8y'"(so)? 249""(so)* 
where 
2x 
g(s) =sz+ -—log s (23a) 
s+1 


and sois the saddle point determined by ¢’(so) =0. 
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Fic. 2. The integrand J(z, x, £) of Eq. (22) as function 
of & for z=4.75 and for the four values of x indicated. 
For the calculation Eq. (24) is used; the full lines are the 
result of taking only the first term, while for the dotted 
lines the other terms given are also used. 


In Table I one finds the results for z=4.75. The 
correction given in the fourth column is com- 
puted from the terms with the third and fourth 
derivatives in (23). One sees that by using the 
main term only the error is already less than 10 
percent over the whole range of x, while with 
the correction it is even less than 0.5 percent. 
For I(z, x, £) one obtains: 


exp [E—2x+yY(5o, So) J 








I(z, x, §)= 
anv —¥) 
8888 4 ssst 3 sstt ssss”  Wssit 
xfs! ne a ie 
16(Yo+¥e)? —-8(¥2 —¥*) 
Wesss ‘a West + SWestt = 5 (Woes + 3West) . 





16(Wss— Yer)? 48(Weet Wer)? 


_ Poon + 2WessW sat pene 3Wet 
‘8 (Wes + st) _— v* ) 
= 3 (Wess — Wet)? _ 
1 6(y? — vy ) (Wss ies Pot) 








+f (24) 
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where 


2é 
¥(s, t) = (s+#)z+-—— 
sti 





+1 
1 1 r(s)P(t) 
+2(x—&) —_+——)+ og —— (24a) 
s+1 t+1 I'(s+t+2) 


and the subscripts denote differentiations after s 
and ¢. Since y(s, ¢) is symmetric in s and ¢, the 
saddle points so and ¢y for the s and ¢ integrations 
will be equal to each other; they are determined 
by ws(So, to) =W2(so, to) =0. The terms with the 
third- and fourth-order differential quotients are 
again correction terms and they have been sim- 
plified by making use of the symmetry of (s, ¢). 
For a given value of x, one now chooses a suitable 
set of values for the saddle points so and one 
then computes the corresponding values of & from 
Ws(So, So) =0 and of I(z, x, £) from (24). Figure 2 
shows the integrand J(£) for a few values of x 
and for z=4.75; the full lines are the result of 
taking only the main term of (24) into account, 
while for the dotted lines the complete expression 
(24) is used. The integration over must be done 
graphically, and in Table II one finds the results. 
One sees that the influence of the correction 
terms in (24) amounts to roughly 3 percent for 
all values of x. Table II also gives the values of 


the fluctuation ¢=(N?)y—(N)w divided by N, 
which are plotted in Fig. 3. The accuracy of 
these results depends mainly on the accuracy of 
the values for (N*),. We believe that with the 
correction terms in (24), the error in (N*), is at 


. 
most 1.5 percent. Since the error in (N)w is 
certainly much smaller, the error in ¢@ will be at 
most of the order of 25 percent. 


TABLE II. (N24, for the Furry model with x =4.75. 


o 


M =2| 1d¢(N?)m=, NN, 
Jo N+M (N2)qy—(N)*ay 





main main ———— % corr. (N2)ay a/N 
x term term N to M corrected corrected 
0 0 1 0 
2.20 38.40 43.97 2.33 
3.95 169.2 181.1 3.25 —3.4 175.4 3.07 
5.67 349.0 366.9 2.60 —3.3 355.6 2.56 
7.50 459.2 480.4 1.57 —3.4 464.8 1.59 
8.48 455.0 476.3 1.16 
9.47 418.6 439.0 1.14 
10.54 348.6 367.1 1.34 —3.6 354.8 
12.8 191.1 204.4 2.00 —3.2 198.1 1.96 
15.2 79.38 87.34 3.01 
17.8 24.74 28.71 3.26 
19.1 12.83 3.19 —2.2 15.18 3.18 


15 48 
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We believe therefore that for the Furry model 
with ionization the fluctuation of N is small and 
roughly equal to twice the Poisson value for all 
values of x. This result is quite different from the 
one reported in I (see table on p. 358), where for 
instance for z=4.75 and near the maximum of the 
multiplication curve ¢/N was found to be equal 
to 9.6. The origin of this discrepancy can be 
traced to the fact that in I the factor I'(s)I'(t)/ 
I'(s+t+2) was considered to be “slowly vary- 
ing.”’ This is not correct, because this factor has a 
pole for s=0 and t=0, and the saddle point 
So=to lies near zero when x is small. We have 
verified that as a consequence the correction 
terms [analogous to those in (24) ] become now 
quite appreciable and since the values of o are 
very sensitive with regard to errors in (N*),,, large 
errors in ¢ can be expected.* We have also verified 
the same fact by only taking the factor 1/st 
into the exponent, considering I'(s+1)I'(¢+1)/ 
I'(s+t+2) as “slowly varying.’’ The correction 
terms then again become quite small, and the 
results for o are in good agreement with the 
values of Table II. This is of importance for the 
actual cosmic-ray problem, since there it is 
quite complicated to take all variable factors into 
the exponent. 

The shape of the curve o/N as a function of x 
is quite curious and we have been unable to find 
a qualitative physical explanation for it. We 
believe that the shape is real and not due to 
possible errors in computation. It probably arises 
from the mathematical form of ¢. One can show 
that for very large x, ¢/N will approach one, but 
whether there will be further maxima in the 
curve is difficult to decide without numerical 
computations. 


§4. THE MODEL WITH q(E, u)=1 


One may look upon the effect of the ionization 
on the fluctuation as due to the fact that because 
of the ionization the particles lose their potency 
for producing pairs. Arley has considered a 
model in which the particles “‘degenerate’’ after 
a certain number of “generations” and he found 
that as a consequence the fluctuation soon reaches 


8 The pole at s=t=0 also explains the difference between 
Fig. 2 and the analogous Fig. 3 in I. By taking all variable 
factors into the exponent no special consideration of the 
region $x is necessary. 
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Fic. 3. The relative fluctuation ¢/N as a function of x for 
the Furry model with s=4.75. @=one term; +=two 
terms. 


the Poisson value (2a). It is therefore perhaps of 
interest to study a Furry-like model for which 
the splitting probability g(E£, ™) is a constant. 
This makes the total splitting probability pro- 
portional to the energy E of the particle, so 
that it will really diminish by the successive 
splitting processes. Furthermore, this model has 
the advantage that the influence of an ionization 
term B0F/dE can be studied exactly, although of 
course the model has very little similarity with 
the actual cosmic-ray problem. 

Without ionization one easily sees that P(N, x) 
will fulfill the equation : 


OP(N,x)/dx= — EoP(N,x)+EoP(N—1,x). (25) 


This is because the total energy is conserved and 
since the splitting probability of each particle is 
proportional to its energy, the sum of the 
splitting probabilities for all particles will be 
constant, and equal to the initial energy Eo, 
when the unit of x is so chosen that g(E£, u) =1. 
The solution of (25), with the initial condition 
P(N, 0) = 6. is given by: 


(Eox)*-! 
(N—1)!" 
so that NW =1-+ ox and (N®)y—(N)w= Ex = N —1. 


As soon as E,x>1 one gets therefore the Poisson 
distribution. 

With the ionization taken into account, we were 
unable to determine P(N,x), but from the 


P(N, x) =exp (— Eox) 
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general formula of §2 one can of course again 
compute the fluctuation. Equation (16a) be- 


comes for our case: 


OF ” OF 
—=—EF+2{ duF(u,x)+B—. (27) 
Ox E dE 


The solution, with the initial condition F(£, 0) 
= 6(E)—E) can be found exactly and is given by: 


F(Eo, E, x) =exp (— Eox+436x?)6( Ey — E— Bx) 
+x exp (— Ex—}Bx’) 


X[2+x(Eo—E—Bx)] (28) 
=0 (E+6x> Ep), 
from which follows: 
7 Eo—Bz 
N(Eo, x) -{ dEF (Eo, E, x) 
=[1+x(E»— Bx) ] exp (—36x*) (29) 


=0 (x>Ep/B). 


One then finds from (20), after a lengthy com- 
putation : 


(N?(Eo, x))w 
= 28x(2+ Eox — 28x") exp (ax) f exp (36t*)dt 
0 


+ (36x? — Eox —1) exp (— 38x”) 
+[462x4—4BEox*+ (Eo? —108)x?+4Eox+2 | 
Xexp (—Bx?); (x<Eo/28) 


= 28x(2+ Eox — 28x") exp (—Ax*) f exp (3(t?)dt 
+ (38x? — Eox —1) exp (— 38x") 


+2(1—x*) exp (8x? —2Eox+Eo?/28) ; 


Eo Eo 
(<2<x<—) 
26 B 
Eo 
=(; («>—). 
B 


(30) 
In Fig. 4 the results are plotted for N and for 
o/N; here t=Eox and 8/E¢ is taken equal to 
0.01; the dotted lines are the results when the 
ionization is neglected. One sees the same peculiar 
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Fic. 4. The average number of particles N and the rela- 
tive fluctuation o/N as a function of t= Eox for the model 
with g(E, u)=1; B/E¢? is taken equal to 0.01. The dotted 
lines show the corresponding quantities when the ionization 
is neglected. 


behavior as in the Furry model case, but the 
deviation from the Poisson value unity is con- 
siderably less. ¢/N presumably begins to ap- 
proach one for ¢ somewhat greater than 30. 


§5. THE COSMIC-RAY PROBLEM 


Before considering the fluctuations we shall 
collect here briefly the results for the average 
energy distributions of the electrons F(E, x) and 
of the photons (E, x). We shall mainly follow 
the notations and method of Landau and Rumer.® 
F and @ fulfill the equations: 


OF E 
—=-—F(E, x) f nr(E, u)du 
Ox 0 
+f duF(u, x)r(u, u—E) 
. 
+2) du®(u, x)y(u, E) 
E 
=Le F(E, x) +Le@(E, x), (31a) 
O® B 
—=-9#(E, » f y(E, u)du 
Ox 0 
+f duF(u, x)r(u, E) 
E 
=Lpo(k, x)+LeF(E, %). (31b) 


9L. Landau and G. Rumer, Proc. Roy. Soc. 166, 213 
(1938). 
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These are the analogues of Eq. (16) for the Furry 
problem. They must be solved with the initial 
condition F(E, 0) =6(Eo—E); &(E, 0) =0, corre- 
sponding to the fact that one electron is falling in. 
The term 80F/0E is omitted since we will take 
the ionization into account by the cut-off method. 
The thickness x is measured in radiation units; 
m(E, u)du is the probability per unit thickness 
that an electron or positron of energy E produces 
a photon of energy between u and u+du; 
y(E, u)du is the probability per unit thickness 
that a photon of energy E produces a positron- 
electron pair of energies u to u+du and E—u to 
E-—u—du; y(E, u)=y(E, E—1u) since no distinc- 
tion is made between electrons and positrons. 
The functions 7(E, u) and y(E, u) have again the 
homogeneity property : 


1 u 1 u 
m(E, u) -—x(-). y(E, )=—1(-) (32) 
E E E E 


and from the calculations of Bethe-Heitler follows 
that: 
4 4 4 4 
m(é)=——-+£&, 7(&)=-&—-E+1. (32a) 
3& 3 3 3 


Because of (32) Eqs. (31) can be solved again by 
the momentum method (see I, §2), and one finds: 


1 
F(Eo, E, x) =— fasE-—e6s, x), | 
2ri 





(33) 
1 
(Ey, E, x) =— fase ~*-19(s, x) | 
2m 
where 
(D—dA)e—** — (D—p)e~** | 
g(s, x) =Ep' | 
u—dX 
(33a) 


@(s, x) = Eo'——(e—* —e~**), 
u—~A 
\=1(A4+D)¥4[(A—D)?+4BC]! (33b) 
Bh) 


and A, B, C, and D are functions of s, given by: 


1 4 
A(s)= f ae(1—#)4(1- 9 = 14) +61] 


5 1 1 
—-+ + ’ 
6 3(s+1) s+2 





wn 
© 
uw 


. 2 8 8 
B(s)=2 f RNG nee, 
0 s+1  3(s+2) 3(6+3) 


| 4 4 1 
C()= f dgea(=—-—__ +, 
0 3s 3(s+1) s+2 

t 7 
D= { ae) =~. 
4 9 
Here C; is Euler’s constant and ¥(s) =I’(s+1)/ 
I'(s+1). The functions A, B, C, \, and yp, together 
with some of their derivatives, have recently 
been given by Rossi and Greisen.” 
For the average number of particles N (Eo, x) 
one then finds: 


Eo 
N( Es, x)= f dEF(Ep, E, x) 


1 D—X(s) 
=— fe—— —e-=4, (35) 
2ni s[u(s) —A(s) ] 


where z=log (Eo/e) and « is the cut-off energy. 
In F we have neglected the term with exp (— yx) 
since for the values of x in which we are interested 
it contributes at the most 0.7 percent to N, and 
generally considerably less. The value of € can 
be related to 8, the ionization loss per unit thick- 
ness, by means of the relation: 


f N (Eo, x)dx = Ey/8, (36) 


which is a strict consequence of the Eqs. (31) 
when the ionization term BdF/0E is added to 
(31a). Introducing (35) (and using the complete 
expression (33) for F) one obtains: 


Ey D "ad 


ds ———, 
8 2m s(AD—BC) 


From (34) it follows that s(AD—BC) has only 
one zero point, namely for s=1. By computing 
the residue there, one finds: 


E./B=[63/(144?+5) Je? 


or «=0.448. 


10B. Rossi and K. Greisen, Rev. Mod. Phys. 13, 240 
(1941). These authors have used slightly more accurate 
forms of the probability functions # and ¥, (32a), in which 
the factor 4/3 is increased by 0.027. With the exception of 
the second derivative, the resulting values of the functions 
differ from ours by 2 percent at the most. 
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The integral over s in (35) must again be com- 
puted by the saddle point method. The factor 
(D— X)/(u—X) can be considered to be ‘“‘slowly 
varying,” but the 1/s must be taken into the 


exponent. One obtains: 


D — X(So) 


sol u(S0) — (50) ] 


N(E£o, x) = 


— 


1 -3 
| 2"| -1"(s0x+—| | e802—2X(s0) | (37) 


So 
where x and So are related by: 


z—xnd' (so) — (1/50) =0. (38) 


§6. CALCULATION OF THE FLUCTUATION 
FOR THE COSMIC-RAY PROBLEM 


To calculate (N?)y one must first generalize 
the results of §2. We have to introduce three new 
distribution functions. F2(F1, Ex, x)dE,dE,2 de- 
notes as before the average product of the num- 
bers of electrons in the energy ranges dE; and 
dE2; &2(E,, Es, x)dE,dE2 denotes the correspond- 
ing quantity for the photons, while H2(i, Eo, x) 
XdE,dE, denotes the average product of the 
number of electrons in dE, and the number of 
photons in dE». The functions F2 and #2 are 
symmetric in E; and £2, and are singular for 
E,=E2; the function H2 is not symmetric but is 
regular for all values of E; and E>. As in (12) one 
can separate off the singular parts of F. and ®, 
by introducing two new functions K2(Fi, Eo, x) 
and J2(Fi, Es, x), such that: 


F,(E,, Es, x) = 6(E1— E2) F(A, x) 
+K.(F), Ea, x) 
?2( Fj, Ea, x) = 6( Fy o E2)*(E,, x) 
+J2(Fi, Ea, x). 


(39) 


Using the operator notation of Eqs. (31) we can 
now write the equations for the regular functions 


TABLE III. (N?)q, for the cosmic-ray case; = 5.67 and x =z. 











o =(N®) py 
N (N2)ay (Nav a/N °Furry %Poisson 
Present 16.12 309.4 49+16 3.0+1 243.7 16.12 
calculation 
Previous 17.0 442 153 9 272 17 
calculation*® 








® W. T. Scott and G. E. Uhlenbeck, Phys. Rev. 57, 1061A (1940). 


G. E. UHLENBECK 


Ko, He, and J, in the form: 


OKs mM 
nen (Le,+ Lee) K2(Fi, Es, x) 
Ox 
(2) 2 
+Lé2H2(Es, Es, x) +L¥\Hs(Es, Es, x) 
+2y(Ei+ Fo, E:)6(Ei+ Ee, x), 
OH, (yy (3) - 
—— = (Le+Le2)H2(Ei, Es, x) 
Ox 
| . (2 —" 
+L2K2(Ei, Ex, x) +L¥;Jo(E1, Ex, x) 
+2(Eit+Ffo, FE) F(E,+ £2, x), (40) 
dS» (3), (3) 
— = (L2,+Lk2)J2(Fi, Es, x) 
Ox 
(4) (4) 
+Le,H2(Fi, Es, x) +Le2H2(Es, Fi, x). 
These equations may be derived from a ‘‘master 
equation,” but they are almost self-evident from 
the physical meaning of the operators L“™; the 
inhomogeneous terms again give the correlation 
between the different energy intervals, because 


of the splitting processes. 
The solution of (40) is straightforward. One 


introduces the momenta 
k(s,t,x)=f  dEdE,EyEs'K2(E;, Es, x) 
0 


and the corresponding functions /(s, t,x) and 
j(s, t,x) formed with H, and Je. These fulfill 
then the differential equations: 


ak/ax = —[A(s)+A (t) ]R(s, t, x) 
+B(t)h(s, t, x) 
+B(s)h(t, s, x) 
+ S(s, t)0(s+4#, x), 
dh(s, t, x)/dx= —[A(s)+D ]h(s, t, x) 
+C(t)k(s, t, x) 
+B(s)j(s, t, x) 
+T(s, t)g(s+t, x), 
dh(t, s, x)/dx= —[A(t)+D h(t, s, x) 
+C(s)k(s, t, x) 
+B(t)j(s, t, x) 
+T(t, s)g(s+4#, x), 
0j/dx= —2Dj(s, t,x) +C(s)h(s, t, x) 
+C(t)h(t, s, x), 


(41) 
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Fic. 5. The integrand for the cosmic-ray problem, 
analogous to the function J(z, x, ¢) for the Furry problem. 
The two terms refer to the two terms of k(s, t, x) in Eq. 
(42) which have been considered. 


where we have introduced the new abbreviations : 


1 
S(s, t) = S(t, s)= 2f dé(1 — &)*E*F(E) 


2 P(s+1)r(¢+1) 
3° «(s+t+1) 
4(s+1)(¢+1) 
fs 
(s+t+2)(s+t+3) 








4ia) 
T(s, t)= i) de (1 —£)*¢%(8) 





ithaca t(4s+t+5) 
3 V(s+t+1) 


These equations must now be solved with the 


initial condition that each of the functions &, h, 
and j vanish for x=0; g and @ are given by (33a). 


The solution is quite cumbersome and we shall 
only indicate the result for k(s, t, x): 


Ey! 
r(s)r(t)r(s+t) 

X [| S(s, HC(s+4[A(s)-—DI DW —D] 

+T7(s, )B(t)[A(s +4) —D][X(s) —D] 

+T(t, s)B(s)[A(s +t) —D]D\@® —D)} 





k(s, t,x) = 


xf deexp {(¢—x)[A(s) +A) ]—EA(s+8)} 


+seven similar terms. (42) 


The seven additional terms are formed by re- 
placing in the first term A(s+?#) by u(s+2), or 
A(s) by u(s), or A(t) by u(t), or by any combina- 
tion of such replacements; r(s) is an abbreviation 
for u(s)—XA(s). 

From k(s, t, x) one finally finds (N*),, according 
to the equation: 


1 esthz 
Jf fasa k(s,t,x), (43) 
(2x1)? stEy*t' 


where as before z=log (Eo/e). The integrals must 
again be performed by the saddle point method. 
A precise calculation becomes quite involved, and 
we have therefore tried to simplify-it by the 
following two considerations. 

1. Only the exponential parts of the integrand 
and the factor 1/st are taken as “rapidly vary- 
ing.’’ From the experience gained in §3, one may 
expect that this will give reliable results. 

2. Only two of the eight terms in R(s, t, x) are 
taken into account, namely the one with X(s), 
A(t), A(s+2) (the first term) and the one with 
A(s), A(z), u(s+2). The reason for this becomes 
clear if one examines the exponential parts of the 
different terms, remembering that r(s)=,(s) 
—(s) is of the order unity or greater. For small 
values of — only the two terms mentioned are of 
the same order of magnitude; all other terms are 
much smaller. For x, some of them will be- 
come comparable with the first term, but then 
this term itself is quite small, so that the error 
made in neglecting the six additional terms all 
together is probably insignificant. 

The calculation has been carried out for 
z=5.67 and only for x=z; with 8B=90 Mev for 
air, this corresponds to an initial energy of 








(N?) ay = N+ 
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11.5X10° ev and to a depth slightly beyond the 
maximum of the multiplication curve N(Eo, x) 
as given by (37)." Table III summarizes the 
result and Fig. 5 shows the integrands as func- 
tions of ~ for the two terms considered. The 
previous calculations referred to in Table III are 
results reported at the Washington meeting of 
1940. At that time we did not know that it was 
important to include the factors 1/s in (35) and 
1/st in (43) in the rapidly varying part of the 


11 The maximum occurs at x=4.8 and has the value 17.5. 


BROXON 


integrand. Just as in I we found as a result far 
too great a value for the fluctuation. The accu- 
racy of the present calculations is hard to esti- 
mate and the limits given are more or less a 
guess. However, it seems sure that also in the 
cosmic-ray case the fluctuations are much smaller 
than the Furry value (2) and of the order of a 
few times the Poisson value (2a). 

One of us wishes to acknowledge the assistance 
of the Horace H. Rackham fund, without which 
much of this work would not have been possible. 
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Relations between 28-day fluctuations of intensity of the 
cosmic radiation and both terrestrial magnetic activity and 
sunspot areas were investigated. Definite pulses, both in 
the magnetic character and in sunspot areas, were found 
to be associated with the primary pulses in the cosmic 
radiation at Boulder, obtained by Chree’s ‘‘superposed- 
epoch” method. They were in general phase opposition to 
the cosmic-ray pulses, but the tip of the magnetic-character 
pulse preceded the tip of the opposite cosmic-ray pulse by 
one day; the lead was three or four days in the case of the 
opposed sunspot pulses. Similar relations were not found 
among secondary pulses, although a 34-day periodicity in 
sunspot-area pulses referred to days selected on the basis 
of cosmic-ray intensity was displayed. Direct application 


INTRODUCTION 


N a recent paper? the author reported a 
statistical investigation of cosmic-ray inten- 
sity fluctuations at Boulder (lat. 40° N; long. 
105°16’ W; alt. 5440 ft.) by Chree’s ‘‘superposed- 
epoch” method of analysis. This provided evi- 
dence for the existence of secondary pulses at 
about 28-day intervals both preceding and sub- 
sequent to the primary cosmic-ray pulses. These 
! Preliminary reports on some portions of this investiga- 
tion were made in a Letter to the Editor, Phys. Rev. 59, 
678.(1941), and at: the Lubbock, Texas, meeting of the 
Southwestern Division of the A.A.A.S., April 29, 1941; the 
Golden, Colorado, meeting of the Colorado-Wyoming 
Acad. Sci., Nov. 8, 1941; and the Detroit, Michigan, 


meeting of the Am. Phys. Soc., Feb. 20, 1942. 
2 J. W. Broxon, Phys. Rev. 59, 773 (1941). 


of Chree’s method to the magnetic character and sunspot 
areas, individually, indicated a 27-day periodicity in the 
former and a 34-day periodicity in the latter. A second 
method of investigation, used by Graziadei, Kolhérster, 
and others, was also employed. This yielded results in some 
respects contradictory to the first. In particular, it indi- 
cated 27-day fluctuations in sunspot areas in phase with 
the cosmic-ray fluctuations and out of phase with changes 
in magnetic character. However, it also indicated the 34- 
day periodicity in sunspot areas for the period of the 
investigation was more pronounced than the 27-day 
periodicity. Among other possibilities, the possible effects 
of sunspots through the agency of their magnetic fields 
were considered. 


secondary pulses represented deviations from the 
mean amounting to about 0.2 percent of the 
general average cosmic-ray ionization rate (cor- 
rected’ for barometric variations) of 38.19 ions 
per cc per sec. in a heavily shielded, high.pressure 


chamber. 


3In the paper of reference 2, this general average was 
incorrectly given as 38.16. It is not supposed that the 
difference exceeds the error of measurement of the absolute 
value of the ionization. However, 38.19 is nearer the 
average of the values used in the statistical investigation. 
Consequently, the upper pairs of curves in Figs. 1 and 2 
of reference 2, and the cosmic-ray curves of Figs. 1, 2, 6, 
and 7 of this paper, are drawn 0.09 percent too high. Con- 
fidence that few errors in the statistical work have gone 
undiscovered is due to the fact that all tabulations have 
been checked, and each step in the computations has been 
performed at least twice except in the case of a few of the 
probable error computations. 
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The purpose of the present investigation was 
to seek an explanation of the recurrence pulses in 
the cosmic-ray intensity. Because of the mag- 
nitude of the period indicated by the pulses, and 
in view of the present state of knowledge of the 
cosmic radiation, one is particularly inclined to 
suspect an association between these fluctuations 
and variations in the earth’s magnetic field and 
perhaps in solar conditions. Evidence of such 
relations has been provided both by theoretical 
and by experimental researches of other inves- 
tigators, some of which will be discussed later. 
Chree* devised his method of analysis for the 
purpose of investigating 27-day recurrences in 
terrestrial magnetic activity. 

Suspecting a relation between the magnetic 
pulses and solar activity, Chree proceeded to 
establish a relation between sunspot areas and 
the primary pulse in terrestrial magnetic vari- 
ables. His procedure has been adopted (and 
extended somewhat) in the investigation of the 
suspected relations in the present instance. 

It is presumed that the data and the explana- 
tion of the procedure and the technical termi- 
nology employed in the work of reference 2 are 
available. However, a very brief review of some 
of these seems necessary. The cosmic-ray data 
consisted of the average (corrected) cosmic-ray 
ionization during the (Greenwich) day for the 
interval from May 25, 1938, to December 1, 
1939, inclusive. There were no cosmic-ray data 
for 14 of the 556 days of this period. For the 
positive-pulse curves, the selected zero days 
were the five days in each full calendar month of 
the interval, during which the average cosmic- 
ray intensity was greatest. For the negative- 
pulse curves, the zero days were the five with 
least average cosmic-ray intensity in each 
month. The primary pulses were the largest 
pulses, occurring of necessity at day number zero 
in the final curves. 


RESULTS 


The terrestrial magnetic data employed in the 
present investigation were the world-wide or 
international magnetic character numbers com- 
puted by van Dijk® from data supplied by 40 to 


4C. Chree, Phil. Trans. Roy. Soc. A212, 75 (1913); 
A213, 245 (1914). 
5G. van Dijk, Terr. Mag. 44, 391 (1939) ; 45, 351 (1940). 


56 magnetic observatories. Van Dijk also selected 
the five days in each month which were mag- 
netically most calm, and the five which were 
magnetically most disturbed. Some indication of 
a relation between the cosmic-ray variations and 
the character numbers (representing magnetic 
activity) was provided by these selected days. 
For the 18 months of the investigation, there 
were 90 (magnetically) calm days and 90 dis- 
turbed days, and also 90 selected zero days of 
high cosmic-ray intensity and 90 of low cosmic- 
ray intensity. Among these, 20 magnetically 
calm days were identical with selected days of 
high cosmic-ray intensity, and 31 magnetically 
disturbed days were identical with selected days 
of low cosmic-ray intensity, appreciably higher 
numbers of coincidences than one would expect 
on the basis of random and unrelated distri- 
butions. In fact, for two different months, 3 of 
the 5 calm days were identical with high C-R 
intensity days, and for five different months, 3 
of the 5 disturbed days were identical with low 
C-R intensity days. On the other hand, there 
were only 6 coincidences of (magnetically) calm 
days with selected days of low C-R intensity, 
and 5 coincidences of disturbed days with high 
C-R intensity days. One instance of two coin- 
cidences in a single month occurred for each of 
these two types. 

In order to find whether geomagnetic activity 
might be associated definitely with the cosmic- 
ray pulses, new tables were formed in the same 
manner as in the investigation of reference 2, the 
only difference being that the character numbers 
were now employed instead of cosmic-ray inten- 
sities. The zero days were the identical days 
selected in the earlier investigation on the basis 
of cosmic-ray intensity. Character numbers for 
all the 556 days were used. 

The primary, positive, cosmic-ray pulse ob- 
tained with selected days of greatest intensity 
for all 18 months is shown by the broken-line 
curve of Fig. 1. The continuous curve represents 
the corresponding magnetic character-number 
values. This shows that a negative character 
pulse, descending 27 percent below the general 
average value of 0.753, is associated with the 
positive cosmic-ray pulse. The magnitude of this 
average of the 556 character numbers shows that 
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Days from Selected Days of High Cosmic-Ray intensity 


Fic. 1. Primary, positive, cosmic-ray pulse, together 
with corresponding pulse in world-wide magnetic character. 
The selected zero-days are the five of greatest average 
cosmic-ray intensity in each of the eighteen months, June, 
1938, through November, 1939. The arrows represent the 
largest of the probable errors of the average values of mag- 
netic character designated in this diagram, for an arbi- 
trarily chosen group of six m days. 


there was considerable magnetic activity during 
the period. The average recorded by Chree and 
Stagg® for all days of the twenty years, 1906 to 
1925, was 0.619. The lowest tip of the magnetic- 
character curve occurs at day, —1, or one day 
preceding the representative day of greatest C-R 
intensity. The shape of the curve at the tip 
indicates the minimum character value might 
have preceded the C-R peak by less than one day 
if a unit of time less than one day had been 
employed. 

Figure 2 shows the primary, negative, cosmic- 
ray pulse obtained with selected days of least 
intensity for all 18 months, together with the 
corresponding magnetic-character-number curve. 
The latter displays a well-defined peak at day, 
—1, rising 44 percent above the general average. 

The curves of Figs. 1 and 2 indicate a rather 
intimate inverse relation between the primary, 
cosmic-ray intensity pulses and terrestrial mag- 
netic activity. To be sure, the character pulses 
do not represent as great deviations from the 
mean as do the primary, magnetic-character 
pulses obtained by Chree with zero days selected 
on the basis of magnetic calm or disturbance. 
However, they are about as large as his first 
secondary pulses observed at 27 days before or 


* C. Chree and J. M. Stagg, Phil. Trans. Roy. Soc. A227, 
21 (1928). 
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after the primary, and are larger than his secon- 
daries of higher order. 

In view of this relation and the nearly equal 
periodicities displayed by the cosmic-ray pulses 
and Chree’s magnetic pulses, it was thought that 
the extension of the procedure to larger day 
numbers might show variations of magnetic 
character corresponding to the secondary cosmic- 
ray pulses. Accordingly, magnetic character- 
number tables were composed and averages 
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Doys from Selected Days of Low Cosmic-Ray intensity 


Fic. 2. Primary, negative, cosmic-ray pulse, together 
with corresponding pulse in world-wide magnetic char- 
acter. The selected zero days are the five of least average 
cosmic-ray intensity in each of the eighteen months, June, 
1938 through November, 1939. The arrows represent 
probable errors of magnetic character values as in Fig. 1. 


computed, for day numbers from = —5 to 135 
for zero days selected from the first fifteen 
months of the period, and for day numbers from 
n=5 to —135 for zero days selected from the last 
fifteen months. These were the identical zero 
days used in obtaining the cosmic-ray, sub- 
sequent-pulse, and previous-pulse curves of 
reference 2. Consequently, any magnetic-char- 
acter pulses yielded should be directly com- 
parable with the cosmic-ray pulses of Figs. 1 
and 2 in that paper. 

The magnetic-character-number curves ob- 
tained in this manner are shown in Figs. 3 and 4. 
Those corresponding to the cosmic-ray, subse- 
quent-pulse curves are shown in Fig. 3, and those 
corresponding to the cosmic-ray, previous-pulse 
curves, in Fig. 4. In either case, curve A cor- 






















































































GEOMAGNETIC AND HELIOPHYSICAL ACTIVITIES 511 
n- responds to the cosmic-ray, positive-pulse curve; magnitude to those at »=—1. However, the 
that is, the zero days for the curves A are those pulsations in the difference curves are smaller 
ial of high cosmic-ray intensity. The curves marked and quite irregular, in general, although two 
3eS B correspond to the cosmic-ray, negative-pulse rather large pulses appear in the subsequent 
at curves. In either case, the ‘difference’ is difference curve at large day numbers. Certainly 
ay obtained by subtracting the ordinates of curve _ the pulses do not display the striking decrease of 
tic B from those of curve A for the same day magnitude with order of displacement from the 
ic- number. primary, characteristic of Chree’s magnetic 
r= The curves do not provide any clear indication _ pulses. 
es of magnetic pulses corresponding to the secon- It was thought that any possible association 
dary cosmic-ray pulses. There are quite large between the secondary cosmic-ray and magnetic- 
fluctuations, in some instances comparable in character pulses obtained in this manner might 
: be displayed by a combination difference curve 
similar to that in Fig. 3 of reference 2. The solid- 
line curve of Fig. 5 was thus obtained, its 
ordinate at any day number » being the average 
of that of the difference curve of Fig. 3 at day 
number » and that of the difference curve of Fig. 
« 4 at day number —2, but in the combination 
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difference curve this average is expressed in 
percent of the general average. The corresponding 
cosmic-ray, combination difference curve is 
represented by the broken line. Although there 
still appears a large magnetic pulse nearly 
opposite in phase to the primary cosmic-ray 
pulse, these combination curves do not show 
such a relation persisting among the secondaries. 
In fact, there is some hint of pulses in phase 
opposition at the lower day numbers, with a 
gradual shift resulting in pulses in phase at day 
number, 112, or the fourth cosmic-ray secondary 
pulse. 

The next step was to find whether a relation 
between the cosmic-ray intensity and solar vari- 
ables might be brought out by this method. Of 
the several solar variables regularly observed and 
recorded by astronomers, it is doubtful whether 
any of them affords a proper measure of the 
heliophysical activities one might consider capa- 
ble of affecting more or less directly the intensity 
of the cosmic radiation at the earth’s surface. If 
one knew the instantaneous rates of discharge of 
electricity from the sun or the variations in the 
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Days from Selected Days of High or Low Cosmic-Ray intensity 





Fic. 5. The solid-line curve represents a combination of 
the “difference” curves of Fig. 3 and Fig. 4. The ordinate 
for amy day number 2, is the average of the ‘difference’ 
value for day number 2, in Fig. 3, and that for day number 
—n, in Fig. 4, expressed in percent of the general average 
character number. The broken-line curve is the curve of 
Fig. 3 of reference 2, and represents the similar combina- 
tion of cosmic-ray difference pulses. 
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solar magnetic moment, he might more optimis- 
tically undertake an investigation of such corre- 
lations. In this case, sunspot areas were chosen. 
The values employed were those supplied by the 
U.S. Naval Observatory in cooperation with the 
Harvard and Mount Wilson observatories, in the 
Monthly Weather Review.’ The numbers listed 
there and used here represent sunspot areas, 
corrected for foreshortening, in terms of one- 
millionth of the sun’s visible hemisphere as unit. 
The areas listed for individual spots and groups 
for a particular day were added to give the total 
sunspot area, irrespective of the positions of the 
spots. These total sunspot areas have been em- 
ployed in this paper as indicators of the helio- 
physical activity. Solar data for 21 days of the 
period of the cosmic-ray measurements were 
missing. Dr. Seth B. Nicholson of Mount Wilson 
Observatory kindly supplied data for 8 of these 
21 days from visual observations. Consequently, 
there were only 13 of the 556 days which were left 
blank. Dr. Nicholson has written that the time 
of maximum sunspot activity was 1937.4. Hence 
the period of the present investigation was during 
a time of decreasing sunspot activity, beginning 
just about one year after the maximum. 

The statistical work was carried out for the 
sunspot areas just as for the cosmic-ray inten- 
sities and the magnetic character numbers, the 
zero days still being the identical days selected 
on the basis of cosmic-ray intensity. The results 
obtained for sunspot areas in this manner are 
shown in Figs. 6-10, inclusive, which correspond 
to the magnetic-character curves in Figs. 1-5, 
inclusive, in order, and also in Fig. 11. In all of 
these the sunspot areas are expressed in percent 
of the general average for all days, which was 
1889.32 units. 

As in the case of magnetic character, Figs. 6 
and 7 show quite large pulses in sunspot area 
in approximate phase opposition to the primary, 
cosmic-ray pulses obtained with zero days for all 
18 months. The tip of the negative sunspot-area 
pulse corresponding to the positive, primary, 
cosmic-ray pulse occurs at day number x= —4, 
and is 14 percent below the average. The peak of 
the positive sunspot-area pulse corresponding to 
the negative, primary, cosmic-ray pulse, occurs 


7 Monthly Weather Review 66 (1938); 67 (1939); see 
tables for individual months. 
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Days from Selected Days of High Cosmic-Ray intensity. 





Fic. 6. Primary, positive, cosmic-ray pulse, together with 
corresponding pulse in total sunspot area. The selected 
zero days are the five of greatest average cosmic-ray in- 
tensity in each of the eighteen months, June, 1938 through 
November, 1939. The arrows represent the largest of the 
probable errors of the average values of sunspot area 
designated in this diagram, for an arbitrarily chosen group 
of six m days. 


at day number x= —3, and is 21 percent above 
the average. These sunspot-area pulses are much 
broader than the corresponding magnetic pulses 
of Figs. 1 and 2, in general agreement with 
Chree.* (In comparing widths of the cosmic-ray 
pulses in Figs. 1, 2, 6, and 7, it should be kept in 
mind that the cosmic-ray curves in all of these 
are drawn 0.09 percent too high.®*) 

Rather surprisingly, the secondary sunspot- 
area pulses of Figs. 8 and 9 are in general larger 
and more clearly defined in relation to the 
primary than are the corresponding secondary 
magnetic pulses of Figs. 3 and 4. They do not, 
however, show the 28-day periodicity of the 
cosmic-ray secondaries, as may be seen by com- 
parison with the abscissae numbered in the 
diagrams. There does seem to be an indication of 
a longer period. This is brought out more clearly 
in Fig. 10, where the lower curve shows the com- 
bination of the sunspot-area difference curves 
obtained in the usual manner, but inverted in 
order to facilitate comparison with the cosmic- 
ray, combination, difference curve, shown above 
it. The secondary sunspot-area pulses in this 
appear to be about as well defined as do the 
cosmic-ray pulses, but as indicated in the 
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Fic. 7. Primary, negative, cosmic-ray pulse, together 
with corresponding pulse in total sunspot area. The selected 
zero days are the five of least average cosmic-ray intensity 
in each of the eighteen months, June, 1938 through No- 
vember, 1939. The arrows represent probable errors of 
sunspot-area values as in Fig. 6. 


diagram they display a periodicity of about 34 
days. 

Figure 11 shows a combination of the sunspot- 
area difference curves which does not correspond 
to combinations made heretofore. Because the 
peaks of the primary sunspot-area pulses oc- 
curred at (or about) three days preceding the 
cosmic-ray primaries, it was decided to reflect 
the difference curve of Fig. 9 horizontally at 
day number n= —3, before combining it with 
the difference curve of Fig. 8. The combination 
was then inverted to give Fig. 11. Oddly enough, 
not only the primary but the first couple of 
secondary peaks and troughs appear sharper in 
this than in the combination curve of Fig. 10. 
The curve of Fig. 11 does not provide as good an 
indication of a 34-day period, although the first 
couple of pairs of pulses fit well into a 35-day or 
36-day period. 

In view of the irregularities of the secondary 
magnetic pulses and the unexpectedly long 
period indicated by the sunspot pulses, it was 
suspected that there might have been extra- 
ordinary abnormalities in the behavior of the 
terrestrial magnetic field and the sunspots during 
the interval of this investigation. Consequently, 
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Fic. 8. Primary and subsequent sunspot-area pulses as 
related to primary cosmic-ray pulses. For curve A, the 
zero days are the five of greatest average cosmic-ray intensity 
in each of the first fifteen months of the period specified in 
Fig. 6. For curve B, the zero days are the five of least 
average cosmic-ray intensity in each of the same fifteen 
months. The arrows represent probable errors as in Figs. 6 
and 7, but in this case the largest for twelve arbitrarily 
chosen m days is represented. 


it was decided to subject the magnetic character 
for the period to precisely the same procedure 
employed by Chree, and then to proceed simi- 
larly with the sunspot areas. The results of this 
are shown in Figs. 12 and 13. 

The zero days used in obtaining Fig. 12 were 
those selected by van Dijk.5 For curve A, the 
zero days were the five magnetically disturbed 
days for each of the first fifteen months, June, 
1938 through August, 1939, as usual in obtaining 
subsequent-pulse curves. For curve B, the zero 
days were the five magnetically calm days for 
each of the same months. 
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Fic. 9. Primary and previous sunspot-area pulses as 
related to primary cosmic-ray pulses. As in Fig. 8, the zero 
days for curve A are the selected days of greatest average 
cosmic-ray intensity, while those for curve B are the selected 
days of least average cosmic-ray intensity, but in this in- 
stance zero days for the last fifteen months of the eighteen- 
months period are used. The arrows represent probable 
errors as in Fig. 8. 


The difference curve is yielded by subtracting 
ordinates of curve B from the corresponding ones 
of curve A. The secondary pulses give fair indica- 
tion of a period of 27 days, in agreement with 
Chree. Their relative irregularities might be 
attributable to the brevity of the time interval 
of this investigation, but Chree and Stagg® found 
a very high degree of regularity of the period 
indicated by the peaks of their pulses, even for 
individual years. They display a marked differ- 
ence from Chree’s magnetic secondaries in that 
all the secondaries remain about equal in mag- 
nitude instead of decreasing with displacement 
from the primary. The ratio of the magnitude of 
a secondary to the primary pulse is comparable 
to the corresponding ratio for Chree’s secondary 
pulses at 54 days from the primary. This con- 
stancy of amplitude of the magnetic secondaries 
of Fig. 12 reminds one of the cosmic-ray second- 
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aries. It will be noted that the secondaries are 
not larger than those obtained with zero days 
determined by cosmic-ray intensity. 

To investigate further the changes in sunspot 
area, new zero days were selected. These con- 
sisted of the five with largest total area, and the 
five with least total area, in each month. In one 
instance the same area was found for a fifth and 
a sixth day. In this case the one nearer the other 
four of the group was selected. Among the groups 
of 90 days selected on the basis of sunspot area 
and the groups selected on the basis of cosmic-ray 
intensity, there did not exist as striking a relation 
as was found among the groups of cosmic-ray 
zero days and the magnetically calm and dis- 
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Fic. 10. The upper, broken-line curve is the cosmic-ray, 
combination difference curve shown in Fig. 5, and in Fig. 3 
of reference 2. The lower, solid-line curve represents a 
combination of the “‘difference’”’ curves of Fig. 8 and Fig. 9. 
The ordinate for any day number 2, is the average of the 
“difference” value for day number 2, in Fig. 8, and that 
for day number —2, in Fig. 9. Positive values are plotted 
downward and negative values upward, thus inverting the 
sunspot-area, combination difference curve. 
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n for "Subsequent’ and (-n-6) for “Preceding” Difference Values 


Fic. 11. Special, sunspot-area, combination difference 
curve, inverted. The ordinate for any day number 2, is 
the average of the “difference” value for day number 2, 
in Fig. 8, and that for day number (—n—6), in Fig. 9. 
Negative differences are plotted upward. 


turbed days. However, 20 selected days of large 
sunspot area were found to be identical with 
selected days of low cosmic-ray intensity, and 19 
selected days of small sunspot area were identical 
with selected days of high cosmic-ray intensity. 
On the other hand, there were 9 coincidences of 
selected days of large sunspot area with selected 
days of high C-R intensity, and 11 coincidences 
of selected days of small sunspot area with 
selected days of low C-R intensity. 

In Fig. 13 are shown the sunspot-area curves 
obtained with zero days from the first fifteen 
months, selected on the basis of sunspot areas. 
As anticipated, the primary pulses of Fig. 13 are 
considerably larger than those of Figs. 8 and 9. 
However, the secondary pulses are rather smaller 
and less regular than those referred to the cosmic- 
ray zero days. In the difference curve there 
seems to be some indication of the 34-day period 
for two intervals, but not thereafter. There is 
some indication of it even after three 34-day 
intervals in the positive-pulse curve A, but the 
irregularity of the negative-pulse curve B 
destroys it in the difference curve. 

In view of the relations indicated by the fore- 
going analysis, it was decided to use in addition 
a very simple method of analysis which has been 
employed by T. and B. Diill,* Sanford,® Grazia- 


~ 8 T. and B. Diill, Virchow Arch. 293, 272 (1934). 
® F. Sanford, Pub. Astronom. Soc. Pac. 47, 180 (1935). 
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Days Subsequent to Magnetically Disturbed or Calm Days. 


Fic. 12. Primary and subsequent magnetic-character 
pulses. For curve A the zero days are the five magnetically 
disturbed days selected by van Dijk for each of the first 
fifteen months. For curve B the zero days are the five 
magnetically calm days selected by van Dijk for each 
month of the same period. 


dei,!° and Kolhérster," for instance. In _ this 
method, if a periodicity of a certain magnitude 
in a certain variable is suspected (and perhaps 
relations among several variables), then the 
whole time of investigation is divided into seg- 
ments of that length and a sort of average or 
typical curve for the variable (or variables) is 


10H. T. Graziadei, Sitz. Akad. Wiss. Wien [2a ] 145, 495 
(1936). 

1 W. Kolhérster, Physik. Zeits. 40, 107 (1939). Kol- 
hérster observed 27-day fluctuations of +0.5 percent, out 
of phase with the sunspot relative numbers, flocculi, etc. 
He- also indicated he was able, by observing changes in 
cosmic-ray intensity, to predict the beginning or at least 
the disposition toward great magnetic disturbances of 
several days’ duration which may seriously interfere with 
radio communication. The writer has not yet been able to 
confirm this. 
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Fic. 13. Primary and subsequent sunspot-area pulses. 
For curve A the selected zero days are the five of greatest 
total sunspot area in each of the first fifteen months. For 
curve B the zero days are the five of least total sunspot area 
in each of the same months. 


obtained for the period being considered. For 
instance, to investigate by this method whether 
there is a 27-day periodicity displayed by the 
cosmic-ray data herein used, the entire interval 
of 556 days is divided into consecutive groups of 
27 days each, beginning with May 25, 1938 as the 
first number-one day. In Bartels’ * harmonic-dial 
method of analysis, the cosmic-ray data for each 
of these periods would be analyzed to two or more 
terms by Fourier analysis, and the amplitudes 
and phase constants corresponding to the 
assumed fundamental (or any particular har- 
monic) designated for successive intervals by the 





By, Bartels, Terr. Mag. 40, 1 (1935). 
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harmonic dial. In the more simple analysis being 
considered here, however, the simple average 
cosmic-ray intensity for all days of the same 
number in the several intervals is determined. 
These averages are then plotted against the day 
numbers to yield a sort of average or represen- 
tative curve which, by its shape, is supposed to 
enable one to decide concerning the “reality’’ of 
the suspected 27-day periodicity of the cosmic- 
ray intensity. 

The results of this treatment are shown in Fig. 
14, which contains curves for the magnetic 
character-figure, the cosmic-ray intensity, and 
the sunspot area, obtained upon the assumption 
of a 27-day periodicity in each of these variables. 
Each curve is drawn twice, as if the representa- 
tive curve were repeated for two cycles, in order 
to facilitate comparison of phase relations. 

The curves of Fig. 14 appear to indicate a 
27-day periodicity in all three variables. The 
cosmic-ray curve indicates average fluctuations 
of the order of 0.2 percent from the average, 
comparable to or larger than the cosmic-ray 
secondaries brought out by Chree’s method. The 
deviations from the mean displayed by the 
sunspot-area curve are of about the same order 
as the secondary pulses yielded by the Chree 
method, while those displayed in the magnetic 
curve are considerably larger. A striking feature 
of the diagram is the similarity of the cosmic-ray 
and the sunspot-area curves. These indicate 
27-day variations of these two variables almost 
exactly in phase. Not only is this true of the large 
variations, but in several instances smaller irregu- 
larities of the same type appear in phase in the 
two curves. Not only was this entirely unexpected 
in view of the generally inverse relation and the 
phase difference displayed in Figs. 6 and 7, and 
the inverse relation between cosmic-ray intensity 
and solar activity found by Kolhérster" by the 
same method, but comparison of the top and 
bottom curves of the figure shows an inverse 
relation between sunspot areas and terrestrial 
magnetic activity, in direct opposition to long 
accepted relations between these variables. 
Although the magnetic-character curve is quite 
irregular (particularly the positive loop) the 
main variations are generally out of phase with 
those in the cosmic-ray curve, as was expected. 


In fact, if the sunspot-area curve could be 
inverted, all the relations among the three vari- 
ables would be less surprising. 

An unsatisfactory characteristic of this method 
of analysis is that it requires the postulation of a 
particular period a priori. If one is not satisfied 
with the results obtained with an assumed 
period, he may start all over with another. Upon 
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Fic. 14. For these curves, the entire period of observation 
was divided into successive intervals of twenty-seven con- 
secutive days, each. The average values of the international 
magnetic-character figure, the cosmic-ray intensity, and 
the total sunspot area, were determined for all days of the 
same number in the several twenty-seven-day intervals. 
These averages are represented by the three curves, 
abscissae representing the order of days in an interval. 
Two complete cycles are included to facilitate comparison. 
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Fic, 15. These curves correspond to those of Fig. 14, but 
in this instance the successive intervals consist of twenty- 
eight days. 


comparison of the shapes of the final curves and 
the amplitudes of the variations they indicate, 
he may decide one is preferable to another, but 
he can hardly be sure that an untried period 
would not be still more satisfactory. Because of 
the relation of the sunspot curve to the others 
when the 27-day period was assumed, it was 
decided to repeat the analysis with other periods. 
Since the Chree method had indicated a 28-day 
period for the cosmic-ray pulses and a 34-day 
period for the sunspot pulses, it was decided to 
repeat the procedure with each of these intervals. 

Figure 15 shows the average curves obtained 


BRIOXON 


for the three variables upon arranging the data 
in groups corresponding to successive intervals of 
28 days. In this the cosmic-ray curve is perhaps 
a little less regular, but regarding both regularity 
and amplitude, it seems about as satisfactory as 
the cosmic-ray curve of Fig. 14. Although the 
fluctuations are still quite large, the magnetic- 
character curve of Fig. 15 is so irregular that it 
can scarcely be considered to provide any evi- 
dence in favor of a 28-day period. The sunspot 
curve of Fig. 15, however, is far more regular 
than that of the preceding diagram, and the 
amplitude of the variations is greater. The phase 
relations between the cosmic-ray and sunspot- 
area variations in the two diagrams are quite 
different. In Fig. 15 the sunspot-area fluctuations 
appear to lead the cosmic-ray fluctuations by 
about 90° or some 6 or 7 days instead of being in 
phase as in Fig. 14. 

In Fig. 16 are shown the results obtained with 
a postulated period of 34 days. Although the 
fluctuations in the cosmic-ray curve are still quite 
large, both this and the magnetic-character 
curve are so irregular that it seems one must 
conclude there is no evidence of a 34-day peri- 
odicity in these variables. However, the sunspot- 
area curve is reasonably regular, and the extreme 
variations are somewhat greater than in either of 
the two preceding diagrams. If the magnitude of 
extreme variations from the mean were the only 
criterion, we should conclude that the method 
provides better evidence for the 34-day, than for 
either the 27-day, or the 28-day period in the 
sunspot-area variations. However, the 28-day 
curve is somewhat more regular than either of 
the others.* 

When the three diagrams are considered 
together, they indicate a 27-day periodicity in 
the magnetic character, a 27-day or a 28-day 
periodicity in the cosmic-ray intensity, and while 
there is some indication of all three periodicities 
in the sunspot area, the 28-day and 34-day 
periods are more strongly indicated in this case. 
Different phase relations between the cosmic-ray 
pe * Note added in proof: Sunspot-area curves for data 
arranged in groups corresponding to successive intervals 
of 29, 30, 31, 32, and 33 days, respectively, have been con- 
structed by an assistant, Mr. S. W. Rasmussen. Of these, 
the 30-day, and the 33-day curve actually display some- 
what greater extreme variations than does the 34-day 


curve of Fig. 16, but are quite as irregular. In fact, none 
is as nearly symmetrical as the 28-day curve of Fig. 15. 
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ata and the sunspot fluctuations are indicated in the recorded intervals varying from 26.63 to 28.5 
s of different diagrams. When we consider the two days. The 27-day (or 28-day) variations in the 
aps statistical treatments together, there appears to cosmic-ray intensity have been announced by 
‘ity be fair agreement in the indication of a 28-day Hess,’® Graziadei,!® Gill,"® Monk and Compton,"” 
as (or possibly a 27-day) period for the cosmic Kolhérster,"' and Forbush.'* A sunspot (sunspot 
the radiation, a 27-day period for the magnetic _ relative numbers rather than areas are commonly 
tic- activity, and a 34-day period for the sunspots. used) period in the neighborhood of 27 days 
t it appears to be generally accepted, although the 
‘Vi- DISCUSSION OF RESULTS 
,0t i . , [es eeuaef.-a 
; lo explain the results obtained by the two fs @ &© § £ & G4 018 
° ° ° | ' . 
ma methods of investigation employed here, and ; ¥ a 
° ° | ! r 
the further, to explain the differences between these , | ; 1 ; Ii 
ase and the results of other investigators, does not hor | ! side ie ! | 120 s 
ot appear to be a simple matter. This may be due i NOOR IN Gs LAT 5 
ute primarily to the fact that the interrelations among nals 
eM the three variables (or four, including the time) ; 
»y are not simple. Perhaps some of the difficulties ! ia ; | 2 
r iy . . . ~ | A 
_ are due to the brevity of the time interval. In $}i003 | 1 | 6015 
; ; s l 
tl this connection, it should be noted that both Fl ! foi | boy 4 i 
“a methods have been applied to data for even F002 1 f 4 7 , * 
. " shorter intervals by some of the investigators : — + hes 
. a . 2) ! 1 
ate mentioned. That there are relations among the — 1) i | 
ani larger fluctuations of the three variables is . ooo Hitt | Ve 
. . . on ° s| vo v7 T 
nat indicated by Figs. 1, 2, 6, and 7, obtained by the 3 W aie 
i. 7 : - 
son first method. If the three variables actually stoos i ¥ ! 
all fluctuate with different periods (particularly if < ‘= ! aw 
. | 
ad one or more has only a temporary or quasi §| ose , 3 Pfs 
° — ° . | 
we" periodicity) then the relations among the second- $ i | 2a 
° ° ° + | 
: ary pulses obtained by this method may be fairly - 3 Ay i 204 
' ! comprehensible. The different results obtained : 2 Tee 2 ! g 
* by the second method upon assumption of dif- ; : - 
“ ferent periodicities, illustrate the fact that con- ! ! nog 
- e . . ! ' -_ 
aa siderable caution should be used in employing 1 ut Saget Gasans ! . 
~ this method, particularly for rather short inter- a ae a. # 2. 7 5 
of : + . - t : t 1004 
vals of time. ; | ! 1a € 
Other investigators have announced peri- ! ! : 
ed °_¢,¢ ° ° . ° \ ' s 
; odicities in the several variables. Not only did ! ee i | gl 
” Chree*:* arrive at the 27-day period in magnetic 1 | | > 
“ activity by his method, but his findings have ! 5 
> ° ! ! ! 
- been (approximately) confirmed by Bartels," iy . 6 . | go” 
. os ' ! 
os T. and B. Diill,* and Benkova,“ among others. oe ee ee 1 
id Benkova states that the 27-day recurrence ,o 2] 2 * S > + i 
se, “ . Days of the Interval. Hs 
tendency of magnetic storms was detected as ' 
a) early as 1858. He found the probability of such _ Fic. 16. These curves correspond to those of Fig. 14, but We 
. ‘ in this instance the successive intervals consist of thirty- iit 
- recurrence to be greatest during years of decreas- four days. A 
als ° ee . we 
. ing activity, which presumably W ould apply to 8 V. F. Hess, Terr. Mag. 41, 345 (1936). 4 
se, the interval of the present investigation. He 16 P. S. Gill, Phys. Rev. 55, 429 (1939). Ai} 
ne- 17 A. T. Monk and A. H. Compton, Rev. Mod. Phys. 11, 4 
ay peti: 173 (1939). nit! 
ae 8 J. Bartels, Terr. Mag. 37, 1 (1932). 18S, E. Forbush, Bull. Int. Union Geod. and Geophys. + 
144.N. P. Benkova, Terr. Mag. 47, 147 (1942). No. 11, 438 (1940). Hh, 
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period varies with time and for individual spots. 
Among those on the subject are papers by 
Sanford,® Alter,'® and Clayton.?® Alter observed 
a sunspot periodicity of 28 days in late 1936 and 
early 1937. While Clayton found an average 
period of 27.26 days, she observed a period of 29 
days in late 1935 and early 1936. Sanford con- 
cluded the sunspot period was 27.25 days. So 
far as the writer is aware, there is no confirming 
evidence for a 34-day sunspot period. Apparently, 
this could correspond to solar rotation only in 
the neighborhood of the poles, well out of the 
sunspot zone. 

The lag of 3 (or 4) days of the cosmic-ray, 
primary pulses relative to the inverse, sunspot 
pulses of Figs. 6 and 7 is not without its counter- 
parts. Chree* found a 4-day lag of the magnetic 
(H-range) pulse corresponding to his primary, 
sunspot-area pulse for the years, 1890-1900, 
inclusive. T. and B. Diill® observed a 5-day lag 
of the magnetic character relative to the sunspot 
relative number.” According to Chapman and 
Bartels,” Maurain found for the forty-one-year 
period, 1883-1923, that days of maximum sun- 
spottedness preceded days of great magnetic 
disturbance (selected by him) by some 2 to 4 
days, the average being about 2.5 days for the 
41 years. Also, according to them, Stagg found 
the time of maximum (projected) sunspot area 
preceded selected days of high magnetic dis- 
turbance by 1.67 to 2.25 days during the years, 
1901-10. Graziadei,'!® employing a period of 27.2 
days and extrapolating from Sanford’s® work, 
concluded that his cosmic-ray maximum must 
have lagged the sunspot-relative-number maxi- 
mum by about 4.5 to 7.5 days. It has been sug- 
gested that electrical particles may be emitted 
from sunspots or flocculi or perhaps other active 
regions only loosely associated with sunspots, 
that these alter the ring current about the earth, 
that this alters the terrestrial magnetic field, and 
that this finally alters the cosmic-ray intensity 





19 —P). Alter, Pub. Astronom. Soc. Pac. 49, 242 (1937). 

20H. H. Clayton, Terr. Mag. 46, 71 (1941). 

21]f the magnetic-character and cosmic-ray curves of 
Fig. ‘14 are accepted, then it might be inferred that the 
times of maximum death rate designated by the curves of 
T. and B. Diill would be times of decreased cosmic-ray 
intensity. 

2S. Chapman and J. Bartels, Geomagnetism (Oxford 
University Press, 1940), Vol. 1, pp. 382-3. 
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at the earth’s surface. If the train of events occurs 
in this order, then some of the phase lags might 
be explained in terms of the time required for the 
transit of the electrical particles, or perhaps by 
the supposition that the particles may be able to 
reach the earth or its ring current only in case 
they start from solar regions with certain 
locations relative to the earth which might cause 
their chief effect to be established some time after 
their culmination, as suggested by Graziadei.'° 
The sun may be able to affect the cosmic radi- 
ation through agencies other than the ring 
current. Vallarta,” Janossy,** and Epstein** have 
shown that the terrestrial intensity of the cosmic 
radiation may be influenced directly by the solar 
magnetic field. Vallarta and Godart,?* assuming 
a solar magnetic dipole of 10* gauss cm’ with its 
axis inclined to the solar axis of rotation, and 
assuming a 27-day period of rotation, arrived 
theoretically at a 27-day variation of appropriate 
magnitude in the cosmic-ray intensity at the 
earth. They do not appear to have stated what 
they assumed to be the angle between the sup- 
posed solar dipole and the sun’s axis of rotation. 
They did write: ““Taking the angle between the 
earth’s dipole and the solar dipole as 6°, we have 
for the amplitude of the 27-day period the values 
given in Table I, .. . It should be emphasized 
that the calculated amplitude is quite sensitive 
to the angle between the two dipoles.’’ Although 
these statements are very definite, it does not 
seem they could have assigned this constant 
value to the angle between these two dipoles, 
since it would actually vary with three periods, 
and the magnitude does not seem appropriate. 
Therefore, the writer is inclined to suspect that 
they actually assumed the angle between the 
solar dipole and its axis of rotation to be 6°. This 
suspicion is strengthened by the fact that Abetti*’ 
has written: ‘‘From observations carried out at 
Mount Wilson from 1912 to 1914 at the epoch 
of minimum solar activity in order to profit by 
relatively tranquil conditions of the photosphere, 


23M. S. Vallarta, Nature 139, 839 (1937). 

*4L. Janossy, Zeits. f. Physik 104, 430 (1937). 

% P. S. Epstein, Phys. Rev. 53, 862 (1938). 

26M. S. Vallarta and O. Godart, Rev. Mod. Phys. 11, 
180 (1939). 

27G. Abetti, The Sun (translated by Zimmerman and 
Borghouts), (D. Van Nostrand Company, 1938), p. 222. 
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Seares deduced the following results: 


i=6.0°+0.4°, 
P=31.52+0.28 days.” 


4 was defined to be the inclination of the solar 
magnetic axis with respect to the axis of the 
rotation, and P, the period of rotation of the 
magnetic axis around the axis of rotation. 
Forbush'® has offered two objections to the 
conclusions of Vallarta and Godart: first, that 
the effect does not vary with latitude as their 
theory demands; and second, that the observed 
variations are not truly or permanently periodic. 
The second of these objections might be met if, 
instead of supposing it to be due to an inclination 
of the sun’s dipole axis to its axis of rotation, we 
attribute the effect to the fields of the rather 
evanescent sunspots, superimposed upon the 
field of the permanent solar dipole with its axis 
supposed coincident with the solar axis of rota- 
tion. According to our supposition regarding the 
assumptions of Vallarta and Godart, the 27-day 
periodicity in the cosmic radiation could be 
regarded as due to a solar dipole-moment com- 
ponent of 10* sin 6°=10* e.m.u. rotating with a 
period of 27 days in the sun’s equatorial plane. 
If, then, a sunspot or a combination of them 
could produce the equivalent of a dipole moment 
of 10** e.m.u. parallel to the equatorial plane, we 
could have the 27-day periodicity derived by 
Vallarta and Godart, but vanishing and reap- 
pearing with the sunspots. Or we might suppose 
the permanent solar-dipole axis to be inclined at 
3° to the axis of rotation so that the component 
in the equatorial plane would be only half as 
great as with the 6° angle. Then the correspond- 
ing sunspot moment would need only to be 
5X10” e.m.u. in order that the 27-day effect 
might be just as computed, at certain times, and 
vanish at others, depending upon the phase 
relations between the sunspot moment and the 
equatorial component of the permanent moment. 
Whether such a notion is worth considering 
would appear to depend upon whether the sun- 
spots could ever be regarded as magnetically 
equivalent to a dipole of such a magnitude. To 
investigate this possibility we must first decide 
upon a method of estimating the equivalent 
dipole moment of a sunspot. In doing this we 
must keep in mind the known characteristics of 


the magnetic field of a sunspot. Spectroscopic 
observations show that this field is most intense 
at the center of the spot and is there directed 
along the radius of the sun or the axis of the spot. 
At the edge of the spot the field becomes vanish- 
ingly small and is directed along the radius of 
the spot, tangent to the sun’s surface. In the 
region between the center and the edge of the 
spot, the total magnetic field intensity is propor- 
tional to (1—r*/a*), where r is the distance from 
the center of the spot to the point and, a is the 
radius of the spot. If we assume these observed 
values of the spot’s magnetic field actually 
represent conditions in a plane (the plane of the 
spot) normal to the solar radius through the 
center of the spot, then we may imagine the 
observed magnetic field of a unipolar spot to be 
represented fairly approximately by the field of 
a dipole with its axis coincident with the solar 
radius through the center and located below 
(nearer the sun’s center) the plane of the spot at 
a distance a/v2 from its center. A dipole so 
located would yield fields in the correct directions 
at the center and the edge of the spot. The total 
intensity in the plane of the spot would not vary 
with the distance from its center as specified 
above, but would be more nearly proportional 
to (1—r/a), and at the edge would still have a 
magnitude nearly 14 percent as great as at the 
center. However, the dipole moment computed 
in this manner might be regarded as providing 
an adequate approximation for computation of 
the order of magnitude of a sunspot dipole 
moment for present purposes. To provide a field 
intensity H at the center of a spot of radius a, 
the dipole moment would need to be v2a*H/8. 

As an example, let us consider one of the 
largest unipolar sunspot groups during the time 
of this investigation, Mount Wilson Number 
6618, with a maximum area of 1939 (millionths 
of the visible hemisphere) on Sept. 23, 1939 and 
a maximum field intensity*® of 7=3600 gauss. 
The spot radius corresponding to its area is 
10°(19)* cm if the sun’s radius is taken to be 
7X10'* cm. Using these values, we find the cor- 
responding dipole moment to be 5.2 X 10* e.m.u. 
This is smaller than the required magnitude, by 
a factor of 10 or 20. If, then, the approximations 





28 Pub. Astronom. Soc. Pac. 51, No. 304, 14 (1939). 























522 .. 


were adequate to yield the order of magnitude, 
and if the conclusions of Vallarta and Godart 
result from the assumption of an‘angle of 6° 
between the solar dipole axis and the axis of 
rotation, it would appear that the sunspots 
could not produce the observed 27-day peri- 
odicity in the cosmic-ray intensity directly by 
their own magnetic fields, for the area of this 
spot group was greater than the average daily 
total-sunspot area during the period of observa- 
tion. Nevertheless, it would be interesting to 
know whether cosmic-ray intensities would 
display a better or worse correlation with sunspot 
magnetic moments*® than with sunspot areas or 
relative numbers. 

There are perhaps other circumstances which 
might give rise to periodic and perhaps related 
fluctuations of the variables under investigation. 
Chapman and Ferraro*® have attributed mag- 
netic storms to a geomagnetic ring current (or 
cylindrical current) which they show can be 
stable at a distance of a few earth radii. This they 
show to be naturally oscillatory in the radial 
direction, and upon assuming reasonable values 
for the variables, they compute natural periods 
varying from 2 sec. to 10 min. Thus any kind of 
impulse which would cause an alteration of the 
radius of the ring current might set up fluctu- 
ations of the terrestrial field with a period quite 
independent of any periodicity of the original 
agent. Evans*' has recently provided an ex- 


29It may be worth noting that the distance from the 
earth to the dipole representing such a sunspot at low 
latitudes, would vary about one-fourth as much on account 
of the sun’s rotation as it would on account of the eccen- 
tricity of the earth’s orbit. 

30S. Chapman and V. C. A. Ferraro, Terr. Mag. 46, 1 
(1941). 

31 F. Evans, Phys. Rev. 61, 680 (1942). 
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planation for a relation between fluctuations in 
cosmic-ray intensity and in the earth’s magnetic 
field, through the agency of currents of slow 
interstellar electrons flowing around the earth 
in spherical shells under the influence of the 
earth’s electric and magnetic fields. The region 
of these shells is estimated to extend about as far 
as the distance to the sun. He suggests the pos- 
sibility of an oscillatory fluctuation in terrestrial 
magnetic field and cosmic-ray intensity. Though 
he considers the effect of the emission of particles 
from the sun, he does not consider the effect of 
a solar magnetic field. It would seem that sun- 
spots might, by means of their magnetic fields, 
produce pronounced effects upon such a system 
and, if it should prove to be naturally oscillatory, 
cause temporary fluctuations in the terrestrial 
field and the cosmic-ray intensity, with periods 
both equal to and differing from the periods of 
the spots. 

Janossy and Lockett® have mentioned objec- 
tions to the theories based on the magnetic field 
of the sun. They write that recent spectroscopic 
measurements of Evershed failed to indicate the 
presence of a magnetic field at the surface of the 
sun, and that Chapman has shown that a mag- 
netic dipole field proceeding from within the sun 
will generate electric currents in the sun’s 
ionized atmosphere, which may limit the out- 
ward extension of the field. 

The writer is indebted to Dr. Seth B. Nicholson 
of Mount Wilson Observatory for supplying 
information regarding sunspots, and acknowl- 
edges NYA assistance, chiefly in checking com- 


putations. 
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The spectra of the configurations JJ's are calculated by using Dirac’s vector model. A com- 
parison is made with the experimental data of Sc I, Y I, and Cu I. 


THEORY 


HE configurations /l’s are of the simplest 
instances of configurations in which a term 
of a particular kind occurs more than once, so 
that the diagonal-sum method does not suffice to 
give complete formulas for the electrostatic 
energies of the different terms. We shall show 
that with Dirac’s vector model we can reduce 
this question to a simple three-vector problem 
which may be solved without writing energy 
matrices. 
The interaction between the electrons / and s 
is! 

Wi.= FotGi, (1) 
where the upper sign holds when the spins are 
antiparallel and the lower sign when the spins 
are parallel. Following Dirac’s vector model we 
may substitute the operator? 


Pi.= —43[1+(e1-e,) ] (2) 
for the double sign, and write 
Wis= Fot+ PGi. (3) 


In the same way the interaction between the 
electrons / and /’ may be written 


Ww=9,t+PuTl 1, (4) 


where ®, and I; are linear combinations of 
Slater’s F*(nl, n'l’) and G*(nl, n'l’); the coeffi- 
cients of these linear combinations are functions 
of 1, l’, and L, and their general formulas were 
given in a previous paper.* 

The electrostatic energies of the configuration 
ll's are the eigenvalues of the operator 


W= F)(nl, n'’s) + Fo(n'l’, n’’s) 
+O,+(PuGitPrsGvtPuT rs). (5) 


1E. U. Condon and G. H. Shortley, Theory of Atomic 
Spectra (Cambridge, 1935), (which we shall denote by 
TAS), p. 197; the index of G is, however, incorrect. 

2P. A. M. Dirac, Quantum Mechanics (Oxford, 1935), 
Section 61, Eq. (35). 

3G. Racah, Phys. Rev. 61, 186 (1942). 
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In the quartet terms every P.» equals —1, and, 
therefore, 


W(*L) = Fot+ Fo’ +%.—Gi-—Gv—-Tlr. (6) 
In order to obtain the eigenvalues of (5) for the 
doublet terms it suffices to calculate the square 
of the operator in brackets: we take into account 
the identities‘ 
(4a)? =3 —2(o,-a) (7) 
andé 
(@4°0) (Ga°O-) + (Ga--)(Ga*G») =2(Gy-e-); (8) 
we obtain 
(PuGitPrGutPuT 1)? 
=G6?2+G,?+T 12+3[1+(¢1-¢,) + (er -@,) 
+(e:-er) )(GiGu+GTr+GrIr); (9) 
owing to the relation® 
2(e,-0,) +2(e, -@,) +2(e,-e7) =4S(S+1)—9 (10) 
and to the fact that for doublets 
S(S+1) =3/4, 
our result reduces to 


(Pi.GitPr.Gut+PuTl 1)? 
=67°+G,?+l?-GiGr-Glt—-GrT, (11) 


TABLE I. The configuration 4s4p3d of Sc I. (F2=315, 
G,=4041, Ga= 1340, G:=375, G3= 26.6.) 











Term Obs. Calc. Diff. 
sp(3P)d *F° 15877 15718 +159 
2p° 16052 16223 —171 
4pe 16132 16099 +-33 
spe 18542 18636 —94 
2 Fo 21063 21001 +62 
sp('P)d 2D® 24955 24859 +96 
2 FO 25665 25751 — 86 
Z 2pe 18808 
y 2%po 24955 aaate 
x 2pe 30662 


4 See reference 2, Section 61, Eq. (33). 

5 This identity follows readily from reference 2, Section 
19, Eqs. (54) and (55). 

6 See reference 2, p. 227. 
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and then 
W(??L) = Fot+ Fo’ +®14+(G?+Gr? 
+1 ,2?-—G.Gy —GT,—GrT 1). (12) 
For the particular case of the configuration pds 
the formulas are 
4P= Fy +7 F.—G,—Ga—Gi—21G3, 
*P= Fo+7F2+[G,?+Ga?+(Git21G;3)’ 
BSG GGG}. 
4D = Fy—7F,.—G,—Git3Gi—7Gs, 
*D = Fy—7F2+[G,?+Gi?+ (3G:1—7G3)* (13) 
aC Sipt (GG NIG,.—IG)?, 
4 F= Fy)+2F,—G,—Ga—6G1—Gs, 
2 P= Fy+2F2+[Gy?+Ge?+(6G1+G;)? 
— GyGa—(Gp+Ga)(6Gi+Gs) }}, 


where’ 
Fy= F°(np, ns) + F°(n'd, ns) + F°(np, nd), 
F,= F?(np, n'd)/35, (14) 


G,=G'(np, n"’s)/3, Ga=G?(n'd, n’’s)/5, 
Gi=G"(np, n'd)/15, G3=3G3(np, n’d)/24S. 

We calculated also in the same way the term 
energies of the configuration d°sp and obtained the 
following results: 
4P=F,)—7F.2—G,, 

*P= Fy —7F.+Gat10G,+[G,?+G2" 
+100G,?—G,Ga—10(G,+G.)G1}}, 
‘D= F,\+7F.—Gy,, (15) 
2D = Fy t+ 7F2+Gat[Gp?+Ga? —GyGa]}, 
4F=F,)—2F,—G,, 
*F= Fy—2F2+Git15G3+[G,?+G.* 
+ 225G;3? —G,Ga—15(Gp,+Ga)Gs }}. 


COMPARISON WITH EXPERIMENTAL DATA? 
Sc I 


If we except the doublets ?P, the remaining 
terms fit very well to the theoretical formulas and 








sp(?P)d 
sp(@P)d | Fo+7F2—G,+4(Gat+Gi+21Gs) 
WeP%)= s%p 1R'v3 
sp('P)d }(Ga—G,—21G3) v3 


where X is the unperturbed value of 4s*4p *P® and 
R’ = R1(4s4p, 4p3d), R= R*(4s4p, 3d4p), (17) 
R* being defined by TAS 8°8. 

7 This definition of G; differs by a factor 3 from that of 
TAS, but agrees with reference 3, Eq. (50). 


’ All experimental data are taken from R. F. Bacher and 
S. Goudsmit, Atomic Energy States (McGraw-Hill, 1932). 
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TABLE II. Sharing of the eigenfunctions between the 
low ?P° terms of Sc I. 





Term sp(8P)d stp sp(P)d 
z 2po 53% 42% 5% 
y 2p? 39% 24% 37% 
x Pe 8% 34% 58% 


the values of the interaction parameters obtained 
by least squares from these seven terms are in 
good agreement with the corresponding values of 
Till.® The results are given in Table I; the 
doublets of the same kind are there distinguished 
by the resultant spin of the electrons 4s and 4p 
instead of that of 4s and 3d, since G, is much 
greater than Gag and I, and, therefore, the spin 
coupling of 4s and 4p outweighs that of 4s and 3d. 
The purity of our correlation is about 95 percent, 
whilst the purity of the correlation to the states 
3d4s of ScII is only 89 percent for ?F and 53 
percent for *D. 

The calculated values of the two ?P terms fall 
almost exactly in the centers of the intervals 
between the three experimental *P° terms which 
are in that region of the spectrum, so that it is 
even impossible to say which of these terms 
belong to 3d4s4p. The classification of the re- 
maining term as 3d?('S)4p ?P® cannot be ac- 
cepted, since this term must be the highest of the 
configuration 3d74p and not the lowest ;° owing to 
its position, it may be assigned without difficulty 
to 4s*4p. 

In order to see if the proximity of 4s*4p *P® may 
produce displacements of about 3000 cm™! in the 
*P° terms of 3d4s4p, we calculated the configu- 
ration interaction between these terms. The re- 
sulting matrix is 


sp sp('P)d 
IR’ v3 }(Ga—G,— 21G3) v3 | 
v4 1R! —2R’ | (16) 


+R’ 2 R’ Fot+7 F2+Gyp- 3 (GatGi+21G;3) 


Introducing the values of the known parame- 
ters we have 

22047. ER’3 203v3_ || 
gR’'N3 x LR” —2R’ , (16’) 
203v3 1R’—2R’ 27855 || 


*G. Racah, Phys. Rev. 62, 438 (1942). Section 9. 
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and the requirement that the three eigenvalues of 
this matrix should agree with the three experi- 
mental doublets *P° gives us three equations for 
X, R’, and R”; the resolution of this system 
gives X = 24223 and 


R’=—2015, R’’=10831, (18a) 
or 


R’ = 8764, R”’ =9535. (18b) 
The linear combination 
S=a°R'(sp, ps) +R'(dp, pd) +2aR'(sp, pd) (19) 


may be considered as the integral G' for the p 
electron and an electron which has the radial 
eigenfunction aR(s)+R(d) ; S must, therefore, be 
positive for every value of a,” and this gives us 
the theoretical limitation 


| R'| <[G*(4s, 4p) -G4(3d, 4p) }! 
=[45G,Gi]!=8528; (20) 


the value (18b) of R’ must therefore be rejected. 
The analogous limitation for R’’ is 


|R” | <(F*(4s, 4p)- F°(3d, 4p) ]!; (21) 


unfortunately F?(4s, 4p) does not appear in the 
formulas of the term energies, but it is reasonable 
to assume that it must be greater than F?(3d, 4p) ; 
since the value (18a) of R”’ does not reach 
F?(3d, 4p) =35F2, we have no reason for re- 
jecting it, and the deviations of the *P® terms 
from the theoretical formulas (13) may be con- 
sidered as explained. 

It is interesting to consider the sharing of the 
eigenfunctions between the three doublets *P°; 
we calculated it with the parameters (18a) and 
the results are given in Table IT. We see that the 
sharing is so strong that it has no physical 


TABLE III. The configuration 4d5s5p of Y I. (F2= 399, 
Gp=3287, Ga=1716, Gi=550, G3=45.5.) 








Term Obs. Cale. Diff. 
dsp ‘4F 15690 15498 +192 
dsp ‘*D® 16855 16583 +272 
dsp ‘*P° 19079 19333 — 254 

ds('D)p 2 F° 21750 22245 —495 
ds(8D) p ?F° 24736 25447 —711 
dp(!D)s 2D° 16098 16188 —90 
dp(*D)s 2D° 24500 24322 +178 
sp(§P)d 2P® 24553 24155 + 398 


sp('P)d 2P° 28035 27527 +508 


10 See reference 9, Section 10. 


TABLE IV. The configuration 3d%4s4p of Cu I. (F:=315.5, 
Gp= 6906, Ga=901, G: = 391.5, G3=95.) 








Term Obs. Calc. Diff. 





sp(@P)d* *P® 22603 22589 —14 
‘Fe 20990 21012 +22 
‘pe 18180 18172 —8 
27 17875 17541 — 334 
spe 16448 16067 — 381 
2p° 15880 16867 +987 
sp('P)d® 2F° 5383 6019 +636 
spe 5290 5667 +377 
2pe 4841 3863 —978 





meaning to assign the different doublets to a 
definite configuration ; it is also remarkable that 
the term y *P°, the value of which differs very 
little from the unperturbed value of 4s°4p ?P°, 
shares only 24 percent of its eigenfunction. 


Ti I 


Two terms of 3d4s4p are unknown; the re- 
maining do not fit to the formulas, since the 
order of ?P and ‘P is inverted. 


zs 


The parameters were calculated by least 
squares and the comparison with the experi- 
mental data is given in Table III. It is remark- 
able that in this case the parameters of the spin 
coupling are of the same order of magnitude, so 
that the coupling of the spins is different for 
every value of L; the purities of the given spin 
couplings are 100 percent for *F and ?P and 97 
percent for *D. 


Zr Il 


The terms of 4d5s5p are so mixed with those of 
4d°5p that it is impossible to fit them to the 
theoretical formulas without considering the two 
configurations together, and also their inter- 
action. 


Cu I 


The intervals between the quartets have 
almost exactly the theoretical ratio, but the 
doublets are perturbed by the terms of the 
normal series and fit only roughly to the formulas 
(15). Owing to this fact, we calculated at first 
F)—G, and F, from the quartets, and then the G* 
from the doublets. The results are shown in 


Table IV. 
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The arc spectrum of bismuth was investigated in three 
orders of the best 30-foot Chicago grating. The light source 
consisted of a heated quartz tube containing bismuth metal, 
filled with a few mm of helium and excited by a 4-meter 
oscillator. At a sufficiently high vapor pressure of bismuth 
only the arc spectrum is excited and is so intense, that not 
only almost all previously known lines but also many new 
lines were studied. Because of the very high quality of the 
grating and the sharpness of the lines, pictures were ob- 
tained which show in the third order an exceptionally high 
resolution of about 400,000, which constitutes the highest 
resolving power ever obtained with a diffraction grating. 
The results of measurements of the hyperfine structures 


for 55 lines are presented in Tables I and IV; the wave- 
lengths for the components corresponding to the transi- 
tions between hyperfine structure levels with highest F 
quantum numbers, and the exact term values for the 
atomic energy levels derived from them are collected in 
Tables II and IV. Many new levels were identified; the 
total number of known levels is increased from 27 to 45. 
Hyperfine structure separations of levels are given in 
Tables II, III, and IV; all values previously reported by 
Zeeman, Back, and Goudsmit for 13 atomic levels were 
corrected and separations for 28 more levels were deter- 
mined. 





INTRODUCTION 


HE arc spectrum of bismuth has been the 
subject of many investigations. The energy 
level scheme was definitely settled by the clas- 
sical work on Zeeman effect and hyperfine 
structure by Goudsmit and Back and by Zeeman, 
Back, and Goudsmit.! There was not much hope 
that a reinvestigation of this spectrum would 
yield enough new results to make the work worth 
while. The present investigation was started in 
collaboration with Dr. F. Bueso-Sanllehi_pri- 
marily as a study of the resolving power of the 
Chicago 30-foot gratings. It was believed that 
the bismuth hyperfine structure patterns would 
constitute a very good source for such a study. 
This expectation was partly fulfilled, since pic- 
tures showing an unusually high resolution were 
obtained. The author, however, found a little 
later that a mercury hydride spectrum (mercury 
deuteride for regular 21-foot gratings) excited 
with high intensity in a mercury arc? constitutes 
in every respect a much more convenient source 
for such a study. The latter source is operated 
very easily; it is ready for use at any moment 
and the spectrum, extending from A2800 to 
4300, consists of many band lines with a very 
narrow structure (isotope effect of mercury) of a 
variable distance. 
1p, Zeeman, E. Back, and S. Goudsmit, Zeits. f. Physik 
66, 1 (1930) and the references there to previous work. 


? A description of such a powerful source is contained in 
the article S. Mrozowski, Zeits. f. Physik 95, 524 (1935). 


The closer examination of the first spectro- 
grams obtained of the bismuth spectrum revealed 
several points of interest. It was evident (espe- 
cially in the case of the line 43024) that some of 
the hyperfine structure separations given by 
Zeeman, Back, and Goudsmit! are not quite 
correct. Therefore a more extensive study of the 
spectrum was made by using a source which 
turned out to be extremely good for this inves- 
tigation. A few discrepancies between early 
works of Nagaoka and Mishima* and Zeeman, 
Back, and Goudsmit! were found to have an 
explanation; the hyperfine structure separations 
reported by Zeeman, Back, and Goudsmit! were 
corrected (13 levels) and separations in 28 more 
levels were determined. The number of estab- 
lished atomic energy levels of bismuth has been 
raised from 27 to 45 and their positions deter- 
mined with high accuracy. A very short account 
of a part of the results obtained was presented 
a year ago at the Chicago Physical Society 
Meeting.‘ 


EXPERIMENTAL 


In the early part of the experiments a water- 
cooled hollow cathode discharge tube was used. 
This type of source causes the inconvenient ap- 
pearance of both the arc and spark spectra at the 


3H. Nagaoka and T. Mishima, Proc. Imp. Acad. 2, 249 
(1926). 

4S. Mrozowski and F. Bueso-Sanllehi, Phys. Rev. 61, 
108 (1942). 
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ARC SPECTRUM 


same time. Because of a low density of bismuth 
vapor, only the strongest lines can be investigated 
and even in those lines very weak components 
often can hardly be detected. A quartz tube con- 
taining a piece of metal and filled with helium 
at few mm pressure excited by a short wave 
oscillator constitutes a much more satisfactory 
source for an arc spectrum of a metal with 
relatively high vapor pressure. Such a source has 
been studied in detail by the author in the case 
of the spectrum of lead® and the reader is referred 
to this paper for the description of the experi- 
mental arrangement and precautions to be taken 
in order to avoid the destructive influence of the 
condensed metal on the walls of the tube. A 
similar arrangement was adopted in this inves- 
tigation; the length of the quartz tube was, 
however, increased to about 25 cm. 

An increase of the vapor pressure of bismuth, 
which is regulated by the temperature of the 
oven surrounding the tube, has a double effect 
on the spectrum. First, the average velocity of 
the electrons is diminished which favors the 
excitation of the lower energy levels; by a suf- 
ficient increase of pressure the first spark spec- 
trum can be completely eliminated. Secondly, 
the rapidly increasing density of the vapor 
produces a great increase in the intensity of the 
weakest lines (or components); the strongest 
lines (or components) being more and more 
weakened by reabsorption. This is the reason 
that, for instance, the weakest components of 
the line 42938 could be investigated and gives 
also an explanation of the presence in the 
spectrum of many weak lines which have not 
been recorded previously by other authors. 

The spectrum was photographed end-on with 
the 30-foot grating spectrograph: \2850-\6500 
in the first, \2228-\4722 in the second, and 
\2228-A3067 in the third order. The slit width 
varied from 12, (strongest lines) to 25u (weakest 
lines). The exposure time ranged from a few 
minutes to about 48 hours. In the case of long 
exposures there was no difficulty in picking out 
a period of time with a sufficiently constant 
pressure (Ap <0.1 inch) ; a greater difficulty arose 
from long period variations of the grating tem- 
perature. Although the functioning of the tem- 


5S. Mrozowski, Phys. Rev. 58, 1086 (1940). 
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Jt 
Nm 
=J 


perature control arrangement was_ recently 
improved, there are still slight variations of 
temperature present which decrease the resolving 
power in the case of long exposures. Moreover, 
recently by using a powerful mercury arc? as a 
source, the author has found by very short ex- 
posures (2-5 min.) that at the time of switching 
on of the automatic heating arrangement a heat 
wave is created which in most of the cases 
decreases a little the resolving power of the 
grating for a period of a few minutes. It seems 
therefore quite certain that the highest resolving 
power obtained in this work (400,000 in the 
third and 300,000 in the second order) does not 
represent the true maximum resolving power of 
the best 30-foot grating (so called No. 3) in use. 

Secondary iron standards were used‘ as a 
comparison spectrum. Because of the width of the 
iron lines and the small density of standards 
(especially in some regions) the wave numbers 
obtained are less exact than should be expected 
for such a high dispersion (+0.02-0.03 cm, in 
some regions even worse). The distances of 
hyperfine structure components are obtained 
naturally with a much higher accuracy amount- 
ing to around +0.003 cm~ in favorable cases 
(completely resolved structures). 


RESULTS 


In Table I are given the results of measure- 
ments of hyperfine structures for lines corre- 
sponding to transitions to the five lowest levels 
of Bil belonging to the configuration 6p*. Lines 
observed only in one order are marked by crosses 
+. The positions of components are given in 10-* 
cm relative to the component corresponding 
to the transition between h.f.s. levels with 
highest F quantum numbers (in most cases this 
is the strongest component). The numbers in 
parentheses indicate the components in order of 
decreasing intensity. At low pressures the 
intensity distribution seems to be in almost all 
cases in agreement with the theoretical formulas 
of Hill; one exception will be mentioned below. 

Enlarged photographs of several lines are given 
in Fig. 1. The first picture (a) shows the com- 
pletely resolved pattern of the line \2627. This 


®See the discussion of this point in the paper S. 
Mrozowski, Phys. Rev. 60, 730 (1941). 
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(a) 42627 


(b) A2864 


TA 


(c) 42938 
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(d) 42993 (e) A3024 


Fic. 1. Enlargements of photographs of hyperfine structures of lines of bismuth obtained in the third order of the 
30-foot grating spectrograph. The picture (e) shows the highest resolution ever obtained with a grating (400,000). 


line has been investigated by Zeeman, Back, and 
Goudsmit! in the fifth order of a 21-foot Rowland 
grating and the structure found by them seemed 
to agree with the predicted pattern satisfactorily. 
But the comparison of their predicted pattern 
with Fig. 1a shows the absence of their fourth 
component from the left (wide gap after com- 
ponent 3 marked by an arrow in Fig. 1); this is 
exactly what should be expected, since according 
to the formulas of Hill this component should 
appear only with extremely low intensity. 

Zeeman, Back, and Gousmit! mention a dis- 
crepancy between their results and those of 
Nagaoka and Mishima*® in the case of the line 
2938. The study of this line at different pressures 
shows that there is no discord at all. Zeeman, 
Back, and Goudsmit used a well-cooled source 
giving an intensity distribution close to the 
theoretical, and were able to observe only the 
strongest components, group (1) and group 
(2)-(3). Nagaoka and Mishima had a source 
with a higher vapor pressure of bismuth and the 
strongest components were weakened by reab- 
sorption. This is seen also from the deviations 
from the theoretical intensity distribution in the 
case of other lines observed by them. In A2938 
they were able to observe even the weakest com- 
ponents. A picture of this line at conditions cor- 
responding to theirs is given in Fig. 1 (c), where 
the main components (1) and (2)—(3) weakened 
by reabsorption are marked with arrows. 

The high resolving power of the grating is 
illustrated also by the pictures of the lines \2864, 
2993, and \3024. The sharpness of components 
in the line \2993 indicates the extreme smallness 
of the total hyperfine structure separation in the 


upper level. The picture of the line 43024 (1-hour 
exposure) is especially interesting, since it shows 
the highest resolution ever obtained with a 
grating. Until now many spectra have been 
investigated in higher orders of diffraction grat- 
ings, but so far as the author is aware, even with 
the best gratings available a resolving power of 
more than 300,000 never was obtained in prac- 
tice. The distance of the components in the close 
doublet in the middle of the pattern (fifth and 
sixth components from the right, marked with 
arrows) corresponds to a resolving power of 
385,000; the nice resolution of this doublet in- 
dicates that the resolving power obtained in this 
case was probably even a little higher, maybe 
around 400,000. The high resolving power of the 
grating can be well appreciated by a comparison 
of the microphotometer tracings obtained for this 
line from pictures taken in the third order of the 
Chicago 30-foot grating and in the fifth order of 
a 21-foot Rowland grating (the last curve was 
taken from the work of Zeeman, Back, and 
Gousmit).! Both curves are drawn with the same 
scale in Fig. 2 (a); the contour observed by 
Zeeman, Back, and Goudsmit is resolved here 
into 11 components. The structure observed is in 
a very good agreement with the earlier results 
of Nagaoka and Mishima,’ but differs consider- 
ably from the pattern predicted by Zeeman, 
Back, and Goudsmit! [their predicted pattern is 
included in Fig. 2 (a) ]. In spite of such an ap- 
parently great discrepancy the predicted pattern 
was brought into agreement with the observed 
by only a slight adjustment of the separations 
reported by Zeeman, Back, and Goudsmit for the 
upper 75/2 and the lower 2D°5,2 levels (see Table 
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II), since this type of a pattern is very sensitive 
to a small change in the ratio of separations. In 
Fig. 2 (b) a comparison of microphotometer 
traces obtained for the line. \2697 in this work 
and by Zeeman, Back, and Goudsmit is pre- 
sented. A very interesting unexplained feature 
of this line is a too high intensity of the fifth 
component from the left, which should be weaker 
than the sixth according to the formulas of Hill. 

In general the data agree with those of former 
investigators very well. The wave-lengths of the 
components with highest F quantum numbers 
were determined with the highest accuracy and 
the results are given in the left part of the 
Table I. From the corresponding wave numbers 
obtained from Kayser’s Tabelle der Schwingung- 
szahlen the exact term values of the energy levels 
with highest F quantum numbers were found and 
are collected in Table II. The terms marked by 
a star are new terms established in this work; the 
levels denoted by (Th) were reported only by 
Thorsen,’ those by (T) were given only by 
Toshniwal.* All the remaining levels were re- 
ported by both of these authors. The numbering 
of the even terms 1 to 19 adopted by Bacher and 
Goudsmit® has been retained in Table IT; a level 
reported by Thorsen (20) and not included in the 
tables of Bacher and Goudsmit (they probably 
regarded it as insufficiently well established) was 
confirmed by the newly found line \3038. A new 
level 23 has been found, which almost coincides 
with a level reported by Thorsen on the basis of 
only one line observed (A2251, no combination 
principle involved). This line does not fit the 
level scheme exactly; on the other hand Tosh- 
niwal’s tables do not include such a wave-length, 
only a slightly different one \2249, which was 
used by him in identifying the level 18. Which of 
these classifications is the right one cannot be 
decided without remeasuring the position of the 
line in question. 

In order to make the data collected in Table II 
as complete as possible the wave-length of all 
other classified lines which could not be deter- 
mined in this work but which were reported by 
Thorsen’ or Toshniwal® are included in Table II 


7V. Thorsen, Zeits. f. Physik 40, 642 (1927). 

8G. R. Toshniwal, Phil. Mag. 4, 774 (1927). 

*R. F. Bacher and S. Goudsmit, Atomic Energy States 
(McGraw-Hill, New York, 1932), pp. 85-86. 
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and put in parentheses to indicate their lower 
accuracy. All newly classified lines are marked by 
stars. Three terms 16, 17, and 18 identified by 
Toshniwal® were added, although several doubts 
can be raised. The existence of the level 16 is 
extremely doubtful, since the lines in question 
for such a low level should show a high intensity ; 
in spite of this none of them was observed in this 
work. Further, if Toshniwal’s line 3239.7 is 
identical with a weak line \3241 observed in this 
work the classification is incorrect, since the 
h.f.s. of the latter line (a triplet) is incompatible 
with the structure expected. For the last level 
18 the classification of the line \2249 is uncertain 
(see above). 

The J quantum number of the level 12 was 
changed from 3 to }, since the line 43093 reported 
by Toshniwal* could not be found on my spec- 
trograms. An uncertainty exists about the J 
value of the level 19, since the structure found 
for the line \2448 seems to fit better with J=3, 
but on the other hand the absence of the tran- 
sition 19-+?D°;,2 seems to require J =}. 

The exact wave-lengths for the three lines 
\2276, \2230, and 2228 are also given in Table 
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Fic. 2. Comparison of microphotometer tracings ob- 
tained from plates taken in the third order of the 30-foot 
Chicago grating (continuous curve) and in the fifth order 
of a Rowland 21-foot grating (broken curve; redrawn from 
the paper by Zeeman, Back, and Goudsmit, reference 1). 
For the line 43024 the pattern predicted by Zeeman, Back, 
and Goudsmit (reference 1) on basis of their measurements 
of other lines is added. 
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TABLE I. Structures of Bi I lines. Transitions X—>6p%. 





Ain A Components in 10-3 cm~'} 

2400.965* —580(10)?, O(1), 242(9), 892(2), 1049(8), 1569(3), 1704(7), 2143(4), 2273(4), 2659(5), 3013(6) 
2448.683* —670(4), —525(3), weak diffuse (5), —270(2), 0(1) 
2496.678* 0(1), 605(2), 1113(3), 1524(4), 1847(5), 2097(6) 
2515.619* —616(4), —466(3), —253(2), O(1) 

2524.470* —762(4), —575(3), —330(2), 0(1) 

2594.771* 0(1), 588(2), 1089(3), 1502(4), 1837(5), 2084(6) 
2600.654* 0(1), 585(2), 1084(3), 1489(4), 1814(5), 2067(6) 
2627.834 —2020(4), —1880(4), —1640(3), —1292(3), —1027(2), —810(4), —560(3), —262(3), 0(1) 
2696.614 —3744(7), —3367(8), —3223(9), —2864(11), —2705(4), —2503(10), —2092(6), —1881(3), —1637(12), —1131(2), —874(S), O(1) 
2730.583* O(1), 125(2), 2007(2), 2101(1) 

2780.450 —779(2), —626(2), —449(2), —250(2), 0(1) 
2781.260* —242(2), 0(3), 1648(4), 1890(1) 

2798.724* O(1), 577(2), 1078(3), 1480(4), 1827(5), 2059(6) 
2809.684 0(1), 546(2), 997(3), 1387(4), 1685(5), 1922(6) 
2863.863* 0(1), 198(2), 365(6), 2080(5), 2252(4), 2362(3) 
2883.642* 0(1), 1886(2) 

2897 .969 —466(6), —259(2), 0(1), 98(4), 256(3), 448(5) 
2926.315* 0(1), 1878(2) 

2938.320 0(1), 58(5), 424(3), 640(2), 798(4), 949(6), 1104(7) 
2960.918* 0(1), 1857(2) 

2988.993 —620(4), —464(3), —258(2), 0(1) 

2993.312 —600(4), —460(3), —255(2), O(1) 

3024.577 —1229(6), —1153(7), —1013(2), —863(8), —642(12), —521(10), —376(4), —290(3), —90(5), O(1), 226(9), 590(11) 
3038.740* 0(1), 1911(2) 

3067 .663 —1000(1), —892(2), —95(2), O(1) 

3076.632 —660(4), —504(3), —277(2), O(1) 

3103.052* 0(1), 1859(2) 

3267 .686* 0(1), 1877(2) 

3397 .302 0(1), 592(2), 1102(3), 1519(4), 1837(5), 2107(6) 
3402.877* 0(1), 588(2), 1089(3), 1501(4), 1833(5), 2075(6) 
3405.326 0(1), 1760(2) 

3510.975 0(1), 568(2), 1058(3), 1465(4), 1782(5), 2026(6) 
3596.123 —1040(6), —573(2), 0(1), 458(4), 841(3), 1320(5) 
3888.166 —179(2), 0(3), 1708(4), 1888(1) 

4121.906 0(3), 714(2), 1879(1), 2597(4) 

+4271.509* 0(1), diff. 285(2), diff. 490(3) 

4308.523 0(1), 1896(2) 

4493.005 0(1), 1850(2) 

4722.373 —1430(2), —1284(3), —1089(5), —459(6), —260(4), O(1) 
16134.570* —1196(4), —910(3), —795(6), —565(5), —384(2), O(1), 110(7) 





t The numbers in parentheses give the qualitative intensities: 1 — strongest, 2 — next strongest, and so on. 


t Observed only in the first order. 
* Structure not investigated before. 


II although their structure could not be resolved 
and only the position of the center of gravity was 
determined. The very good agreement of the 
wave numbers obtained and the term values 
determined from other lines shows that the cor- 
rections for the index of refraction of air given 
by Meggers and Peters and used in the prepa- 
ration of tables by Kayser are correct within the 
limits of accuracy of my measurement. The more 
recently reported formulas for the correction, 
having a smaller coefficient in the term A~™, 
seem to be inadequate for the short wave-length 
region (A2500-A2000).'° The small difference 
between the wave numbers and the term values 


10 Cf, the discussion of this point in Rev. Mod. Phys. 
14, 65 (1942). 


amounting to +0.03 cm™ is opposite in direction 
to the expected shift of the center of gravity 
caused by the h.f.s. and also to the shift required 
by the newer formulas. Moreover, the real dif- 
ference is probably even a little higher, since the 
wave-lengths were determined from the third 
order by use of second order Fe standards; the 
correction introduced for the departure of the 
temperature and pressure from normal was found 
for the lines around \2500 to be rather a little 
too small than too big. 

The h.f.s. collected in Table I gave the h.f.s. 
separations of levels presented in Tables II and 
III. In the last table the less certain values are 
marked by stars. For the ground state the con- 
siderable uncertainty is caused by the incomplete 

















ARC SPECTRUM OF BISMUTH Bi |! 531 


TABLE II. Even terms and their combinations with the 6p* terms. 
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Even 7 

terms 
ty 2 

2s 

35/2 

43/2 

Si/2 

61/2 

7o/ 

83/: 

93/2 
105/2 
19i/2 
116/2(Th) 
125/2 
20:/2( Th) 
213/2? 
22772 
231/2 OF 3/2 
241/2 OF 3/2 
133/2(T) 
145/2(T) 
2517 
153/2(T) 
l6sn(T) 
17s/2(T) 
1832(T) 


Over-all 
h.f.s. 

Separation 
in cm™! 


s. =single. 


6p 


32588.66 
43912.43 
44816.91 
44865.13 
45915.61 
47373.61 
48491 .31 
49461.59 
$1019.18 
51158.47 
§2244.70 
53878.82 


53965.95 


56330.04* 
56569.45 


57075.66 


57606.14* 
§8273.21 
(43208) ? 
(64020) 
(66103) 
This work 


Z. B. G 
(reter 
ence 1) 


4°, 
0 
3067 .663 
32588 .66 


2228.207 
s. 44865.15 


(2110.31) 
47371.3 


(2061.70) 
48488.1 


(2021.21) 
49459.3 


(1959.48) 
§1017.3 


(1953.89) 
51163.2 


(1913.6) 
52241 


(1855.9) 
53865 


(1852.3) 
53970 


(1832.3) 
54558 


(2313.8) 


43205.6 


2D 3/2 
11418.77 
4722.373 
21169.89 
3076.632 
32493.66 


2993.312 
33398.13 


2988.993 
33446.36 


2897 .969 
34496.84 


2780.450 
35954.82 


2696.614 
37072.56 


2627.834 
38042.82 


2524.470 
39600 .34 


2515.619 
39739.69 


2448.683 
40825.90 


(2354.49) 
42459 


(2349.10) 
42556.5 


(2317.70) 
43133.4 


(2251.73) Th 
44396.6 


(2214. 


11) 
45150.7 
(2189.58) 
45656.6 


(2164.1)? 
46194. 


(2133.68) 
46852.3 


(3144.6) 
31791.4 


—0.606 


—0.606 


v.w. =very weak. diff. =diffuse. 
+ Observed only in the first order. 
* New terms and new lines. 


® Ds 
15438.42 


3510.975 
28474.01 


3402.877 
29378.50 


3397 302 
29426.71 


3024.577 
3305 2.88 


2938.320 
34023.15 


2809 .684 
35580.73 


2798.724 
35720.06 


2600.654 
38440.38 


2499.52)? 
39995 .6 


2496.678 
40041.15 


(2430.45) 
41132.1 


2400.965 
41637.13 


(2333.79) 
42835.6 


(3599.9) 
277704 


(2057.68) 
48582.8 


(1973.08) 
50665 .7 


2.082 





2.030 


* 


2P°i/2 
21661.81 


4493.005 
22250.61 


4308 .523 
23203.31 


4121.906 
24253.80 


3888.166 
25711.80 


3596. 


3405 .326 
29357 .37 


3267 .686 
30593.90 


3103.052 * 
32217.02 


(3093.58) T? 
32315.7 


3038.740 . 
32898.83 


2960.918 . 
33763.47 


2926.315 . 
34162.70 


2883.642 * 
34668 .23 


2863.863 
34907 .64 


2781.260 * 
35944.33 


2730.583 
36611.40 


(2360.1) 
42358.1 


(2249.38) 
44442.9 


Over-all h.f.s. 
| Separation in cm 


2P°3/2 Z. B. G. 
33165.01 This work (reference 1) 
0.828 0.830 
0.056 ~O.1 
(8579.74) +0.005 <0.1 
11652.2 
(8544.54) —0.020 ~—0.1 
11700.2 
(7840.33) —0.715 —0.708 
12751.1 
(7036.15) 0.180 0.20 
14208.4 
3.113 3.166 
16134.570 1.430 1.415 
16296.56 
(5599.41) 0.160 very narrow 
17854.1 
0.017 
0.060 
0.025 
| <+0.005 
~ —0.04 
~0.04 
—0.015 
| 0.012 
| 
+4315.700? * | 3S +0.005 
v.w. 23164.72 | 
14271.509 * | —0.477 
23404.37 | 
] 
4181.1002_* —0.930 
v.w. 23910.44 
diff. 
4090.310 | 0.242 
s. 24441.15 
—0.265 


(3239.73) 
30854.8 


(3035.18) 
32937.5 


0.250 
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TABLE III. Hyperfine structure. 
Term Separation in 1073 cm™! 
F=2 3 4 5 6 
| | | 

4S°s/2 — 55* — 75* — 95 
2D° 3/2 —146 -—200 —260 
2D ° sie 245 332 416 500 590 
2P°s/0 60* 80* 110 

75/2 375 500 620 748 870 

83/2 382 478 570 
13s/2 —110 — 169 — 198 
145/2 —118* —154* —190* —226* —242 
153/2 — 60* — 80* —125 








* Relative to the"level“with highest F value. 
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resolution of the line 43067; only the separation 
of the levels F=6 and 5 could be taken directly 
from the data of Table I and the two others were 
evaluated approximately. The value of the 
over-all separation is therefore also uncertain; 
however, it is clear that the total separation is 
not smaller than 0.2 cm~ and that consequently 
Zeeman, Back, and Goudsmit! greatly under- 
estimated its value. Two separations in the level 
15 were evaluated roughly in a similar manner 
from the line 42730. Also for the level 14 only 
one separation of sublevels could be obtained 


TABLE IV. Odd terms (6p2nx) and their combinations with the term 1)/2. 





Line 
Term d v 
1°3(Th) (41125.3) (11711.1) 
8536.6 
2° 3/2 (42941)* (9657) ° 
10352 
ts 50599.67* +5550.621 - 
18011.01 
4° 3/2 53468.70* 14787.928 ” 
20880.04 
5° 53708.49* 14733.566 ° 
21119.83 
6°:2(T)  53893.70 74692.417 
21305.04 
7°s/2 54570.59?* 74547.920 ° 
21981.93 
8° 55393.67* 4383.774 . 
22805.01 
9°(T) 56088.38 4254.175 
23499.72 
10°12 56275.96* 4220.486 ° 
23687 .30 
11° 56316.47* 4213.280 - 
23727.81 
7” 56908.97* 4110.634 ss 
24320.31 
13° 57309.64* 74044.008 ° 
wW. 24720.98 
14° 57593.72* 73998.064 " 
Ww. 25005.06 
iy 57800.77* +3965.229 43 
‘ v.w. 25212.11 
16° 64711.14* 3112.185 “4 





t Observed only in the first order. 
* New terms and new lines 


32122.48 


7 H.f.s. in 10-8 em>! Separation in cm- 
0(1), 817(2) ~0.015 
0(1), 804(2) ~0.040 
0(1), 848(2) ~—0.035 
0(3), 560(1), 830(2), 1390(4) — 0.560 
0(1), 1022(2) ~ —0.250 
both very diffuse 
0(1), 825(2) +0.008 
0(1), 833(2) +0.010 
0(3), 306(2), 830(1), 1134(4) — 0.305 
0(1), 820(2) ~0.012 
0(1), 835(2) ~—0.012 
0(1), 820(2) ~0.012 
0(1), 817(2) ~0.015 
0(1), 870(2) ~—0.06 
0(1), 840(2) ~— 0.020 
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from the line \2401; the others were calculated 
by approximate application of the Landé interval 
rule (the total separation being known from the 
over-all width of the line \2401 minus the over-all 
separation of the level ?D°s,2). Finally the two 
smaller separations in the level ?P°3;2 were ob- 
tained indirectly from only one line (A6134) and 
are therefore of a low accuracy. 

The high accuracy obtained in the case of the 
levels 2D°; 2, 2D° 5/2, 75 2, and 83/2 makes a study 
of interval rule deviations possible. However, 
these deviations have already been studied for 
the line \4722 and for several lines of Bi II and 
Bi III by Schiiler and Schmidt" with great 
precision and the quadrupole moment of the 
nucleus Bi**® was calculated. Since the results of 
the present work for the level 2D°3,2 almost 
coincide with theirs, no new determination of the 
quadrupole moment from the other three levels 
mentioned above has been attempted. 

Of the odd terms corresponding to tue excita- 
tion of at least one electron (6p*nx), only three 
were known previously. Thirteen more were 
detected in the course of this investigation. All 
these terms combine with the even term 1, and 
were easy to identify in view of the wide doublet 
structure of the term 1,. They are all collected 
in Table IV. The term 1° was discovered by 
Thorsen? and was denoted by Bacher and 
Goudsmit® as 20°; the term 2° was identified by 
the author on the basis of combinations. The 
term 7° and the classification of the line \4548 
are subjected to some doubts in view of the large 
deviation of the doublet separation from the 
value 0.828 (term 1,). An _ unclassified line 
\3659.862 was observed with the following 
structure: —1824(2), —1637(3), —1378(5), 
—597(6), —333(4), O(1); the separations of 
levels involved are: a doublet level +1045, and 

1H. Schiller and Th. Schmidt, Zeits. f. Physik 99, 717 
(1936). 
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TABLE V. Lines not involving 6p levels. 





Linet 
ny v Classification 

7335.01* 13629.5 13 —2° 

7838.70 12753.7 11 —1° 

8210.83 12175.7 9°—2 

8761.54* 11410.4 10°—4 

8907.81 11223.0 9°—4 

9828.8 10172.8 9°—5§ 
10106.1 9892.3 9 -1° 
10301.7* 9704.5 7°—4 
11073.2 9028.4 6°—4 
11555.5* 8651.5 4°—3 
11994.5 8334.9 8 -—1° 
12166.5* 8217.1 10 —2° 
22554.2* 4432.6 6°—8 








t All wave-lengths taken from Toshniwal (reference 8). 
* Newly classified. 


a fourfold level #187, 262, +333. However, 
since the components of \4548 are very diffuse 
the doublet separation of the level involved is 
certainly considerably different from the value 
1022 and both these lines certainly do not have 
a level in common. 

In Table V are collected all remaining lines 
classified as transitions not involving one of the 
6p? levels. The small number of such lines 
(Tables IV and V in all 29 lines) indicates that 
our knowledge of the arc spectrum of bismuth is 
still fragmentary. It seems that more substantial 
progress can be expected at present from a study 
of the spectrum in the infra-red and short ultra- 
violet regions only. 

In the course of this work several forbidden 
lines emitted from the four metastable levels 
(configuration 6p*) were observed. A closer study 
of this group of lines is intended in the near 
future. 

The author is very much indebted to Professor 
R. S. Mulliken for his kind interest in this work 
and to Dr. F. Bueso-Sanllehi for the fruitful col- 
laboration in a considerable part of the experi- 
ments. 
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The discharge in fast counters has been localized by numerous devices. The reduction of the 
field in the counter by surface charges accumulated on an insulator has been shown to interrupt 
the spread of the discharge. This implies strong photon absorption in the gas of the fast counter. 
The efficiency of localization has been investigated as a function of the concentration of alcohol 
in an argon-alcohol filling. It is shown that localization is more complete at low values of over- 
voltage. The effectiveness of various gas fillings is given, among which are a number of inorganic 
gases which show both localization and fast counter action. 


REVIOUS studies of fast G-M counter action 
have indicated that the discharge initiated 
by an ionizing particle spreads very rapidly along 
the length of the counter wire. Stever! has re- 
cently shown that the discharge can be con- 
strained by placing a glass bead on the wire. The 
same effect has been accomplished by Ramsey’ 
by separation of the cathode into segments. This 
localization of the discharge has increased the 
versatility of the G-M counter, and a number of 
applications immediately suggest themselves; the 
determination of particle direction and range are 
perhaps the most important. For if the con- 
straining device is completely effective in local- 
izing the discharge, then a primary ionizing 
particle must produce ions in each section of a 
segmented counter which shows a discharge. 
The present paper reports on some of the 
properties of this localization phenomenon. The 
values given apply only to the particular tube 
used by the authors. Since G-M counter prop- 
erties depend upon many factors, the data should 
perhaps be considered a guide as to what may be 
expected with other tubes. 


APPARATUS 


The tube used in most of these experiments was 
designed to permit direct visual observation of 
the wire and to lend itself easily to internal 
changes. The cathode was made of ordinary 
copper wire screen. The inside of this screen 
cylinder was painted with a heavy coat of 
Aquadag and the edges were rolled outward in 
order to reduce the background rate. Other 


1H. G. Stever, Phys. Rev. 61, 38 (1942). 
2W. E. Ramsey, Phys. Rev. 61, 96 (1942). 


details of its construction are shown in Fig. 1, 
Tube A. The position of the quartz disk in the 
figure is typical of other devices used. 

The means used to determine whether or not 
localization took place were (1) direct visual ob- 
servation of the wire, and (2) the pulse sizes on 
the oscilloscope screen. It was possible to observe 
corona flashes along the wire for an individual 
pulse only when the tube had such a filling as to 
yield a slow counter. When operating as a fast 
counter, the corona was visible only when a high 
counting rate was attained. A small gamma-ray 
source held near one end of the tube produced a 
the segment 


seemingly continuous glow in 


nearest it when fairly complete localization was 


TUBE A 


TO VACUUM 
MANIFOLD 














S33 BRASS C2) PYREX 


Fic. 1. Design of counters used to study discharge 
characteristics. Tube A, for study of localization phe- 
nomenon; Tube B, for study of effect of surface charge 
collected on walls of tube. 
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G-M COUNTERS 


TABLE I. Effects of objects on spread of discharge in fast and slow counters. 


No. Object Dimensions 
1 Quartz disk 25-mm diam., 2 mm thick 
2 Plicene disk 20-mm diam., 2 mm thick 
3 Quartz tubing 1 cm long,4 mm i.d. 


6 mm o.d. 


4 Quartz tubing 1 cm long, 1 mm i.d. 
2 mm o.d. 
5 Bakelite ring 8 mm thick, 9 mm i.d. 
30 mm o.d. 
6 Bakelite ring 8 mm thick, 15 mm i.d. 
30 mm o.d. 
7 Bakelite ring 8 mm thick, 19 mm i.d. 
30 mm o.d. 
8 Brass ring 8 mm thick, 19 mm i.d. 
30 mm o.d. 
9 Split cathode l-cm separation 
10 Split cathode 2-cm separation 


taking place. This was helpful in the ambiguous 
case when only one pulse height was observed on 
the oscilloscope screen. It may be interesting to 
note that the glow produced by the gamma-ray 
source in the fast counter was of uniform intensity 
over the part of the wire that is covered. If each 
pulse did not spread the full length of the glow, an 
inverse square law of illumination might have 
been expected. 

The size of the pulses viewed on the oscilloscope 
furnished the more direct means of finding the 
extent of spread of a discharge along the wire. 
This is true because the total charge in a single 
pulse is directly proportional to the effective 
length of the counter. Hence, if a localizing device 
separated the counter into unequal parts, sepa- 
rate pulse sizes identified each segment and a 
third size, approximately equal to the sum of the 
other two, indicated complete spreading. 


TABLE II. Comparison of localization due to a split cathode 
in an ordinary and an extended wall counter. Cathode 
separation, 2 cm. Overvoltage 100 volts. 
Percent localization 
Alcohol vapor 
concentration 


Straight wa! Extended wall 





100% 95 85 
40% 90 75 
20% 85 15 
15% 65 3 
10% 35 3 

5% 20 2 
2% 8 0 


Percent 
localization 








Vapor Pure 
Mounting filling argon 
Through small hole in center 97 0 
Rigidly at center of wire 95 0 
Slung on wire at mid point 95 0 
Slung on wire at mid point 96 0 
Fitted on inside of cathode 90 
Fitted on inside of cathode 87 - 
Fitted on inside of cathode 70 0 
Fitted on inside of cathode 0 0 
Cathode cut perpendicular to 0 0 
its axis at midpoint 
Cathode cut perpendicular to 65 0 


its axis at midpoint 





LOCALIZING DEVICES 

The effects of a number of devices on the 
localization of the discharge are shown in Table I. 
Each was used with a filling of 9 cm argon and 1 
cm ethyl alcohol, or with pure argon. When com- 
plete spreading or complete localization did not 
occur, the percentage of single sized pulses as 
observed on the oscilloscope was used to measure 
the effectiveness of the localizing device. These 
percentages were found to be a critical function 
of the overvoltage, the voltage above threshold 
at which the tube was operated. Therefore in this 
experiment the overvoltage was kept constant at 
140 volts. Table I shows that localization occurs 
only when a fast counter filling is used. The 
previously mentioned corona glow test always 
indicated complete spreading when pure argon 
was used. 

It is possible to explain the localizing action 
in fast counters on the basis of a reduction in field 
strength over a part of the length of the counter, 
even though in some cases the opacity of the 
object to radiation may have been effective. In 
the case of the split cathode, reduction in field 
intensity comes about in an obvious way. When 
any insulator is placed in the counter it acquires 
a surface charge which acts to reduce the field 
intensity in the active part of the counter. A ring 
slipped into the cathode acquires positive charge, 
and a bead on the anode acquires negative charge. 
The brass ring slipped into the cathode, similar 
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Fic, 2. Curves showing effectiveness of a Bakelite ring 
(No. 7, Table I) and a split cathode (No. 10) on localization 
of counter discharge for various mixtures of alcohol and 
argon. 


in dimension to one of the insulating rings, had no 
localizing effect. It may be important for a 
variety of experiments that so simple a device as 
an insulating ring, which does not obstruct the 
body of the counter, serves well to segment the 
tube. 

Even in the case of the split cathode, surface 
charge on the uncovered glass wall plays an im- 
portant role. A tube constructed with an ex- 
tended glass wall (Fig. 1, Tube B) so as to reduce 
the influence of surface charges gave markedly 
different results than the ordinary tube with 
straight walls (Fig. 1, Tube A). Table II shows 
that localization is consistently more complete in 
the ordinary tube for various concentrations of 
ethyl alcohol in argon. These data were taken at 
constant overvoltage. 

Figure 2 compares the effectiveness of a Bake- 
lite ring (No. 7, Table I) and a split cathode in an 
ordinary tube (No. 10) in producing localization. 
The Bakelite ring, whose surface charges are 
nearer the central wire, is the more effective. The 
differences between the curves are in degree 
rather than kind. The particular coordinates 
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chosen lend themselves most easily to reading 
data from the curves. Since pulses of double 
height indicate the spreading of the discharge 
through the reduced field region, it is a low value 
of double to single height pulses that indicates a 
high degree of localization. 

If reduction of field intensity over a part of 
the length of the counter is effective in interrupt- 
ing the step by step mechanism which spreads 
the discharge along the wire, it implies that the 
‘“‘mean free path”’ of the agent responsible for the 
spreading is of the order of magnitude of the 
length of the reduced field region. Photons* are 
considered responsible for the propagation of the 
discharge. They undoubtedly produce abundant 
electrons or ions to build up the discharge until 
it is quenched by the positive ion space charge.‘ 
The localization phenomenon apparently depends 
upon a strong photon absorption by the added 
vapor or gas. Localization fails to occur when 
photons penetrate the reduced field region and 
reestablish the discharge. Since the new ions are 
formed in the gas by a strong photon absorption, 
they are undoubtedly far more abundant than in 
the case of the slow counter. For a certain dis- 


TABLE III. Type of counter action produced by 
different gases and vapors. 











Type of Per- 
Base gas counter cent 
Pressure Vapor or gas added action _locali- 
9 cm Pressure 1 cm Formula produced zation 
argon alcohol C:H;OH fast 95 
argon acetylene C2He fast 93 
argon ethylene C.H, fast 96 
argon methane CH, fast 50 
argon alcohol C.H;OH fast 95 
oxygen Oz 
argon water vapor H:O fast 50 
argon carbon dioxide CO, fast 2 
argon carbon disulfide CS. fast 99 
argon carbon CCl, fast 99 
tetrachloride 
air alcohol C2H;OH fast 90 
oxygen alcohol C.H;OH fast 98 
hydrogen alcohol C.H;OH fast 96 
carbon =: CO, fast 10 
dioxide 
argon oxygen Oz slow 0 
argon bromine Br» slow 0 
argon — A slow 0 
oxygen — Or» slow 0 
hydrogen — H. slow 0 
air — slow 0 


3W. E. Ramsey and E. L. Hudspeth, Phys. Rev. 61, 95 
(1942). 

4C. G. Montgomery and D. D. Phys. 
Rev. 57, 1030 (1940). 


Montgomery, 
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charge current, each new ion or electron produces 
a smaller total number of ions by collision since 
there are more new ions formed in the fast 
counter than in the slow counter. This may 
account for the lower voltage at which the fast 
counter operates and the fact that the space 
charge is able to extinguish it. 


EFFECT OF COUNTER FILLING 
AND OVERVOLTAGE 


The above considerations suggest several prop- 
erties that the localization phenomenon should 
show. Unless the reduced field region is very long 
compared to the photon half-value length, some 
photons should penetrate this region and localiza- 
tion should not occur for every discharge. The 
length of the region reduced below threshold field 
values, and consequently the efficiency of localiza- 
tion, should depend upon the overvoltage. The 
range of the photons should depend upon the 
nature and the concentration of the added vapor 
or gas. 

Figure 3 shows at the same time the effect of 
alcohol concentration and overvoltage. The co- 
ordinates are the same as used in Fig. 2. The 
extent of localization changes extremely rapidly 
at the lower alcohol concentrations. At the higher 
concentrations it is nearly 100 percent and de- 
pends only very slightly upon concentration. It 
is also apparent that the localizer is most effective 
at low overvoltage. 

Table III ennumerates the results obtained 
with a variety of counter filling combinations. 
These data were taken at voltages representing 
the middle of the plateau in each case, and the 
localizer used was the Bakelite ring (No. 7, 
Table I). The criterion for a fast counter was the 
existence of a usable plateau with one megohm in 
series with the counter. It is interesting to note 
that several inorganic gases show fast counter 
action. The experiment with carbon dioxide indi- 
cates that fast counters are obtainable even 
though there is only a small tendency to produce 
localization. However, localization was observed 
in all the fast counters. The shortest plateaus 
were about 40 volts with argon-water and argon- 


vi 
w 
~ 


G-M COUNTERS 



























































I 
3 q 
” fl! 
4 0 
Fe | 
ql 
WwW 
yt 
o i 
z 
a " 
" 
2 rT OVERVOLTAGE 
wm <tk! * 50 VOLTS 
om oN e 100 VOLTS 
3 nl e 150 VOLTS 
ei 
| 
Sf 
oO 
kb 
x 
_” = 
° 
o* 6 
° v 
al ° 
-3 
0 20 40 


PERCENT OF ALCOHOL VAPOR 


Fic. 3. Dependence of localization on the overvoltage 
at which the counter is operated. 


carbon dioxide fillings, and about 60 volts with 
the argon-carbon disulphide filling. 

Of the gases used, only monatomic and dia- 
tomic gases and their mixtures showed slow 
counter action. Gases whose molecules contain 
three or more atoms and at least two different 
atomic species gave fast counters. For example, 
air, pure argon, and argon-oxygen mixtures gave 
slow counters, while pure alcohol, carbon dioxide, 
and mixtures of these with argon gave fast 
counters. Bromine showed slow counter action, so 
molecular weight alone is apparently not suffi- 
cient to produce a fast counter. Presumably the 
larger number of degrees of freedom and the lack 
of atomic symmetry in the polyatomic molecules 
permit more abundant rotational and vibrational 
levels in these molecules and may account for 
their greater photon absorption. However, at- 
tempts to correlate any detailed photo-ionization 
mechanism with our results were completely 
unsuccessful. 
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The Wave Function for the Negative Hydrogen Ion 


RALPH E. WILLIAMSON 
Yerkes Observatory of the University of Chicago, Williams Bay, Wisconsin 
(Received October 23, 1942) 


A six-parameter wave function for the negative hydrogen ion has been derived, and several 


tests for the goodness of a wave function have been applied to it. 


N connection with some astrophysical prob- 
lems,! the writer has recently used the meth- 
ods of Hylleraas to derive an accurate wave 
function for the negative hydrogen ion. The trial 
wave function assumed was of the same form as 
that used by Hylleraas? in predicting the energy 
of the ground-state of helium within the limits 
of observational error. It is found that if 71, 72, 
and rj. denote, respectively, the distances of the 
two electrons from the nucleus and from one 
another, then 
v= Nexp (—a(ri+r2))[1+Bristy(r2—11)” 
+ 6(rotri)+e(ret+ri)?*+irie], (1) 
with 
N =0.073965, 
a=+0.70120, 56= —0.14387, 
B= -+0.26017, e= + 0.00436, (2) 
y=+0.07225, ¢=—0.00896, 


provides a “‘sixth’’ approximation to the wave 
function. The wave function currently used for 
H~ is that of Bethe* and Hylleraas:4 


V=WN exp (—a(rit+re)) 
X(1+6riety(re—r1)?], (3) 
where 
B= +0.308, 
y= +0.119. (4) 


N=0.0623, 
a=+0.770, 


All quantities in Eqs. (2) and (4) are expressed in 


atomic units.® 
On the basis of the new wave function, the 





'R. E. Williamson, Astrophys. J. 96, 438 (1942). 
* E. A. Hylleraas, Zeits. f. Physik 54, 347 (1929). 
3H. Bethe, Zeits. f. Physik 57, 815 (1929). 

*E. A. Hylleraas, Zeits. f. Physik 60, 624 (1930). 

‘®In an earlier paper by P. Starodubrowsky (Zeits. f. 
Physik 65, 806 (1930)), an attempt has been made to im- 
prove the Bethe-Hylleraas wave function. However, 
Starodubrowsky estimated only the decrease in the energy 
of the ground state due to the successive addition of various 
terms to the wave function, and did not explicitly obtain 
the wave function itself; also his final estimate of E is 
—0.52628, which is to be compared with E= —0.52646 
derived from (1). (Cf. Eq. (5).) 


energy E of the ground-state of the H™ ion is 


found to be 


E= —0.52646 atomic unit. (5) 


This value is even slightly lower than the mini- 
mum value of E= —0.52642 found by extrapola- 
tion from other two-electron systems.® Since the 
difference between the two is smaller than the 
correction for motion of the center of mass, and 
of the same order as the relativity-correction,’ it 
seems that the wave function (1) represents the 
ground state of H~ almost within the limits of 
accuracy of the Schrédinger equation itself. We 
can safely put for e, the difference between the 
true and the computed values for the energy 
level, 
e~10- atomic unit. 

The accuracy of the new wave function there- 
fore appears to be quite high. It is consequently 
of some interest to apply to this wave function 
the alternative tests of accuracy which have been 
suggested ,* with a view to discriminating between 
them. 

Now, the smallness of the quantity 6°, where 


= [ (- Ev ar, (6) 


and y is the trial wave function, has been advo- 


TABLE I. Tests of accuracy applied to wave 
functions for H~. 


Wave 

function e & €/& E 
(1) 10-5 0.00302 0.004 — 0.52646 
(3) 0.0011 0.00316 0.4 — 0.52531 


6 Hylleraas, Zeits. f. Physik 60, 624 (1930). 

7H. Bethe, Handbuch der Physik (Berlin, 
XXIV, Part 1, p. 363. 

8D. H. Weinstein, Proc. Nat. Acad. Sci. 20, 529 (1934); 
J. H. Bartlett, Jr., J. J. Gibbons, Jr., and C. G. Dunn, 
Phys. Rev. 47, 679 (1935); H. M. James and A. S. Coolidge, 
Phys. Rev. 51, 860 (1937); A. A. Frost, J. Chem. Phys. 10, 


240 (1942). 
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cated by Weinstein,’ Bartlett,’ and more re- 
cently by Frost® as a criterion of the goodness of 
the wave function. On the other hand, James and 
Coolidge’ indicate that this is not a very reliable 
test of the over-all accuracy of the wave function. 
It is now possible to confirm the conclusions of 
James and Coolidge on the basis of the wave 
functions for H~ now available. Table I gives 
the values for e, &, €/5°, and E for the two wave 
functions (1) and (3). It is seen that the difference 
between the exact and the approximate wave 
functions is not indicated by any significant or 
appreciable change in &. The quantity ¢€/5? shows 
a marked decrease for the exact function, in ac- 
cordance with the expectations of James and 
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Coolidge.* However, in this case, the decrease in 
«/* is almost entirely due to the decrease in e«, 
since 6 changes by only about 4 percent. The 
quantity ¢€/& will show different variations than 
those of « when we are dealing with compara- 
tively inaccurate wave functions; for, then 6? may 
also change appreciably. However when the 
highest accuracy is demanded, the Ritz principle 
provides the simplest and the best criterion. And 
in all events the smallness of the quantity & 
appears to be a very poor criterion for the general 
over-all accuracy of a wave function. 

The writer is deeply grateful to Dr. S. Chan- 
drasekhar for his advice and encouragement in 
the course of the work. 
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Diffusion Rate of Carbon in Iron-Cobalt Alloys 


R. SMOLUCHOWSKI 
Research Laboratory, General Electric Company, Schenectady, New York 
(Received October 26, 1942) 


Studies of rates of various transformations and reactions in iron alloys seem to indicate that 
certain elements have a definite influence on the rate of diffusion of carbon in y-iron. In par- 
ticular, cobalt shows a pronounced accelerating influence. By direct measurement of the rate 


of diffusion of carbon in iron-cobalt alloys this conclusion is confirmed and an increase in the 
diffusion coefficient established. An addition of 4 percent cobalt decreases the activation 
energy from 32,500 cal. down to 30,000 cal. corresponding to an almost twofold increase of the 
diffusion rate at 1 atomic percent carbon concentration. At higher carbon contents the ac- 
celerating influence of cobalt is smaller. This property of cobalt is unusual especially in view 
of the fact that both manganese and nickel are known to have practically no influence on the 


diffusion of carbon. 


NE of the fundamental mechanisms of vari- 

ous reactions occurring in the solid state is 
diffusion. Although its importance, especially in 
applied physics of metals, was early recognized, 
there are few studies based on the modern theory 
of solids. Only recently Huntington and Seitz! 
have shown that a diffusion of a vacancy is the 
most probable mechanism in the case of self- 
diffusion of copper and gives a reasonable value 
for the activation energy. Less is known about 
interstitial diffusion in which the diffusing atom 
does not replace atoms of the solvent but pene- 
trates between them. 


1H. B. Huntington and F. Seitz, Phys. Rev. 61, 315 
(1942). 


A most widely studied example of interstitial 
diffusion is diffusion of carbon in iron,?~ which 
is of great practical importance. The rate of 
diffusion seems to have an important role in de- 
termining the rate of various transformations in 
the crystal lattice. One of these is the decomposi- 
tion of the solid solution of carbon in y-iron. 
Carbon is much more soluble in the face-centered 
y-iron than in the body-centered a-iron. On 


21. Runge, Zeits. f. anorg. allgem. Chemie 115, 293 
(1921); G. Tammann and K. Schonert, Stahl und Eisen 
42, 654 (1922); A. Bramley and A. J. Jinkings, Carnegie 
Schol. Mem. 15, 127 (1926). 

3M. Paschke and A. Hauttmann, Arch. Eisenh. 9, 305 
(1935). 

*C. Wells and R. F. Mehl, Trans. A. I. M. E. 140, 279 
(1940). 
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cooling, y-iron transforms into a-iron and ejects 
excess carbon in the form of the carbide Fe;C. 
We will not discuss here the details of this trans- 
formation which is not yet fully understood al- 
though great progress has been made toward this 
goal.® An interesting phenomenon is the influence 
of other elements on the rate of this transforma- 
tion. Small additions of Cr, Mo, W, Ni, Mn, etc., 
have an appreciable retarding influence. The only 
known exception is cobalt which increases the 
rate. It has been thought that the retarding ele- 
ments might decrease the rate of diffusion of 
carbon thus slowing down the transformation. 
Wells and Mehl‘ have studied directly the in- 
fluence of nickel and manganese on the diffusion 
rate of carbon in y-iron. They found almost no 
influence up to several percent of added element 
and a slight increase of the rate at higher concen- 
tration (see Fig. 2). Thus the mechanism of the 
influence of these two elements on the ya 
+Fe;C reaction must be looked for in another 
direction. There remains, however, the very in- 
teresting case of cobalt which has a pronounced 
and unique accelerating influence on the trans- 
formation described above. This influence seems 
to be related to the increase in the depth of 
decarburization of iron-cobalt alloys with in- 
creasing cobalt content.® Cobalt is in many ways 
an exceptional element in iron alloys. Besides the 
well-known rise in Curie point, it has the unique 
property of leaving the lattice constant of iron 
almost unchanged up to a few percent concen- 
tration, in spite of the difference in atomic radii. 
It has also practically no influence on the tem- 
perature of the a—y transformation up to several 
percent. The aim of the experiments described in 
this paper was to make a direct and quantitative 
study of the presumed exceptional influence of 
cobalt on the rate of diffusion of carbon in y-iron. 


EXPERIMENTS 


The experimental technique was that of weld- 
ing together two rods, one with, the other without 
carbon, and analyzing the layers parallel to the 
interface after an appropriate heat treatment. 
The same method was used by Paschke and 


5 F.C. Hull and R. F. Mehl, Trans. A. S. M. E. 30, 381 
(1942) and earlier papers there listed. 

6 E. Houdremont and H. Schrader, Arch. Ejisenh. 5, 
523 (1932). 


Hauttmann’ and by Wells and Mehl‘ and seems 
to give the best results. 

Alloys were prepared in an induction furnace. 
They contained 0, 1.98, and 3.91 atomic percent 
cobalt. Other constituents were silicon (0.3 
atomic or 0.15 weight percent), manganese (0.2 
atomic) and carbon (0.5 atomic or 0.1 wt. per- 
cent) which were about the same in all the alloys. 
Source of carbon was a carbon steel (5.5 atomic 
or 1.2 wt. percent) with approximately the same 
impurities. The specimens were about 2 cm in 
diameter and 4 cm long. The surfaces to be 
welded were ground very carefully perpendicular 
to the axis of the specimen. A few samples were 
polished after grinding, but it turned out that 
polishing gave surfaces less uniformly flat than 
can be obtained by careful grinding. This is 
especially true near the edges which are fre- 
quently rounded off by polishing. The welding of 
ground surfaces was also much better and 
quicker. Welding was done in air in a standard 
spot-welding machine using about 5000 amperes 
under 80 kg per cm? pressure for a few seconds. 
The short time involved in this process allowed 
us to disregard the decarburization and oxidation 
which would result in a longer procedure. The 
fact that the total carbon content did not change 
during the welding process was checked by 
analysis of the immediate neighborhood of the 
weld before heat treatment. The small amount of 
diffusion which occurs during welding is entirely 
negligible compared with that obtained in the 
subsequent treatment.‘ The welded specimens 
were heat treated for 72 hours at 1000° or for 63 
hours at 1190°C. To prevent any change in 
carbon content the heat treating was carried out 
in argon. Then the outside layer of about 0.25 
mm was removed and thin layers parallel to the 
weld turned off and analyzed. The layers were 
alternately 1.2 mm and 0.2 mm thick. The posi- 
tion of the layer was measured with a microm- 
eter with an accuracy of about 0.05 mm. Nor- 
mally the concentration distribution curves were 
obtained from the analysis of the thicker layers. 
However, if there was any doubtful point or a 
large gradient of concentration, then the inter- 
mediate thin layers were analyzed by a micro- 
analytical method and in this way additional 
points obtained. This scheme proved to be very 
useful and labor-saving. The micro-analytical 
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Fic. 1. Carbon concentration 
curve in iron obtained after 72 
hours at 1000°C. 


method was an adaptation of the conventional 
micro-carbon-hydrogen technique. Samples of 50 
to 100 mg were burned in a furnace connected 
to a standard combustion train.’ 

A typical concentration distribution curve is 
shown in Fig. 1 in which the carbon concentra- 
tion drops from 5.5 atomic percent down to 0.5 
atomic percent within a distance of about 25 mm. 
Only values obtained from the thick layers are 
indicated on the diagram. 


COMPUTATION OF THE DIFFUSION 
COEFFICIENTS 


From a plot of concentration versus distance 
the diffusion coefficient can be calculated by 
means of the Grube*® or Boltzmann-Matano® 
method. These two methods permit the calcula- 
tion of the diffusion coefficient at a given concen- 
tration of the diffusing element. This is important 
since the diffusion coefficient of carbon depends 
on the concentration ;' e.g., at 1000°C it varies 
from 2.6X10~-7 cm? per sec. at 0.45 atomic per- 
cent C to 4.4X10-7 cm? per sec. at 4.5 atomic 
percent C. Since the differences between values 
obtained with the two methods were small, only 

7A detailed description of this method will be published 
shortly in Anal. Ed. Ind. Eng. Chem. by E. W. Balis, H. 
Liebhafsky, and E. H. Winslow, to whom the author is 
indebted for carrying out these analyses. 

Pos Grube and A. Jedele, Zeits. f. Metallkunde 19, 438 


°L. Boltzmann, Wied. Ann. 53, 959 (1894); C. Matano, 
Jap. J. Phys. 8, 109 (1933); 9, 41 (1934). 
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the Grube method, which is simpler, was finally 
used, although the other one is theoretically more 
correct. 

In calculating the diffusion coefficient, account 
had to be taken of the fact that cobalt was pres- 
ent only on the low carbon side of the weld. Thus 
the solvents on both sides of the weld were not 
identical. In order to show how this fact affects 
our calculations, let us consider first the simpler 
case when both solvents are the same. The Grube 
procedure is then as follows: Let Co be the origi- 
nal concentration of the diffusing atoms on one 
side of the interface (see Fig. 1) and C the con- 
centration at a distance x from the interface, at 
time ¢. Let us assume also that the original con- 
centration on the other side of the weld was zero. 
Then the diffusion coefficient corresponding to 
that concentration is given by 


D=x"/4a’*t, (1) 

in which “‘a”’ satisfies the relation, 
&(a) = (Cyo—2C)/Co, (2) 
and # indicates the Gauss error function. The 


concentration at the interface, C;, is then equal 
to Co/2, which is the condition of the applica- 


bility of the Grube method. In other words, the 
diffusion coefficient at the interface must be the 


same, independently from which side we ap- 
proach. If the two solvents are not identical and 


have in general different diffusion coefficients for 
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TABLE I. Diffusion coefficient of carbon in iron-cobalt 
alloys at 1 atomic percent carbon concentration. 
DX 107 cm? sec.7!. 
3.91% Co 


T 0% Co 1.98% Co 





1273° abs) =: 3.15 4.5 6.4 Ist method 
(4.3) 4.5 (6.0) 6.7 2nd method* 
4.5 6.55 mean 
1463° abs 16.4 22.8 29.8 ist method 
(21.8) 23.1 (26.7) 30.9 2nd method* 
22.95 30.35 mean 


* Parentheses indicate values before correction, see text. 


the same concentration of the solute, then there 
is a discontinuity of the diffusion coefficient at 
the interface and C;#}Cpo. If the distribution 
ratio x of the solute in the two solvents is differ- 
ent from unity, then there is also a discontinuity 
in the concentration itself at the interface. For 
a diffusion coefficient independent of concentra- 
tion in a two-solvent system the concentration 
curves are given by :'° 


C,/Co=1-— {«(Dz2)*/(«(D2)'+(D,)')} 

X {1+(x%/2(D,t)')}, 
C2, Co= { «(D,)*/(«(D2)'+(D,)4)} 

X {1—(x«/2(Dst)*)}, 


where Cp is the original concentration in solvent 
I and C, and C; are the concentrations in solvents 
I and II, respectively ; D; and Dz are the corre- 
sponding diffusion coefficients and «x the distribu- 
tion ratio of the solute at equilibrium. We shall 
set x equal to unity. One reason for this assump- 
tion is that there is no discontinuity in concentra- 
tion observed at the interface. Another reason is 
the result of a direct experimental study of the 
equilibrium of carbon in an Fe+CFe+Co+C 
system. A sample with 4 percent cobalt and 
another without cobalt, both containing 3.6 
atomic (0.80 weight) percent carbon, were welded 
together and kept for several days at 1000°C. No 
change in the carbon distribution occurred. It 
was, therefore, concluded that within the limits 
of experimental error carbon is in equilibrium at 
least at temperatures and times here involved. 

If we wish to calculate the diffusion coefficients 
in a two-solvent system we have to make certain 
changes in the Grube calculation and derive some 
new formulae. First of all, from any one of Eqs. 
(3) the following expression for the ratio of the 


(3) 


ie W. Jost, Diffusion und Chemische Reaktion in festen 
Stoffen (Steinkopff, Dresden, 1937); R. M. Barrer, Diffu- 
sion in and Through Solids (Macmillan, New York, 1941). 
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diffusion coefficients in the two metallic solvents 
at the interface (i.e., for x =0) is easily obtained: 
D,(C;) Co—Ci\? 
ma8), 
D,(C;) Ga 
where C; is the concentration at the interface. 
Computation of the diffusion coefficients for x =0 
can be done in two ways. The first method con- 
sists in substituting in (2) the value 2C; for Cp. 
This is based on the fact that the concentration 
distribution on one side of the interface in a 
given solvent does not depend upon the solvent 
on the other side, provided C; is always the same. 
In other words the fact that D, is different from 
D, appears only in a linear factor in Eqs. (3). 
In this way this case is brought back to the usual 
case of one solvent on both sides of the weld, and 
the computation proceeds according to the Grube 
or Matano method mentioned above. The second 
method of making the necessary correction is to 
compute the change in the concentration curve 
brought about by the difference between the two 
solvents. For a given value of x the ratio of the 
concentrations for D,;#D»2 and D,=Dz is easily 
found from (3): 





C2(D,= Dz) 
F,=— — =4/(1 —(D, ‘D;)* ] 
C2(Di+ D2) 
for x>0O, 
or (5) 
Ci(D,:=D2)—Co 
Fe — = }({1+(D./D:2)'] 
C,(D, + D2) —Co 
for x<0. 


Formula (5) cannot be applied directly since in 
our case D, is not known. One has to apply suc- 
cessive approximations: Suppose that in solvent 
II, at a distance x’ from the interface, the concen- 
tration is C; then by applying the normal Grube 
formulae (1) and (2) a certain coefficient D»’ is 
obtained. Knowing from other experiments the 
diffusion coefficient D,, at the same concentra- 
tion, we compute the correction factor F2’ from 
formula (5). Multiplying all measured values of 


TABLE II. Diffusion coefficient of carbon in iron-cobalt 
alloys at 2.2 atomic percent carbon concentration. 
DX 107 cm? sec.*!. 


0% Co 1.98% Co 


1273° abs r 
1463° abs 21.9 


3.7 
18.8 
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C in solvent II by this factor gives a new concen- 
tration curve which corresponds approximately 
to the distribution in the case when both solvents 
are identical. The distance x’ at the concentration 
C is shifted now to x” and another value D.” ob- 
tained. A new factor F” could be introduced, but 
in most cases the process is convergent rapidly 
enough so that only one correction is necessary. 

Finally, we have to consider the case x=0 in a 
one-solvent system where formula (1) is not ap- 
plicable. The diffusion coefficient can be calcu- 
lated in that case from the concentration gra- 
dient dc/dx. From the concentration curve for a 
system with one solvent: 

Co 
eet” [1 —$(x/2(Dt)')], 


we readily obtain the necessary relation 


dc\* 
D..0= (1 Ani)( C./ ) ‘ (6) 
Ox 


This formula has the same theoretical limitations 
for a concentration dependent diffusion coeff- 
cient as the Grube method. 


EXPERIMENTAL RESULTS AND DISCUSSION 


In connection with the application of our 
formulae to actual experimental data one point 
needs mentioning. It has been assumed that the 
concentration of cobalt remained constant during 
the heat treatment. This is not exactly true be- 
cause not only carbon but also cobalt diffuses 
across the weld. There are no direct measure- 
ments known of the diffusion coefficient for co- 
balt in iron. However, we can assume that it is 
of the same order as for nickel and manganese 
which are known." Putting D=2X10-" cm? 
sec.-! we can compute that within 0.2 mm from 
the weld, the concentration of cobalt did not 
change by more than 0.1 percent. This is entirely 
negligible. 

In Table I are given the calculated diffusion 
coefficients for the various alloys at 1 atomic 
percent carbon concentration, as computed by 
the Grube method, and both methods of correct- 
ing for the absence of cobalt on the high carbon 
side of the weld. The first conclusion to be drawn 
from these data is the pronounced influence of 
cobalt on the diffusion rate of carbon in iron 


329 (1941). 
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Fic. 2. (Left) Diffusion coefficients of carbon in iron- 
cobalt alloys, Full and dotted lines correspond to 1 and 
2.2 atomic percent carbon, respectively. The line for Mn 
and Ni was obtained from measurements of Wells and 
Mehl (reference 4). Open circles indicate values calculated 
from Dushman-Langmuir equation. 

Fic. 3. (Right) Activation energies of the diffusion of 
carbon in iron-cobalt alloys. Open circles indicate values 
calculated from Dushman-Langmuir equation. 


(see Fig. 2). An addition of about 4 percent cobalt 
increases the diffusion coefficient twofold. The 
influence seems to be less pronounced at higher 
temperatures. This accelerating influence of co- 
balt is in accord with the qualitative conclusions 
of other experiments on the rates of transforma- 
tion and decarburization in the iron-cobalt- 
carbon alloys discussed above. 

Table II gives the diffusion coefficients as 
computed from the concentration of carbon 
(about 2.2 atomic percent) at the interface, ac- 
cording to Eq. (4). The values for alloys without 
cobalt were computed from the concentration 
gradient according to Eq. (6). Here the influence 
of cobalt is much less pronounced than at 1 
atomic percent carbon. The influence is so much 
smaller that the dependence of the diffusion 
coefficient of carbon on carbon concentration is 
changed. In pure iron-carbon alloys the diffusion 
coefficient increases with increasing carbon con- 
tent. In iron-cobalt alloys the situation changes 
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and for cobalt contents higher than about 2 per- 
cent an increase in carbon content decreases the 
rate of diffusion. This is clear from Fig. 2. Wells 
and Mehl’s experiments‘ show that the influence 
of nickel and manganese on the diffusion rate 
increases with carbon content. Thus also in this 
respect cobalt differs from these two elements. 

The dependence of the diffusion coefficient on 
temperature can be expressed in the known way 
by the Arrhenius equation 


D=A exp (—Q/RT), (7) 


where Q is the activation energy, R the gas con- 
stant, and A, a constant. This constant according 
to Dushman and Langmuir” can be expressed 


in the form 
A=Q8/Nh, (8) 


where WN is the Avogadro number, h is Planck’s 
constant, and 6 is the length of the “‘free path” 
in the atomic picture of diffusion mechanism. 

By use of the mean values of the diffusion 
coefficients obtained at 1000°C and 1190°C (see 
Table I), the activation energy at 1 atomic per- 
cent carbon was calculated from Eq. (7). The 
results are shown in Fig. 3 by full dots. The 
activation energy appears to decrease with in- 
creasing cobalt content, which agrees with the 
observed increase in the rate of diffusion. 

The experimental data can be analyzed also 
from the point of view of the Dushman-Langmuir 
relation, which is known to show a remarkably 
good agreement with experiment, in spite of lack 
of strict theoretical basis. Especially the interpre- 
tation of the magnitude 4 is difficult in practical 
applications. Probably the best use of Eq. (8) is 
to consider A as proportional to Q. However, it 
turns out, that with 6 assumed equal about 1.8A, 
i.e., half of the lattice constant a, of iron, Q 
agrees with experiment. This is the only justifica- 
tion for this assumption since the shortest dis- 
tance between atoms or between equivalent 
interstitial vacancies is about 2.5A, which gives 
Q too large by about 2000 cal. This seems to be 
outside of the experimental error. On the whole, 
Q is very insensitive to the exact value of 4, 
which is one reason for the success of the formula 
(8). For the same reason, it is not necessary to 
take into account the change in lattice constant 
with composition. This is especially true in the 


2S. Dushman and I. Langmuir, Phys. Rev. 20, 113 (1922). 
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case of cobalt which leaves the iron lattice (a-iron 
at least) almost unchanged at the concentrations 
here involved. The Q values calculated from the 
Dushman-Langmuir relation (with 6=4a,) are 
indicated in Fig. 3 by open circles. They are 
based on the value of the diffusion coefficient at 
1000°C, and fall remarkably close to the experi- 
mental values although the arbitrary choice of 6 
has to be kept in mind. With the Q values from 
Eq. (8) the diffusion coefficients at 1190°C can 
be computed. These are indicated on Fig. 2 by 
open circles. The same decrease in Q with in- 
creasing cobalt content can be obtained from 
measurements at 1000°C alone. Making the 
plausible assumption™ that A in Eq. (7) does not 
change much, we have 


D,/D2=exp (—AQ/RT), (9) 


where D, and D; are the diffusion coefficients in 
two different solvents. For an increase in cobalt 
content from 0 to 4 percent Eq. (9) gives a de- 
crease in Q of 1850 cal., in sufficient agreement 
with data plotted on Fig. 3. 

Similar considerations and calculations of the 
Q values for diffusion coefficients can be applied 
to the data for 2.2 atomic percent carbon. How- 
ever, the lower accuracy of these coefficients 
(Table II) does not permit us to draw many con- 
clusions. With the exception of the high cobalt 
alloy for which the diffusion coefficient at 1190°C 
is too high, probably due to an error in analysis, 
the Q values show a small decrease (about 500 
cal.) with increasing cobalt content. This lies, 
however, within the limits of the experimental 
error. 

The decrease in the influence of cobalt with 
increasing carbon content is interesting. The 
point is that, so far as we know, cobalt leaves the 
iron lattice practically unchanged. This and other 
“exceptional” properties of cobalt seem to indi- 
cate that the facts here observed may have an 
“electronic” rather than a purely “‘crystallo- 
graphic” origin. A discussion of the probable 
mechanism of the peculiar influence of cobalt on 
the diffusion of carbon, should be postponed until 
data for other elements are known. 

The author wishes to express his sincere thanks 
to Mr. E. R. Parker for frequent stimulating 
discussions and advice. 


18 The author is indebted to Dr. S. 
remark. 
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The Polarizability of H- 


Louis R. HENRICH 
Yerkes Observatory of the University of Chicago, Williams Bay, Wisconsin 


(Received October 23, 1942) 


The value for the polarizability of H~ is computed from both third and sixth order Hylleraas 
type wave functions. It is found to be 11.58 X 10-* and 14.63 x 10-, respectively, in c.g.s. units. 


HE polarizabilities of helium and ionized 
lithium in an electric field have been com- 
puted by Hassé! and Baber using the variation 
method with unperturbed wave functions of the 
Hylleraas? type. Essentially the same method is 
applied here to compute the polarizability of H-. 
The equation in which we must make the 
energy a minimum is 


Ef var= [&, (grad, Wytdr+2 f Vora (1) 


In this equation 
E=Eote; P=W(l+e); V=VotV’; (2) 


where 
Vo = —Z ‘ty —-Z/tet 1/rie 


; (3) 
V’= — F(x,+-%2). 


The quantities are in terms of atomic units, with 
the energy unit being the ionization energy of 
hydrogen from the ground state. The distances of 
the two electrons from the nucleus are 7; and 72, 
respectively ; x, and x2, etc., refer to the Cartesian 
coordinates of the respective electrons measured 
from the nucleus, and 7,2 is the separation of the 
two electrons. Ey is the energy of the unperturbed 
ground state ; € is the change in energy due to the 
presence of the field, F. Yo is the unperturbed 
wave function for the ground state and (1+ ¢) 
represents the factor by which the unperturbed 
wave function is to be multiplied in order to ob- 
tain the perturbed wave function. The function ¢ 
is assumed to be of the form 


g=A(x1+%x2) + B(xy714+Xefe). (4) 
1H. R. Hassé, Proc. Camb. Phil. Soc. 26, 542 (1930); 
27, 66 (1931); T. D. H. Baber and H. R. Hassé, Proc. 
Camb. Phil. Soc. 33, 253 (1937). 
2 E. A. Hylleraas, Zeits. f. Physik 54, 347 (1929). 


Introducing Eq. (2) in Eq. (1) we have 


Ef wi(+o)%dr 
= [+e Dd: (grad; Yo)2dr 
+ fv > ; (grad; ¢)*dr 
+2 (1 + )Wo ¥; (grad; g- grad; po)dr 


+2f Volt toydrt2f yet 9dr. 
(5 


In solving for the unperturbed wave function 
the solution was made with the condition that 


Eo { wotdr= { ¥. (grad; vo)*dr +2 f VoWo*dr. (6) 
In addition we find by integrating by parts that 
fe >; (grad; ~o)*dr 
+ fv >: (grad; po: grad; g")dr 
= = [doe > i Awodr. (7) 


Now introducing Eqs. (6) and (7) into Eq. (5) we 


find that 
cf 1420+ evetdr= — fval2e+e% 
X CX AwWot+(E—2Vo)po jdt 
+ [ve >; (grad; ¢)*dr 
+2f Vvot(1+2e+6°)dr. (8) 


Now y> is not an exact solution of the unper- 
turbed wave equation. However, with a third- or 
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546 I R. HENRICH 
TABLE I. Constants and polarizabilities. 
He H- 

Units Order a A/F B/F a A/F B/F 
Atomic 3rd 2.62 0.2859 0.3657 156.3 1.679 1.151 
Atomic 6th 2.72 0.2888 0.3727 197.5 1.418 1.340 
C.g.S. 3rd 0.1942 x 10-** _— 11.58 10-* —- ~~ 
C.g.S. 6th 0.2016 x 10-** — — 14.63 X 10-*4 ~- -- 





sixth-order Hylleraas type wave function, the 
error will probably not be large if we set the first 
integral on the right-hand side in Eq. (8) equal to 
zero. 

In addition the Hylleraas type wave function 
is expressible as Yo(r1, 72, 712). Consequently, in 
terms of Cartesian coordinates we have 


Wo(x1, Xo, +++, 2) =Wo(—xX1, —X2, +++, —Z2). (9) 


The functions g and V’ contain factors of x; or 
x2 so that we may expect certain integrals to 
vanish in the present case. In particular, 


f evtar=0; [ Viverdr=0; J Viveear=o. 
. (10) 


Also, for small values of the field strength F we 
may neglect the term ¢/ ¢’o"dr. 
Consequently, Eq. (8) becomes 


, f do'dr= f vo? Es (grad; g)'dr-+4 f Vebidr. 
( 


On performing the differentiations, etc., we 
find that 


>: (grad; ¢)*dr=2A? 


x xe? 
+24B(ntre+ —-— 
YT) To 
+B? {r2+re2+3(xyr2+x27)}, (12) 
and 
— V'g=AF(x1?+2x1x2+%x2") 
+BF}xy?r714+x1%0(ritre) +x2"r2}. (13) 


Then it is found convenient first to express the 
integrands in terms of the variables (71, 72, 712). 
Since the unperturbed function yo is a function 
only of these variables, we see from symmetry 


that 

Jf eenetdr=a fet tystetastyvotdr, (14) 
and 

J sxobear= bf Gouvetyiyetsee)votdr (15) 
Hence, we find 


five >i (grad; o)'dr= fv? (grad f)*dr, (16) 


where 


(grad f)?=2A?+ Ee I OAT, 


Also, we find 


af V'ebotdr=4 f Pyotdr, 


(18) 


where 


AF , 
—— (2(r1? +127) — rie} 


BF 2 
bint trb ta (rit re) (rita? riz) }. (19) 


We now introduce the Hylleraas variables; i.e., 


we let 


S=rntnre; t=fe-—71; U=N}2 (20) 
Using these we have 
(grad f)?=2A?+ 8A Bs+ B*(s*+??), (21) 


and 
—4P= 464A F(s?+??— 1’) 
+2BF(s*+2st?—su*). (22) 





‘F 


51 
40 


(14) 


(15) 


19) 


.e., 


20) 


21) 





THE POLARIZABILITY OF hk 


The Hylleraas type wave function used is 
Yo=e-™(1+cu+dt?+es+fs?+gu*). (23) 


We see that all of the functions that we need to 
integrate are symmetrical in r; and r.—or involve 
only even powers of ¢—so that, if the general 
integrand is g(r, ’2, 712), we have 


fet re, ruddr= at f asf du 


xf G(s, t, u)u(s?—f)dt. (24) 


0 


These can all be broken down into terms of the 
following type {(f-fe"*s*t’u'dsdidu, where 


f f f e~2"*s2p>ycdsdidu 


1 i (a+b+c+2)! 7 


(b+1)(b+c+2) 


(2n)statd+e 


1 
=—___{¢, d, ¢] (25) 


(2n)' 3+a+b+c 


Then, if we express the volume factor separately, 


we have 
ae ff femsetur-u(s?—edsdida 
2r? 
=—_—_—— {a,b,c}, (26) 
(22) Stetbte 
where 


(a, b,c} =[a+2, b,c+1]—[a, b+2,c4+1]. (27) 


A 
—_— 
~ 


We now introduce the following substitution 
Cina f situryuide, (28) 


whence Eq. (11) may be written as 
Coove = 2A *Coo0 + 846A BCioo 
+ B*(C200+ Co20) — 46.4 F(C200+ Co20 — Cooz) 
— §BF(Csoot+2Ci2o—Cio2). (29) 


Applying the Ritz method we must minimize 
the « with respect to A and B (i.e., de€/0A =0, 
de/8B=0). The solution of the resulting two 
equations gives the values of A and B in terms 
of F—and so gives the perturbed wave function. 

If a represents the polarizability of the system, 
we have 


«= — }aF’. (30) 


To convert a@ to c.g.s. units from the atomic units 
used it is necessary to multiply by 0.7410 
X 10-75( = h*/2mie*). 

Solutions were made with both third-order and 
sixth-order wave functions for He and H-.* The 
following polarizabilities and constants found for 
the perturbed wave functions are given in 
Table I. 

I wish to thank Dr. S. Chandrasekhar for sug- 
gesting this problem. 


3 The constants for the sixth-order wave function of H 
were kindly supplied by Mr. Ralph E. Williamson in 
advance of publication. 
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The Coefficients of Thermal Diffusion of Neon and Argon and 
Their Variation with Temperature 
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In the case of a mixture of two isotopic molecules, D, the coefficient of ordinary diffusion, and 
Dr, the coefficient of thermal diffusion, are related by the equation: Dr/D =aC,C2, where C; and 
C; are the relative concentrations of the molecules, and a@ is the thermal diffusion constant. 
Experimental values of a for neon and argon in seven different temperature intervals from 
90-720°K are given. For both gases it appears that a varies linearly with the logarithm of the 
absolute temperature. These results are compared with values obtained from viscosity data and 
measurements made on binary mixtures of the noble gases. There is poor agreement with values 
calculated from the Sutherland and Lennard-Jones 9,5 models. In particular, no negative values 
of a are observed in the neighborhood of the critical temperature of argon as the 9,5 model 


predicts. 


INTRODUCTION 

UPPOSE we have a system consisting of two 
bulbs joined by a connecting tube. If the 
system is filled with a mixture of two gases, and 
if one of the bulbs is maintained at a higher 
temperature than the other, then it will be found 
upon analysis that the relative concentration of 
the heavier component will, in general, be greater 
in the cold bulb than in the hot bulb. In the case 
of isotopic molecules, the change in concentration 
at equilibrium can be expressed by means of the 

following equation: 


grad C,= —(aC,C2/T) grad T, (1) 


where a is the thermal diffusion constant; C; is 
the relative concentration of the heavy com- 
ponent; C2 is the relative concentration of the 
light component; T is the absolute temperature. 

An experimental determination of the thermal 
diffusion constant of neon for three different 
temperature ranges was made by Nier.! His re- 
sults showed a definite variation of the coefficient 
with temperature. Jones? has made an extensive 
theoretical investigation of the problem and has 
proved that such an effect is to be expected for 
both the Sutherland and the Lennard-Jones 9,5 
molecular models. One of the many interesting 
implications of the latter model is the possibility 
of negative values for the coefficient at tempera- 
tures in the neighborhood of the critical tem- 
perature. 


1A. O. Nier, Phys. Rev. 57, 338 (1940). 
2 R. Clark Jones, Phys. Rev. 59, 1019 (1941). 


The following considerations initiated the work 
to be presented here. In the first place, it was 
desirable to measure a over a large number of 
different temperature intervals so as to deter- 
mine, if possible, the relationship between a 
and TJ. In the second place it was possible to 
test the predications of Jones’ theory, as regards 
the negative coefficient, in the case of argon 
whose critical temperature is 151°K. Neon which 
has a very large value of a serves as a good com- 
parison gas. 


I. APPARATUS 


The twin-bulb method described by Nier* was 
employed here in unmodified form. To establish 
the desired temperature gradients electric fur- 
naces were used in the high temperature ranges 
and constant temperature baths at the low tem- 
peratures. The temperature of the connecting 
tubes and stopcocks was kept constant by 
means of a specially constructed circulating water 
bath. 

The mass spectrometer used for analyzing the 
gas samples was similar to one described by 
Nier‘ for routine gas analysis. A 60° magnetic 
analyzer was used. Positive ions were formed by 
electron impact of the gas at low pressure and 
the ion currents were measured with an elec- 
trometer tube amplifier. The voltages for the 
various components of the apparatus were sup- 
plied by specially designed electronic regulating 


3A. O. Nier, Phys. Rev. 56, 1009 (1939). 
4A. O. Nier, Rev. Sci. Inst. 11, 212 (1940). 
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devices. No storage batteries were used. This was 
a unique feature of the mass spectrometer. 


II. RESULTS 
A. Argon 


Argon has three isotopes® of mass 36, 38, and 
40. The percent abundances are 0.307, 0.061, and 
99.632, respectively. The concentration of A** is 
so small that it can be neglected. The initial 
pressure in the system was 40 cm Hg. The hot 
and cold bulbs were allowed to come to equi- 
librium twelve times, and after each equilibrium 
condition was attained, the hot bulb was evacu- 
ated. 

To be assured of equilibrium conditions the 
relaxation time ¢,- must be known. An approxi- 
mate calculation shows that this is given by the 
following formula: 


1 
bre = , ( 


(—)(= 
= ow, 


where A is the average cross section of the con- 
necting tube; D is the coefficient of ordinary 
diffusion ; / is the length of the connecting tube; 
V is the volume of one of the bulbs; 7) is the 
temperature of the hot bulb; 7, is the tempera- 
ture of the cold bulb. Since D is given approxi- 
mately by 


bt 
— 


D=1.4n/p, (3) 


where 7 is the viscosity of the gas and op is its 
density, it is seen that ¢,, is directly proportional 
to p. Thus the relaxation time decreases with each 
run. In all cases a time five times the relaxation 
time was allowed for each run, with thirty 
minutes for the minimum time. By running two 
different samples of gas under identical condi- 
tions but with different relaxation times, an 
experimental check was obtained. 

For accurate analysis of the gas in the mass 
spectrometer the background under the A*® peak 
must be zero, or one must be able to correct 
for it. Also a sufficient number of readings of 
the relative abundance must be taken to ensure 
a small probable error. Background may be due 
to poor resolution in the instrument causing a 





5A. O. Nier, Phys. Rev. 50, 1041 (1936). 


TABLE I. Values of the coefficient of thermal diffusion of 
argon in seven different temperature ranges. 








Temp. 
range°K TK A(A%*/A%)> a Rr 


90-195 129 0.94% 0.00315 +10% 0.07+10% 

90-296 154 2.38% 0.0709 + 5% 0.154 SY 
195-296 238 1.94% 0.0116 + 5% 0.254 5% 
195-495 300 4.68% 0.0146 + 5% 0.314 SY 
273-623 400 5$.22% 0.0182 + 5% 0.394 5% 
455-685 555 3.32% 0.0218 + 5% 0.474 5% 
638-833 720 2.73%° 0.0250 + 5% 0.534 5% 





« Ten repeated runs; in all other cases 12. 
+ This number is the average of a total of sixteen readings on two 
different samples. In some cases three different samples were used 


tailing out of the high intensity peak, as well as 
to residual impurities in the instrument itself or 
in the sample to be analyzed. A** and A“, whose 
relative mass difference is ten percent, were 
completely resolved by the mass spectrometer. 

Equation (1) may be rewritten in the form: 

(T1/To)+1 
a=[(C2°— C2) nC2C, | en - . (4) 

In (71/T 0) 
where 7; is the temperature of the hot bulb; 7» 
is the temperature of the cold bulb; C2 is the 
concentration of A** in the cold bulb; C2° is the 
concentration of A** in the normal argon; C, is 
the concentration of A*® in normal argon; m is 
the number of runs. Since (C2°—C2)/C2 repre- 
sents the fractional change in the relative abun- 
dance of the argon isotopes in the concentrated 
argon sample relative to normal argon this ratio 
can be obtained from measurements of the rela- 
tive abundances A**/A* in both the concentrated 
and normal argon samples. The procedure em- 
ployed was first to make ten determinations of 
the ratio for the concentrated sample, then to 
introduce the normal argon sample and make ten 
determinations for this; averaging these ratios, 
taking the difference and dividing by the average 
ratio of the normal argon one obtains the factor 
(C2°—C2)/C2. Ten separate values were obtained 

in this way and a grand average taken. 

This average, however, cannot be substituted 
directly in (4) since a correction has to be made 
for the gas in the connecting tube which is not 
at the temperature 7». When the volume of the 
connecting tube and the volumes of the bulbs 
were taken into account, the results given in 
Table I were obtained for argon. The values in 
column three represent the average of at least 
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two determinations. In most cases three separate 
determinations were made. 

It is obvious that @ increases with the tem- 
perature. However, since a is measured over a 
temperature interval we are faced with the 
problem of assigning a proper temperature to 
the measured value. Harrison Brown® has shown 
that the proper temperature 7, is given by 


T,=(T1—To)/T:To-In (T1/T»). (5) 


This is an adequate approximation for small tem- 
perature differences. 

When the first measurements were made it was 
observed that a varied linearly with In 7,. This 
is shown by the lower curve in Fig. 1. The em- 
pirical equation for this curve is 


Rr=0.25 In (T,/86.9). (6) 
Here Rr is the ratio of the experimentally deter- 
mined value of a to that calculated for hard 
spheres by means of Enskog’s formula: 


a=105/118(M2—M,)/(M2+M)), (7) 


where the M’s are the relative masses of the two 
molecular species. 

The results presented in Table I can be com- 
pared with values obtained for binary mixtures 
of the noble gases and those deduced from vis- 
cosity data. Table V gives the results obtained 
by Atkins, Bastick, and Ibbs’ from measurements 
made in binary mixtures of the noble gases. In 
all cases the temperature range was 273—373°K. 
This corresponds to 7,=314°K. The value of Rr 
for A-Ne is 0.54; that for A-Kr is 0.19. Since the 
hardness of argon is between that of neon and 
krypton it would be expected that the value of 
A-A would be between 0.19 and 0.54. From Fig. 1 


TABLE II. Comparison of various values of Rr for argon. 





R 





Rr (c) T 

; Rr R7(n) (Suth.) (L-J 9,5) 

T, (Meas.) (Viscosity) C =142° e/k =124° 
129 0.07 —0.44 
154 0.15 0.04 —0.31 
238 0.25 0.18 0.27 0.01 
300 0.31 0.34 0.38 0.17 
400 0.39 0.41 0.51 0.33 
555 0.47 0.47 0.60 0.47 
0.53 0.53 0.68 0.53 
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6 Harrison Brown, Phys. Rev. 58, 661 (1940). 
7B. E. Atkins, R. E. Bastick, and T. L. Ibbs, Proc. Roy. 
Soc. A172, 142 (1939). 
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Fic. 1. The dependence of Rr for neon and argon on 
the temperature. 


it is seen that the value for Rr is 0.325, which is 
in good agreement. 

If the viscosity of a gas varies as 7", it can be 
shown? that 

Rr=1.7(1—n). (8) 
Trautz® has given the experimentally determined 
values of » for argon as a function of tempera- 
ture. From these and Eq. (8), values of Rr(n) 
were calculated and compared in Table II with 
those measured by thermal diffusion. The agree- 
ment is much better than was expected. 

We are now in a position to compare the results 
with the predication of the Sutherland and 
Lennard-Jones 9,5 models. According to Jones? 
Rr can be written in terms of the constant C 
which appears as the Sutherland constant for 
viscosity data and this results in the expression 


1—0.98(C/T) 


Rr= . (9) 
1+0.92(C/T) 


TABLE III. Values of the coefficient of thermal diffusion of 
neon in seven different temperature ranges. 








Temp. range°K T°K A(Ne?2/Ne2°)¢ at5% Rr+5% 
90-195 129 3.92%> 0.0162 0.39 
90-296 154 6.08%° 0.0187 0.44 

195-296 238 2.62%" 0.0233 0.55 
195-490 298 5.47%" 0.0254 0.60 
302-645 432 5.33%° 0.0302 0.71 
460-638 545 2.10%° 0.0318 0.75 
621-819 712 2.34%" 0.0346 0.82 


« Five repeated runs. 
» Four repeated runs. 
¢ The number given is the average of a tota! of ten determinations on 


two different samples. 
8M. Trautz and H. Binkele, Ann. d. 
(1930). 
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Viscosity data may be fitted rather well by a 
value of C=142°K. With this value of C the 
values of Rr(c) are given in Table II. The agree- 
ment with experiment is poor. 

In the case of the 9,5 model the agreement is 
no better than was obtained by the simpler 
Sutherland model.* However, in this case a more 
reasonable value of (e/k), (124°), was used. 
In the case of the Sutherland model (€/k) = 823°K 
which is about five times the critical temperature. 


B. Neon 


> 


Neon has three isotopes Ne*®, Ne®', Ne* with 
percent abundances’® of 90.00, 0.27, and 9.93, 
respectively. What was said concerning argon 
applies here equally well except for the following 
considerations: (1) Since the value of Rr for 
neon is much greater than for argon the separa- 
tion per run is larger. Only four or five runs are 
required in order to obtain an appreciable separa- 
tion. (2) For this reason the initial pressure of 
the gas can be much lower (22 cm Hg). Also 
neon has a smaller density and greater viscosity 
than argon. Therefore, the relaxation time is 
much smaller. (3) No background difficulties in 


TABLE IV. Comparison of various values of Rr for neon. 








Rric) Rr(9,5)¢ 
Rr Rr(n) (Suth.) (L-J 9,5) 
Tr (Meas.) (Viscosity) > =60° e/k =42.5° 
129 0.39 0.38 0.36 
154 0.44 0.46 0.45 
238 0.55 0.61 0.62 
244* 0.56* 0.59 
298 0.60 0.76 0.66 
333* 0.63* 0.62 
423* 0.69* 0.64 
432 0.71 0.81 0.71 
498* 0.73* 0.64 
545 0.75 


712 0.82 0.85 0.73 








* These values are obtained from the graph in Fig. 1. 


TABLE V. Values of Rr for mixtures of the noble gases over 
the temperature range 0-100°C. 











Ne A Kr Xe 
He 0.80 0.65 0.63 0.59 
Ne 0.54 0.51 0.43 
A 0.19 0.17 
Kr 0.08 











* A complete description of the calculations is too lengthy 
to be presented here. The reader is advised to refer to 
Jones’ paper, reference 2. 

1 Vaughan, Williams, and Tate, Phys. Rev. 46, 327 
(1934). 
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analysis were encountered with neon. The isotope 
ratio measured is ~ 1/10 as compared with 1/300 
for argon. (4) Since the percent separation was 
larger than in the case of argon the average of 
only five differences was taken instead of ten. 
The results for neon are given in Table IIT and 
in Fig. 1. The empirical relation is 


Rr=0.25 In (T,/26.6). (10) 


These results check excellently with those ob- 
tained by Nier' for 129, 154, and 426°K. His 
values of Rr for these temperatures are 0.39, 
0.44, and 0.71, respectively. 

From the results obtained by Atkins, Bastick, 
and Ibbs we see that Rr for Ne-Ne should have 
a value between 0.80 and 0.54. This corresponds 
to a value of T, equal to 314K. From Fig. 1 we 
obtain the value 0.62 which is consistent with 
their results. 

In Table IV we have compared the experi- 
mental values of Rr for neon with those calcu- 
lated from viscosity measurements, the Suther- 
land model with C=60°, and the 9,5 Lennard- 
Jones model. 

From the empirical equation for argon it is 
seen that Rr is zero at 86.9°K. By the Lennard- 
Jones model this should occur at 220°K. Hence, 
negative values of Rr should be obtained for 
129°K, 154°K. Experimentally the values of Rr 
at these temperatures are 0.07 and 0.151. 

In Fig. 1 the effect of the greater hardness of 
neon is clearly demonstrated. The significance of 
the fact that both curves have the same slope is 
not immediately evident. However, since neon 
and argon are both spherically symmetrical and 
monatomic it might be expected that curves for 
all gases of this type will have the same slope, 
the difference in hardness being taken care of by 
the ‘‘characteristic temperature.” 
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Cloud-Chamber Track of a Decaying Mesotron 


RALPH P. SHutTtT, S. DEBENEDETTI, AND T. H. JOHNSON 


Bartol Research Foundation of the Franklin Institute, 
Swarthmore, Pennsylvania 


October 7, 1942 


N a previous letter! we described a track of what was 

evidently a negative mesotron at the end of its range 
in the gas of the cloud chamber. As with most of the other 
five or six published photographs showing mesotron stop- 
page in the gas there was no evidence in our picture that 
the mesotron had disintegrated and until now a single 
photograph by E. J. Williams? has constituted the only 
indisputable cloud chamber evidence for mesotron decay. 
It is therefore of interest to put on record another track 
of a mesotron disintegration recently obtained during first 
tests of a high pressure cloud chamber. 

The details of the chamber will be described in another 
article and it will suffice to state that the chamber is 
cylindrical, 30 cm in diameter and 10 cm deep. It is de- 
signed to operate at an internal pressure of 200 atmos- 
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Fic. 1. Photograph of a mesotron disintegrating at C. 
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Fic. 2. Decaying mesotron of Fig. 1, enlarged. 


pheres but the present photograph was obtained at a 
pressure of 70 atmospheres in a mixture of argon, n-propyl 
alcohol, and water vapor, without magnetic field or counter 
control. The expansion ratio decreases as the pressure is 
increased and at 70 atmospheres good tracks were obtained 
at a volume ratio of 1.05. 

Figure 1 shows the entire field of the chamber, and Fig. 2 
is an enlarged view of the track of particular interest. 
What is evidently a mesotron enters the chamber at A, 
at B the track is definitely heavy, and at C a conversion 
takes place in which the mass of the mesotron disappears 
and a lightly ionizing particle of high energy, presumably 
an electron, flies off at an angle of 85° with the direction of 
the primary mesotron track. 

It is interesting to observe the apparent change in 
density of the mesotron track 5 cm from its end at B. 
According to the Bethe-Bloch formula the residual range 
of a mesotron of mass 180 m is 6 cm after the ion density 
in the track has become four times that of a ray of mini- 
mum ionization. Thus the apparent change in density 
occurring over a relatively short distance at this pressure 
is probably not an effect of uneven illumination. The 
curvature of the mesotron track is presumably caused by 
multiple scattering, and may be used for an estimate of 
the mass of the mesotron. According to Williams*® p,/R 
=1.3(m/Z)', where ps is the average radius of curvature 
of the track over the first half of its range R, m is the mass 
number of the particle, and Z the atomic number of the 
scattering gas. Examination of the track shows that pz, is 
between 2 and 6 times the residual range which places the 
mass of the particle between 40 and 400 electron masses. 
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The lightly ionized track of the disintegration particle 
shows no deflection, and taking into account the high 
pressure of the gas we may conclude that its energy exceeds 
3X10? ev whereas the density of the mesotron track indi- 
cates that its final energy was less than 106 ev. 

In conclusion the picture seems to represent a slow 
particle of intermediate mass converted into a fast particle 
of small mass, the type of disintegration which has been 
postulated to account for certain anomalies in the absorp- 
tion of cosmic rays in the atmosphere. Since this is the 
first track found in our experiments out of a total 2 10®cm 
of mesotron tracks examined we have an estimate that the 
probability of disintegration is 5X10-7 per cm. In our 
high pressure chamber we can have an equivalent path 
length of 6000 cm for each ray and with the long time of 
sensitivity realized we may have several rays in each 
photograph. Thus we may expect to find events of this 
type occurring once in something of the order of every 
one hundred photographs at maximum pressure. 

Acknowledgment is made of the generous support given 
for the high pressure cloud-chamber project by the Car- 
negie Institution of Washington. 

1T. H. Johnson and R. P. Shutt, Phys. Rev. 61, 380 (1942) 


2? E. J. Williams and G. E. Roberts, Nature 145, 102 (1940). 
. J. Williams, Phys. Rev. 58, 292 (1940). 
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The Beta-Ray Spectrum of Antimony (124) 


E. B. HALES AND E. B. JORDAN 
University of Illinois, Urbana, Illinois 
November 12, 1942 


EUTERON bombardment of antimony gives rise to 

a large number of radioactive isotopes. Sb'™ and Sb"* 
are formed from the two stable isotopes Sb™ (56 percent) 
and Sb"3 (44 percent) by a (d, p) reaction. The periods 
associated with these two radioactive isotopes have been 
identified by producing them in a number of different 
ways.! The average of the values reported for the half- 
life of Sb! is 61 hours, and that of Sb!‘ is 60 days. Both 
isotopes emit hard gamma-rays. 

In addition to the radio-antimony isotopes, a number of 
tellurium isotopes are also formed by deuteron bombard- 
ment; Te!! is formed by a (d, 2m) reaction, and has a half- 
life of 125 days.? There is some evidence’ that it is a posi- 
tron emitter giving rise to several gamma-rays.** Other 
tellurium isotopes which are probably formed at the same 
time are those having mass numbers 122, 124, and 123; 
the first two by a (d, ”) reaction and the latter by a (d, 2) 
process. Although these isotopes occur as stable in nature, 
there is some evidence* that two of them, Te! are 
formed in excited states by the deuteron reaction, yielding 
gamma-rays having energies of 82.0 kev and 88.3 kev. 
Other gamma-rays which are emitted by the separated 
tellurium fraction and whose origins are not, definitely 
known, have energies of 136,157.3, 210.8, and 615 kev.’ 
A period of 30 days has been assigned to the composite 
gamma-ray activity. 

The present investigation was undertaken in coopera- 
tion with Professor Goldhaber. A source produced by 
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deuteron bombardment of antimony was prepared for him 
by Dr. Kamen of Berkeley. A fraction of this source, in 
the form of a powder scraped from the surface, was at- 
tached to a copper foil, 0.001” thick, by pouring on a very 
small ‘amount of Vinylite dissolved in methyl ethyl ketone. 
The surface density of the source, including all materials, 
was estimated to be 30 mg/cm? 

More than seventy days after the sample was bom- 
barded, measurements were started on the continuous 
(negative electron) spectrum of Sb"*. The beta-ray activity 
of Sb’ was, by this time, reduced to approximately one 
part in two hundred million. The only other long-lived 
activities present were those due to a 30-day gamma-ray® 
and to the 125-day Te™ activity. Since the latter decays 
either by K-electron capture or by positron emission, there 
were no beta-rays due to these activities except possibly 
those due to internally converted gamma-rays. In addition, 
measurements made at the momentum intervals specified 
by 3170 gauss-cm (0.57 Mev) and 5680 gauss-cm (1.26 
Mev) gave half-lives of 59 days and 624 days, respectively. 

All relative numbers of counts were reduced to a stand- 
ard time on the assumption that the 60-day half-life given 
by Livingood and Seaborg is correct. 

Figure 1 is a plot of the relative number of counts per 
momentum interval against the corresponding momenta 
specified in gauss-cm. The probable errors in the relative 
numbers are given approximately by the lines through the 
points. A greatly magnified drawing of the end-point region 
is shown in the upper right-hand corner. Figure 1 clearly 
indicates that the continuous spectrum consists of two 
components. 

The end points were determined both by inspection and 
by the Fermi plots. The two sets of values are in good 
agreement. Allowance was made for the fact that the last 
beta-rays to enter the counter were grazing the inner slit 
jaw. The upper end point did not include a small tail on 
the observed curve. The end points assigned were 3800 
+100 gauss-cm and 9720+200 gauss-cm. These probable 
errors are due largely to the magnetic field calibration. 
They would be about one-half as large on the basis of 
internal consistency alone. These end points correspond 
to energies of 0.74+0.03 Mev and 2.45+0.07 Mev. 
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Fic. 1. The beta-ray spectrum of Sb". The curve is plotted 
with arbitrary ordinates. 
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The difference between these end points, 1.71+0.07 Mev 
agrees with the known value of the energy of the hard 
gamma-ray of Sb!*, Using the absorption method, Mitchell, 
Langer, and McDaniel‘ obtained a value of 1.82+0.05 Mev 
for this gamma-ray. More recent measurements by Scharff- 
Goldhaber and Klaiber® gave a value of 1.75+0.02 Mev. 
The internal conversion lines indicated at points A and 
B (~0.23 Mev, 0.61 Mev) correspond to two of the in- 
ternally converted gamma-rays reported by Kent and 
Cork’ as emanating from the separated tellurium fraction. 
1 J. J. Livingood and G. T. Seaborg, Phys. Rev. 55, 414 (1939). 


2G. T. Seaborg, J. J. Livingood, and J. W. Kennedy, Phys. Rev. 57, 
363 (1940). 
3C. V. Kent and J. M. Cork, Phys. Rev. 62, 297(A) (1942). 
4A. C. G. Mitchell, L. M. Langer, and P. W. McDaniel, Phys. Rev. 
57, 1107 (1940). 
5G. S. Klaiber and G. Scharff-Goldhaber, Bull. Am. 
7 (1942). 
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Spin Inversion Processes and 
Nuclear Spectroscopy 


Guipo BECK AND J. L. RODRIGUES MARTINS 
Laboratorio de Fisica, da Universidade de Coimbra, Portugal 
November 11, 1942 


T has been shown recently! that the resonance scattering 
due to compound states of Li; has to be described by a 
superposition of s; and p; waves. We are thus led to con- 
clude that spin inversion processes can play a considerable 
role in nuclear collisions and that orbital momentum 
cannot be considered to be an integral of motion. 

In the case of a two-body system, the simultaneous 
occurrence of partial states belonging to different values 
of angular momentum has to be attributed to interaction 
terms of the type 

g(r) (or) Cer) /r?, (1) 


which may be either ordinary spin-spin forces or exchange 
forces. They imply that the ground state of the deuteron 
becomes a mixture of a 4S and a *D state. The numerical 
value of the quadrupole moment, measured by Kellogg, 
Rabi, Ramsey, and Zacharias,’ indicates that the ratio of 
the *S and the *D amplitude is of the order 50 : 1. We can 
conclude from this fact that the interaction term (1) 
amounts only about to one-tenth of the forces responsible 
for nuclear binding. 

We have examined the question whether small forces of 
the type (1) can account for spin inversion processes 
occurring with considerable probability. In the case of a 
two-body problem rigorous solutions can be derived if we 
put in (1) 

g(r) = h®«/8x?mr’, (2) 


where ¢ is a numerical parameter measuring the intensity 
of the spin-spin coupling. Within the energy range in 
which we can expect marked resonance phenomena, 
0<2x2R/X<1, the behavior of the resonance scattering 
depends essentially on the ratio of (2rR/A)4‘ and ¢, where 
Ai is the change of spin momentum during an inversion 
process. 

In the case of extremely large wave-lengths, (27R/d)4' Ke, 
the spin-spin coupling has to be considered to be a strong 
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one. Spin inversion leads like other inelastic collision 
processes, to a broadening of resonance levels. Even if 
inside the compound system the spin inversion probability 
is very considerable, only s-neutrons enter and leave the 
compound system with appreciable probability. Resonance 
scattering of slow neutrons will, therefore, always show 
spherical symmetry. In the opposite case, (2rR/A)4'>«, 
the spin-spin coupling is weak. Nuclear resonance levels 
behave according to the value of orbital momentum to 
which they mainly belong. If, however, (27R/A)4i~e, 
narrow resonance levels can exist, for which scattering 
waves of different orbital momentum are simultaneously 
excited and which give rise to a neutron scattering depend- 
ing on the scattering angle. Using our simplified model, 
which does not take into account features of the many- 
body problem, and introducing Rabi’s value of the quad- 
rupole moment, we find that this effect can already occur 
for energies of the order of 10° ev. Though this numerical 
value may not be very reliable for compound states of a 
many-body system, we may still conclude that spin in- 
version processes may occur, as in the case of Lis, with 
considerable probability even if they are due to small 
forces. 

A detailed account on this investigation will be given 
by one of us (R. M.) in a paper to appear in the Revista da 
Faculdade de Ciéncias de Coimbra. 

Cahiers 


'G. Beck and Tsien San-Tsiang, Phys. Rev. 61, 379 (1942); 


de Physique (June, 1942). 
2 J. M. B. Kellogg, I. I. Rabi, N. F. 
Phys. Rev. 57, 677 (1940). 
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Note and Erratum: Evidence of a Periodic 
Deviation from the Schottky Line. II.! 
T. E. Pxuipps 
Noyes Chemical Laboratory, University of Illinois, Urbana, Illinois 
November 8, 1942 
I N a Letter to the Editor Nottingham? referred to an 
apparent disagreement between the results of Seifert 
and Phipps? and those of Turnbull and Phipps! regarding 
the algebraic sign of the temperature dependence of ‘‘the 
amplitude of the deviations’? from the Schottky line, and 
stated that his own observations indicated that this ampli- 
tude is independent of the temperature. It is probable that 
an error in an equation in reference 1, which will be cor- 
rected at the end of this letter, is in large part responsible 
for Nottingham’s criticism. Beyond the correction of this 
error, however, it seems to the writer that a sharp re- 
definition of terms and a survey of the experimental results 
are needed. 
In applying the Schottky equation, 


logio i =logio io+ (4 E}/2.303kT), (1) 


to experimental data, the intercept, logio 79, may be deter- 
mined by the method of least squares (as in reference 3) or 
chosen arbitrarily (as in reference 1). For a particular 
experimental point a quantity A logio i may be defined by 
the equation, 


A logio i=logio i—logio lo— (e§E4/2.303kT). (2) 





an 
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It is apparent that A logioz is the deviation of the experi- 
mental point from a line having the Schottky slope and 
an intercept, logio 70. Since A logio i has been found experi- 
mentally* to be a periodic function of FE}, it becomes con- 
venient to assign the symbols A logio zm and A logio i, to 
adjacent extreme values (whether maximal or minimal) of 


, 
the deviation, and to define a quantity 4(m,,) by the 
relation, 


5¢m,n) = | (A logio im — A logio tn) |. (3) 


aw (abbreviated hereafter to 6’) may be called “the 
oscillation of the A logio 7, E* function in the E? interval, 
(m,n).” It is easily seen that 65’ is independent of the 
arbitrarily chosen intercept, logio io. The statement made 
in reference 3, pp. 657-658, namely, “‘the heights of the 
maxima and minima . . . decrease . . . with increasing 
temperature” could be restated thus: ‘“dé’/dT has a 
negative sign.” ; 

The quantity 6 employed in reference 1 is defined by 
the equation, 


a 100) [Foes aaa |. (4) 
logic tm logio tn }| 

It can readily be shown that 8 m,,) (abbreviated hereafter 
to 5), though insensitive to the choice of the intercept, 
logio Zo, is not entirely independent of it. For this reason it 
is regrettable that the quantity 6 should have been intro- 
duced; and it is suggested that 6’ be employed in future 
discussions. 

With regard to the temperature dependence of 6, the 
data of references 1 and 3 agree that dé/dT is positive. 
With regard to the temperature dependence of 6’, the 
data of references 1 and 3 indicate, but do not establish 
conclusively, that d5'/dT is negative. In Table I appears a 
survey of dé’/dT values, calculated from representative 
data of reference 3 (namely, that of Fig. 3), and of refer- 
ence 1 (namely, that of Figs. 2 and 3). 

The following comments may be made upon Table I. 


Tasxe I. Survey of values of dé’/dT. 








(d8'/dT) X 108 
Ref. 3, Fig.3 Ref.3, Fig.3 ‘Ref. 1, Fig. 2 Ref. 1, Fig. 3 
Ent Ent L.S.2 Schottky? Schottky? Schottky? 
200 266 —24(41.0) —03(40.7) — 2.0 (45.8) 
266 348 —4.9 (+1.2) —7.7 (+1.8) —15. (+3.0) —10. (+2.0) 
348 «505 «= —7.6(41.4) —1.7(+3.8)" — 28(440)* — 5.4 (41.5) 





505 752 nelle —38. (+12.) _ 


@ LS. indicates that slope and intercept were calculated by least squares. 

+ Schottky indicates that the Schottky slope was employed, with an arbitrary 
intercept. 

* Data badly scattered. 


1. In the £! interval (200, 266), large experimental error 
renders the sign and magnitude of d5’/dT uncertain. 2. In 
the intervals, (266, 348) and (348, 505), (except for the 
two cases starred in Table I, in which the data were badly 
scattered) the sign of dé’/dT is consistently negative. An 
average value of dé’/dT is —7X10-® deg.—'. 3. In the 
interval (505, 752) the single very large negative value for 
dé’/dT may point to an important effect of field upon 
this coefficient. 
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In conclusion, it would appear that more accurate data, 
especially at high fields where 3’ is large, will be necessary 
to decide whether dé’/d7T is zero* as the observations of 
Nottingham? show; or whether it is negative (and possibly 
field-dependent) as the observations of references 1 and 3 
indicate. 

The second equation on page 667 of reference 1 should 
be corrected to read 


dé= eo [+0 logio D,+(MqdV)/(2En') ] 
10g 10 @m 
100 . — 
+——-[+0 logio Dz +(MqdV)/(2E,') }. 
logio tn 


1D. Turnbull and T. E. Phipps, Phys. Rev. 56, 663 (1939). 

2W. B. Nottingham, Phys. Rev. 57, 935 (1940). 

3R. L. E. Seifert and T. E. Phipps, Phys. Rev. 56, 652 (1939). 

‘It is assumed that Nottingham’'s “amplitude of the deviations” 
corresponds to @’. 





Gamma-Rays from Sc* 


C, E. MANDEVILLE 
The Rice Institute, Houston, Texas 
November 24, 1942 


HE radioactive isotopes formed by the reaction 

Ca+deuterons have been investigated extensively, 
and the particles emitted by them have been studied in 
detail.!~* Walke! has ascribed with certainty a half-period 
of 44+1 hr. to Sc*’, and his absorption experiments have 
indicated that Sc** emits two 8-ray spectra, one intense 
and of maximum energy 0.50 Mev, the other weak and of 
maximum energy 1.4+0.1 Mev. He reported also a y-ray 
of energy 0.9+0.1 Mev which appeared to be related to 
the two 8-ray spectra. More recently, G. P. Smith‘ has 
examined the 8-ray spectrum of Sc** by means of a spec- 
trometer and has found one continuous distribution of 
maximum energy 0.640+0.004 Mev. No conversion 
electrons arising from the y-ray reported by Walke! were 
observed nor was there any evidence of the weak disinte- 
gration electron distribution of maximum energy 1.4+0.1 
Mev. 

Sc** has been produced by Ca+deuterons and a further 
study made of its gamma-radiation. Compton recoils of 
the y-rays from Sc*® have been observed with a y-ray 
spectrograph which has been previously described.’ The 
momentum distribution of the recoil electrons is given in 
Fig. 1. By comparison with earlier results,® it is seen that 
the y-rays are monochromatic. Curve A was obtained on 
receipt of the radioactive sample and about 3 days after 
bombardment. Curve B was obtained 103.2 hr. after 
curve A, and curve C 175.5 hr. after curve A. The peaks 
of curves A, B, and C are in the ratio 134/27.0/8.5 to one 
another, and the half-life calculated from the ratios is 
44+1.5 hr. It is to be noted that the shape of the curve 
remains unchanged with time. This suggests that only 
one half-period is present. The y-rays from Ca isotopes 
of long half-lives formed by Ca+deuterons were not 
present in sufficient intensity to be measured by the 
spectrograph. The counting rate of the G-M counters of 
the spectrograph was checked throughout the time of the 
experiment with a radium source and was found to vary 
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from Sc*8, 


by an amount less than the statistical probable error of 
the points on the curve of Fig. 1 

The end point of the Compton electron spectrum of 
Fig. 1 corresponds to a y-ray energy of 1.35+.03 Mev 
It is probable, therefore, that one mode of disintegration 
of Sc*® would consist in the emission of a negative electron 
spectrum of maximum energy 0.640 Mev followed by a 
y-ray of quantum energy 1.35+.03 Mev, a more precise 
value than that of Walke’s absorption method. In addition 
to information which this result gives with regard to the 
disintegration energy of Sc**, an excitation level in the 
Ti*® residual nucleus is established at 1.35+.03 Mev 
Pollard® has reported a level at 1.1 Mev in the Ti** nucleus. 
The absence of the electrons of a converted y-ray in the 
case of Smith’s experiment might be attributed to a low 
internal conversion coefficient arising from the energy of 
the y-ray and a small spin difference between the level at 
1.35 Mev and the ground state of Ti**. The presence of 
this y-ray indicates that the existence of the weak §-ray 
spectrum of high maximum energy reported by Walke! 
is possible and of maximum energy about 1.99 Mev. 

It is a pleasure to thank Professor A. L. Hughes and 
the cyclotron group of Washington University, St. Louis, 
Missouri, for the preparation of the radioactive material. 
Their effective cooperation made possible this work. 
Experiments of a similar nature are in progress. 

1H. Walke, Phys. Rev. 57, 163 (1940). 

2H. Walke, F. C. Thompson, and J. Holt, Phys. Rev. 57, 177 (1940). 

3D. R. Elliott and L. D. P. King, Phys. Rev. 60, 489 (1941). 

4 Gail P. Smith, Phys. Rev. 61, 578 (1942). 


5C, E. Mandeville, Phys. Rev. 62, 309 (1942). 
6 E. Pollard, Phys. Rev. 54, 411 (1938). 





The Diffusion Length of C Neutrons in Water 


G. R. GAMERTSFELDER AND M. GOLDHABER 
Department of Physics, University of Illinois, Urbana, Illinois 
November 25, 1942 


HE diffusion length / of C neutrons in water, which 
is the average distance C neutrons will diffuse in 
water before being absorbed, has been obtained previously 


THE EDITOR 


in two different ways. It has been computed either from 
the data given by Amaldi and Fermi! for paraffin or with 


the help of the equations 
l=d(n/3)! 

where A is the mean free path for scattering of C neutrons, 

n the average number of collisions a C neutron makes 


before being absorbed, and + the mean lifetime of slow 
The average 


and m\=@dr, 


neutrons (all quantities measured for water). 
velocity # of thermal neutrons, is assumed to be 2.2 105 
cm/sec. 

Recent computations? have given /=2.5 cm, as deduced 
essentially from the early work of Amaldi and Fermi, and 
1=2.22 cm as deduced from the work of several authors 
on \ and 7.8 

It seemed to us desirable to determine / in a somewhat 
more direct manner. For a point source of C neutrons, 
or at some distance from a spherical source, the density 
N of C neutrons in water at a distance R from the center 
of the source, should be given by the diffusion equation 


NR=const. Xe~/4, (1) 


if the scattering is assumed to be isotropic.‘ Using this 
equation, we have measured /, with the help of a “virtual’”’ 
source of C neutrons in the following way. 

Inside a large water tank, at a distance of 6.35 cm from 
a 90 mg Ra-a-Be source, a graphite sphere of 1”’ diameter 
was suspended. The relative slow neutron density Ni was 
then measured as a function of the distance R from the 


center of the sphere on the far side from the Ra-a-Be 
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The root-mean-square error of RN increases from about 5 percent for 
small values of R to about 30 percent for the larger values of R 


neutrons we used 
” high, 


source. As a detector for the slow 
small ionization chamber }” in diameter, and 13 
lined with boron carbide and connected to a linear ampli- 
fier. A similar series of measurements was carried out with 
Cd (0.5 mm thick) covering the sphere, the slow-neutron 
density in this case being N2 N=Ni—N: should then 
vary according to the diffusion Eq. (1) for C neutrons. 

In Fig. 1 NR has been plotted on a semilog scale as a 
function of R. From the slope of the straight line fitted 
to the experimental points we obtain a value /=3.1 cm. 
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A correction, taking account of the fact that our source 
is not quite spherically symmetrical and that our detector 
has a finite size, should reduce this value slightly. A round 
value /=3 cm with a probable error of 10 percent may be 
concluded from our results. This is somewhat larger, but 
barely outside the respective limits of error, than the 
value computed from the work of Amaldi and Fermi for 


paraffin. 


1 
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According to a calculation kindly carried out by Dr. P. 
Morrison, the diffusion Eq. (1) will not be appreciably 
affected by assuming a slight deviation from isotropic 
scattering for C neutrons in water. 

1 E. Amaldi and E. Fermi, Phys. Rev. 50, 899 (1936). 

2H. Feeny, C. Lapointe, and F. Rasetti, Phys. Rev. 61, 469 (1942). 

3See J. H. Manley, L. J. Haworth, and E. A. Luebke, Phys. Rev 
61, 152 (1942) for the most recent work on +r and for references to 


earlier work on A and r. 
4E. Fermi, Ric. Sci. 7, 13 (1936). 
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MINUTES OF THE MEETING OF THE NEW ENGLAND SECTION HELD AT 
HARTFORD, CONNECTICUT, OCTOBER 24, 1942 


HE twentieth regular meeting of the New 

England Section of the American Physical 
Society was held at Trinity College, Hartford, 
Connecticut, on October 24, 1942. At least 
forty-five members were in attendance. The in- 
vited papers were as follows: 


H. G. HouGutTon (Department of Meteorology, Massa- 
chusetts Institute of Technology): The Training of 
Meteorologists for the Armed Services. 


A. P. R. WapDLuND (Trinity College): Physics at Trinity 
College. 


Nora M. MouHLer (Smith College): Intra- and Extra- 
Curricular War Courses at Smith College. 


HERMAN FESHBACH (Massachusetts Institute of Technology): 
Applications of High Voltage Electrons in Nuclear 
Research. 


Abstracts of six of the eight contributed papers 
are given below. The others will appear in the 
American Journal of Physics. A vote of thanks 
was extended to the members of the Physics 
Department of Trinity College for their generous 
hospitality. 
At the business meeting the following officers 
were elected for the current year: 
Chairman: GLADys A. ANSLOW (Smith College). 
Vice Chairman: A. P. R. WADLUND (Trinity College). 
Secretary-Treasurer: MILDRED ALLEN (Mount Holyoke 
College). 
Members of Program Committee: 
C. E. BENNETT (University of Maine), 
W. W. StTIFLER (A mherst College). 
MILDRED ALLEN 
Secretary-Treasurer 


ABSTRACTS OF CONTRIBUTED PAPERS 


1. Anisotropic Solutions of Colloids. Lars ONSAGER, 
Yale University.—The solutions of certain colloids com- 
prised of highly asymmetrical particles—plates or rods— 
are known to form anisotropic phases at remarkably low 
concentrations. For tobacco mosaic virus (rods), isotropic 
solutions containing 2—3 percent virus are in equilibrium 
with anisotropic phases containing 3-4.5 percent, respec- 
tively, according to the amount of electrolyte present. 
This phenomenon can be explained as a result of repulsive 
forces by the observation that the mutual co-volume of 
two swarms of parallel rods (or plates) is roughly propor- 
tional to the sine of the angle between their orientation, 
and larger than the volume of the particles by a factor 
which is proportional to the asymmetry. The case of rods 
is particularly simple in that the virial coefficients of order 
higher than 2 in Mayer’s expansion are small, and a quanti- 
tative theory is possible. The computed ratio of concen- 
trations at equilibrium is 1.34. The predicted osmotic 
pressure of the anisotropic phase is nearly proportional 
to the concentration, in fact, slightly greater than 3cR7/V. 


2. An X-Ray Diffraction Study of Liquid Carbon Tetra- 
chloride. A. EISENSTEIN, University of Missouri.—X-ray 
diffraction patterns of liquid CCl, at 27°C have been ob- 
tained using both the photographic and Geiger-Mueller 
counter methods. These patterns have been investigated 
to a value of 1.1 sin 6/A and show seven diffraction peaks 
in this region. The corrected patterns were subjected to a 
Fourier analysis to obtain the electron density distribution 
function for the liquid. This function may be interpreted 


as indicating that each C atom is surrounded by four Cl 
neighbors at 1.85A, and that each Cl atom in turn has 
three Cl neighbors at 2.95A; thus confirming the tetra- 
hedral molecular structure that has been predicted from 
electron diffraction studies. In addition, each Cl atom has 
2-3 Cl neighbors at 3.9A, these neighbors belonging pre- 
sumably to neighboring molecules. A further electron con- 
centration at 6.2A seems to indicate that at that distance 
each CCl, molecule finds 4 nearest molecular neighbors. 
By an approximate method of calculation, this distance is 
found to increase to 7.3A when the liquid is heated to 
240°C under a pressure of 35 atmospheres. Using the equa- 
tion of Hildebrand and Wood, the energy of vaporization 
has been calculated from the distribution function and 
found to be in reasonably good agreement with the 
measured value. 


3. Thermal Diffusion Separation of Isotopes. W. W. 
Watson, R. Srmon, AND D. L. WorERNLEY, Yale Uni- 
versity.—Details of several variations of multi-stage ther- 
mal diffusion apparatus as well as improvements in the 
small Nier-type mass spectrometer are presented. Results 
and calculations for the separation of the isotopes of 
carbon, argon, and other elements have been obtained. 
An experiment is described for detecting a negative thermal 
diffusion constant, using ammonia gas at low temperature. 


4. The Disintegration of Co®. Martin DeuTscH AND 
Lioyp G. Ettiott, Massachusetts Institute of Technology.— 
The radiations from the disintegration of the two radio- 
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active isomers of Co® were studied by beta-ray spec- 
trometer and coincidence techniques previously described.' 
The 5-year isomer was produced with sufficient specific 
activity by bombarding cobalt metal with deuterons. It 
decays with the emission of a simple negatron spectrum 
of maximum energy 0.300+0.006 Mev followed by the 
successive emission of two gamma-rays of energy 1.10 
+0.02 Mev and 1.30+0.02 Mev, respectively. The 10.7 
minute isomer was obtained by exposing K;Co(CN)¢. to 
slow neutrons and precipitating the radioactive recoil 
atoms with KOH and H:,O:. We find that the direct 
transitions to Ni®, reported by Nelson, Pool, and Kur- 
batov? constitute only 10 percent or less of the disintegra- 
tions. The negatron spectrum has a maximum energy of 
1.50+0.15 Mev and is probably followed by the emission 
of a single gamma-ray. At least 90 percent of the dis- 
integrations proceed by an isomeric transition (probably 
to the 5-year level) with the emission of a gamma-ray of 
0.056+0.003 Mev or a corresponding conversion electron. 


1 Phys. Rev. 60, 470, 544 (1941); 61, 686 (1942); 62, 3 (1942). 
2 Phys. Rev. 62, 1 (1942). 


5. Further Study of Edge Tones. ARTHUR TABER JONES, 
Smith College-—Two types of edge tones are here distin- 
guished: a first type at lower wind velocities, and a second 
type at higher velocities. The first type is the one that has 


] 


wn 
vo) 


heretofore received most attention. In it gradual changes in 
wind velocity or in slit-edge distance are accompanied by 
changes in pitch that are continuous throughout each of 
several “‘stages,’’ the successive stages being separated by 
discontinuities in pitch. In the second type the changes in 
pitch are continuous throughout, and may extend over 
several octaves. The phenomena in the first type are not 
as simple as is suggested by the descriptions usually given. 
A number of possible equations connecting frequency N 
and slit-edge distance 4 are examined to see which fits best. 
For a slit width of 0.7-0.8 mm it is concluded that the most 
satisfactory of these equations is Nk*=j7U, where U stands 
for the velocity of the wind at the slit, s=1.00, 1.14, 1.22, 
1.43, respectively, and 7=3.9, 11.8, 24.0, 6.8, respectively, 
for stage I, stage II, stage III, and the second type of 





edge tones. 


6. Crystal Statistics. Lars ONSAGER, Yale University. 
The partition function for the Ising model of a two- 
dimensional ‘“‘ferromagnetic” has been evaluated in closed 
form. The results of Kramers and Wannier concerning the 
“Curie point’ 7. have been confirmed, including their 
conjecture that the maximum of the specific heat varies 
linearly with the logarithm of the size of the crystal. For 
an infinite crystal, the specific heat near T=T- is propor- 
tional to —log |(T—T-.)|. 
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MINUTES OF THE MEETING OF THE 


New YorK STATE SECTION HELD 


AT ALBANY, NEW YorK, NOVEMBER 7, 1942 


HE Fall Meeting of the New York State 
Section of the American Physical Society 
was held at the New York State College for 
Teachers and the Albany High School on Satur- 
day, November 7, 1942. An attendance of over 
100 members and visitors was recorded. 
The following program of invited papers was 
presented : 
Address of Welcome. JoHN M. SAyLes, President, New 
York State College for Teachers 
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Recent Developments in Fluorescent Materials. G. R. 
FonpDA, General Electric Company. 


Electron Diffraction Studies. R. P. JOHNSON, General 
Electric Company. 
Creep of Metals. SauL DusHMAN, General Electric Com- 
pany. 
Million-Volt Industrial X-Ray Unit. E. E. CHARLETON, 
General Electric Company. 
W. R. FREDRICKSON 
Secretary 
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MINUTES OF THE MEETING OF THE OHIO SECTION HELD 
AT CoL_umBus, On1I0, NOVEMBER 14, 1942 


HE ninth meeting of the Ohio Section of the 

American Physical Society was held in the 
Mendenhall Laboratory of Physics, The Ohio 
State University on November 14, 1942, in 
Columbus, Ohio. Twelve papers were presented, 
six being contributed papers and the latter six 
were devoted to a short Symposium on Physics 
in the Glass and Ceramic Industries. The com- 


plete programme follows. 

Action was taken that the Ohio Section co- 
operate with the American Physical Society and 
join that body in a late winter or early spring 
meeting now being planned to convene in this 
region. 

LEON E. Situ, Acting Secretary, 
Denison University, Granville, Ohio 


ABSTRACTS OF CONTRIBUTED PAPERS 


1. The Effect of Depth of Focus on the Zone of Clear 
Single Binocular Vision. HENRY WILLIAM HOFSTETTER, 
The Ohio State University.—If the range of accommodation 
is measured for various amounts of convergence and the 
range of convergence is measured for various amounts of 
accommodation, the limits of the ranges so obtained can 
be plotted on a system of coordinates. The data so plotted 
circumscribe a zone on the graph which may be called the 
zone of clear single binocular vision. The shape of this 
zone approximates a parallelogram. This zone is usually 
obtained by a subjective method which involves the de- 
tection of a blur to indicate changes in accommodation. 
Because the eye has a depth of focus within which blurring 
cannot be detected, the changes in accommodation that 
take place within the range of relative convergence are 
not ordinarily detected. By making an allowance for the 
estimated depth of focus in subjectively obtained zones, 
it is possible to explain the usually obtained deviations 
from a true parallelogram. These deviations from the true 
parallelogram have often been interpreted to have direct 
physiological correlates. This investigation indicates that 
these deviations are dependent upon the presence of a 
subjective depth of focus. 


2. The Permeability of Iron at Radiofrequencies. 
ALvA W. SMITH AND FREDERICK P. DicKEy, The Ohio 
State University.—The permeability of iron wires, of diam- 
eters 0.005-0.110 cm, was determined in the frequency 
range of 250 kc to 6 mc. A modification of the heterodyne 
method used by Wait! made it possible to make measure- 
ments on the specimens in a magnetic field of constant 
r.m.s. value and practically free from superposed d.c. 
fields. The apparent permeability of iron was found to 
decrease as the frequency increases in all of our specimens. 
Wait reported an increase in apparent permeability with 
increase in frequency. We observed no anomalous changes 
in permeability, as has been reported by some observers. 
Our results agree with Wait’s in regard to the existence 
of anomalous changes in permeability. The variations 


of permeability with wire Size and frequency can be ex- 
plained, at least qualitatively, as due to eddy current in 
the test specimens. 

1 Wait, Phys. Rev. 32, 967 (1928). 


3. A Permalloy Magnetometer. E. H. JoHNson, Kenyon 
College-—For the purpose of locating buried scrap iron, a 
simple magnetometer was developed. It consists essentially 
of a vertical permalloy rod, near the upper end of which 
a magnetic compass is mounted so as to be movable along 
a horizontal scale. The distances of the compass from the 
rod indicate the variations in the intensity of the magnetic 
field in the vicinity of the lower end of the rod. Several 
laboratory uses also are suggested. 


4. A Problem in Electrostatics. L. H. THomas, The 
Ohio State University.—The solution of the two-dimensional 
problem of the field between an infinite plane and a per- 
pendicular semi-infinite thick plate separated from it by 
an arbitrary gap is well known. If the gap is small com- 
pared to the thickness of the plate, the field in the gap is 
nearly uniform, but the field intensity becomes infinite at 
the edges of the plate. It is possible to modify the shape 
of the end of the plate so as to make the field intensity 
constant over it, in which case the intensity everywhere 
else is less. The solution of this problem is contained in 
the conformal transformation x to Z given in terms of 
elliptic functions by 


x=am(w, k) 
Z=ksn(w, k)+E(w, k), 


which can be obtained by the method of Kirchhoff for 
treating similar hydrodynamical problems involving free 


stream lines. 


5. A Comparison of the Absorption Spectra of Some 
Biological Stains with the Spectra of Sections of Various 
Tissues Stained with Them. V. M. ALBers, The C. F. 
Kettering Foundation for the Study of Chlorophyll and 
Photosynthesis, Antioch College, Yellow Springs, Ohio AND 
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RoBeERT E. STOWELL, The Barnard Free Skin and Cancer 
Hospital, St. Louis, Missouri.—The absorption coefficients 
of eosin in ethyl alcohol, aqueous azocarmine G, aqueous 
picric acid, aqueous analine blue, aqueous acid fuchsin, 
and aqueous orange G staining solutions have been meas- 
ured throughout the visible region. The absorption spectra 
of 10u and 20, sections of gelatine, tendon, blood plasma, 
_and thymus have been measured after staining with the 
various stains and their spectra compared with the spectra 
of the staining solutions. The spectrum of aqueous basic 
fuchsin and feulgen stain reduced by formaldehyde have 
also been measured and compared. It was found that the 
spectra of the stained sections are very similar to the 
spectra of the original staining solutions, indicating that 
the staining process is a physical one not involving any 
chemical reactions with the tissues. In the case of the 
feulgen reaction it has generally been assumed that the 
product obtained on reduction with formaldehyde is the 
original basic fuchsin. The spectrum measurements indi- 
cate that this reaction product has a different absorption 
spectrum from the basic fuchsin solution with absorption 
coefficients more than 150 times as great in the region of 
maximum absorption. These data indicate that the reaction 
product formed in the feulgen reaction is not the original 
basic fuchsin. 


6. The Cooperative Study in General Education. GwiLyM 
E. Owen, Antioch College-——The paper will describe how 
approximately twenty colleges and universities, mainly in 


PHYSICAL 


SOCIETY 


the middle west, but also including the University of 
Denver and Mills College, California, have organized to 
study their common problems. In particular, the paper will 
describe the cooperative activities of the physical scientists 
of these institutions in their search for solutions for some 
common educational problems. 


7. The Chemistry and Physical Chemistry of Glass; the 
Physics Problems and Properties Involved in Its Manu- 
facture. ARTHUR S. Watts, The Ohio State University. 


8. The Microscopic Approach to the Knowledge of the 
Properties of Glass. W. J. McCauGuey, The Ohio State 
University. 


9. Some Aspects of the Physical Treatment of Glass. 
Henry H. Biau, The Federal Glass Company, Columbus, 
Ohio. 


10. Viscosity as a Governing Property in Glass Manu- 
facture. H. R. LILLIE, Physical Laboratory, Corning Glass 
Works, Corning, New York. 


11. Some Physical Procedures Used in Ceramic Re- 
search. G. A. BOLE, Engineering Experiment Station, The 
Ohio State University, Columbus, Ohio. 


12. Colorimetry in the Service of Ceramic Technology. 
IsAy BALINKIN, The University of Cincinnati. 
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